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SOME EXTENSIONS OF ‘RAYLEIGH’S PRINCIPLE’ 
By R. V. SOUTHWELL (Cambridge) 
[Received 4 December 1952] 
SUMMARY 

Retaining Lord Rayleigh’s postulate, that any permitted displacement can be 
expressed in a series of ‘normal modes’, this paper waives his restriction to systems 
of which the eigenvalues all have one sign. A simple modification of his approximate 
expression for an eigenvalue allows the wider class of problems to be treated similarly ; 
and used in conjunction with ‘intensification’ (the term used here for the iterative 
process due to Schwarz, Vianello, and/or Stodola) it permits an imposition of double 


limits, also yields estimates of other (unwanted) components in an approximation 
to the gravest mode. 

This last feature will have value—in problems which cannot be solved exactly— 
if two or more of the gravest eigenvalues are nearly or exactly equal. Brief reference 
is made to an example of that kind, relating to the elastic stability of a flat-plate 


cantilever. It will be treated in a separate paper concerned with approximate 


computation of modes of ‘buckling’ (collapse). 


Introduction 

LorpD RAYLEIGH’s presentation (1) of his well-known ‘principle’ relates to 
problems in which the eigenvalues (squares of the natural frequencies) are 
all positive; and strictly speaking assumes the freedom to be finite, since 
it postulates the validity of expansion in terms of the eigenfunctions 
(‘normal modes’). This paper, waiving the restriction to one-signed eigen- 
values, retains his postulate regarding such expansion. Rigorous proof, as 
entailing particular scrutiny of each equation, is not compatible with its 
purpose, which is generalization. Moreover, in effect the freedom is made 
finite in numerical treatment (e.g. relaxation), when derivatives are replaced 
by finite differences. 

The particular class of problems which the paper contemplates is con- 
cerned with the stability of flat plates sustaining edge-thrusts: in it, the 
eigenvalues define the ‘critical’ intensities of that loading, the normal 
modes are the modes of flexural displacement (w) which the critical loadings 
can maintain. When small, w is governed by an equation of the type 


V(w) = AT (w), 


in which V and 7’ are symmetric operators, not necessarily positive. One 
of the worked examples is a problem of this kind. 


But the methods have much wider application, and three extensions of 
Rayleigh’s method are believed to be new, viz. 
(i) ‘Vaisey’s theorem’ (section 4), 
(ii) an alternative (ug) to ‘Rayleigh’s estimate’ Ap, obtained by a use 
of one ‘intensification’ (section 7), 


[Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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(iii) two forms of ‘restrictive inequality’ (sections 11 and 13) whereby 
double limits are set to the numerically smallest eigenvalue (A,) or 
to another eigenvalue nearly equal to A). 
The double limits have value as putting definite bounds to the error of a 
‘Rayleigh estimate’: (1) because when the error is acceptably small com- 
putation can be stopped, and (2) because, when A, has been determined 
closely, more can be learned about the associated mode w,. 

1. For a system defined by a finite number of coordinates, the equations 
which govern a normal mode associated with an eigenvalue A have forms 
typified by éV/ea,—ACT /2a, = 0 (k = 1, 2...., N), (1) 
in which a, is a representative coordinate and V and T are homogeneous 
quadratic functions of the a’s. For a continuous system (N = oo) the 
N equations (1) are replaced by a differential equation of the form 

dw—Ad’w = 0, (2) 
which holds at every point. In this paper # and # will denote any self- 
adjoint linear operators such that 

I (say) = ... | | w dw dedy..., I'/A (say) =... | | ww dedy... 
are homogeneous quadratic functions of the ‘displacement’ w, correspond- 
ing with V and T, respectively, in (1); and whereas in Lord Rayleigh’s 
problems both V and T were essentially positive, here only I is required to 
have that property, the other integral I'/\ may take either sign. In consequence 
the eigenvalues can have both signs, but must still (ef. section 2) be real: 
those which are positive will be denoted by 


those which are negative by 
both series ascending in order of numerical magnitude. (The significance of A 
will appear in section 5.) 

With (2) are associated ‘boundary conditions’ which express the operation 
of restraints. The number of the independent variables z, y...., is not 
restricted, but in practical applications is normally one or two, and is 
taken as two in this presentation of fundamental theory. For brevity the 
same symbols will be employed both for dw, %'w, { { w dw dady and for their 
approximations in terms of finite differences. 

2. Since (2), being linear, does not define completely the amplitudes of 
the normal modes (eigenfunctions), we ‘normalize’ them by the additional 
requirement 


( [ w), Pw, dady = 1 = r, [ [ w;, Ow, dady (3) 
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(which implies no additional limitation, and is consistent with the postulate 
that J is positive). In virtue of the conditions of restraint, any two of them 
satisfy ‘conjugate relations’ of the types 
( ( w, dw, dxdy = 0 = ( [ w, dw, dady | 
al ac # i,). (4) 
|| w,3’w, dady = 0 | | w, O'w, dady | 
These relations serve to establish the reality of the eigenvalues. For the 
existence of a complex eigenvalue A-+7B, associated with a complex mode 
+i W’ (say), would imply the existence of an eigenvalue A —iB associated 
with a mode W—iW’; and when these modes are identified with w, and w,, 
unless B = 0 (so that A, = A,) we have, according to the first of (4), 


({ (WoOW+W’ OW’) dxdy = 0, 


therefore V} 0) = W’ (since J is essentially positive), whether all of the 


eigenvalues have one sign or not, 


Rayleigh’s principle 

3. The postulated validity of expansion in normal modes (section 1) 
requires that, when the A’s have appropriate values, the series 

w= A,w,+A,wet...+A,u,+... 
> (A;.w;,,), Say,T (5) 

will represent any ‘mode of displacement’ (w) which does not violate the 
conditions of restraint, and may be differentiated. From (5), in virtue of 
(3) and (4), it follows that 


I | | w dw dady = > (AP) | 
7 , (6) 
l'=A || ww dady AY (AB/A,)| 
consequently A, as defined by 
Ap = AL/I' = > (Ap)/>d (AB/A;,) (7) 
has a stationary value when the mode w is ‘normal’ so that all of the A’s 
are zero except A, (say), because then Ap, = A,, so 
b(1'/A) €Ap/CA, = A,(1—AT/A, I’) = A,(1—Ap/A,) 
0 for every k. 
This is ‘Rayleigh’s principle’, which (cf. section 1) he propounded as 
holding for systems whose eigenvalues are all positive. It holds also when, 
as here, the eigenvalues have both signs but are all real. 


\ normal mode will be termed ‘positive’ or ‘negative’ in conformity with the sign 


of its associated A. > will denote a summation extending only to positive modes, typified 


by w,; } a summation extending only to negative modes, typified by w_ ; > a summation 


extending to every normal mode. 
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Vaisey’s theorem 

4. Lord Rayleigh (1) argued, further, that A, as defined in (7) always has 
a value intermediate between the least and the greatest of them (or above 
the least, when the system has infinite freedom), consequently Ap as 
computed for a guessed form w (here termed the ‘Rayleigh estimate’) 
may be regarded as an upper limit to A,. This theorem fails when /’)\ 
can take either sign or vanish, for in the last event A,;, has an infinite value, 
But it is easy to establish a theorem—due to Miss G. Vaisey—that A » cannot 
lie between the least positive and least negative eigenvalue (A, and d_,). 

[f A, and A_, typify the positive and the negative eigenvalues (cf. section 1), 
forms equivalent to (6) are 


I = § (Aj) + > (A?) = P+Q, say. 


K l 
I'/A >, (Ak r,) 4 2 (A? /A_;) = R-S, say, 
and as defined in (7) 
Ap = Al/I' = (P+Q)/(R—S), 
where P, Q, R, S are essentially positive quantities. Lord Rayleigh’s 
argument can be employed to prove that 
P>,R, Q > j1_,|S = —2_,S, 
and it follows that 
(R—S)?(Ap—Aj)(Ap—A_1) 
{(P—A, R)+Q+A, S}{P+ |A_,|R+Q—|A_,|S} 
> 0 
This establishes the theorem. 


‘Intensification’ 

5. The practical importance of Rayleigh’s principle lies in his use of it 
(section 4) to obtain an upper limit to A, by simple computations based on 
a guess regarding the form of w,. In most problems of the sort which he 
considered the form can be guessed with considerable accuracy, and then 
in virtue of the stationary property (section 3) the resulting estimate of A, 
is extremely close. But in the wider class considered here little confidence 
can be put in a guess regarding w,: then, it may be improved by a process, 
here termed ‘intensification’, which is variously attributed to Schwarz (2), 
Vianello (3), and Stodola (4). This derives from w, successively, other forms 
Wy, Wyss etc., satisfying the edge-conditions and 

du, = Ak'w \ 
Fwy, = AP'w, }, 
, etc., 


in which A is any convenient multiplier. 


(9) 
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Provided that w, is a component of the guessed form w, when every eigen- 


value is positive, then wy, wyy...., etc., will approximate more and more 
@YVS Nas , ° ° ° 
al closely to w,; when the eigenvalues have both signs they will approximate 
r above i ‘ 5 "he , . 

to w, or w_, according as A, or |A_,| is the smaller. For if w is expressible 

“ in the series - 
imate’) | , w= > (A, w,), (5) bis 
. then i - : 
n I'/A Wy AS (A; w,/A,) 
> value : : > ) 2 
Wry \? YS (A; w,/Az) }> (10) 


» Cannot 


, etc. 


tion ]) | 80 the gravest mode becomes successively more predominant in wy, 


‘Intensification’ used in conjunction with Rayleigh’s principle 
6. Temple (5) and Temple and Bickley (6) have used intensification in 
8) | conjunction with Rayleigh’s principle; and Koch (7) did the same thing 
in effect, by propounding as an upper limit to A, 
\ | | w dw dxdy/ || w dw, dxdy = || w dw decdy/ | w 3'w dady 


} 


elgn s 


in the notation of this paper), for the ratio on the right is Ap. By this 
procedure, when the eigenvalues are all positive, the accuracy of Rayleigh’s 
estimate isimproved: for, as derived from w, it is given by (7), and as derived 
from wy, it is 

(Az); = Al’ /I”, say, 
where r 2 [| wy dw, dudy = A? > (AZ/dZ) (11) 


; 


1” = A3 | [ wy dw, dady = A® ¥ (ABR) 
se of it (Ande Ap = & (ARAB) > (AZ/A,)/LS (AR/AB) S (AR} 
sed 01 | — } [ A? A2(A> AS) 1+A;4)]/> (AR/AR) S (AR) (12) 


d ther (S extending to all pairs such as r and s) 


| when every A is positive (otherwise the co-factor of some 


A? A®, in (12), will be negative). 
rocess 


rz (2), An alternative to the ‘Rayleigh estimate’ 
forms 7. All of Rayleigh’s arguments apply to 

uz (say) \2 0/1" = > (AR)/> (AR/23). (13) 
whether all of the eigenvalues have one sign or not, because every A? > 0. 


( That is to say 9 >) ‘ 
si uz > the smaller of A? and A? ,, 


also, as in (12) 


(u&); < we, when (u2), is derived from wy. 
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Finally, by (7) and (13), 
mda = PITT" = (5 ABATE HM). & (ab 

[ A? A2(A>4—A>1)?]/¥ (Az/AR). > (AR) (14) 


r,s 
< 1, whatever the sign of the eigenvalues. 
























Thus for every system intensification should invariably reduce |ug|, an 
estimate which imposes a closer limit than Ap, on the magnitude of the 
numerically smallest eigenvalue. Consequently a fatlure of \us| to decrease is 
unambiguous evidence of a computational error. Vaiscy’s theorem (section 4) 
is not required. 

A price in labour is paid for these advantages, since w must be intensified 
(to obtain w;) before a value can be attached to J”; but it is justified by the 
resulting information, in the complex problems which this paper contem- 
plates. The indication of error is especially valuable, since it shows when 
the attainable accuracy is limited: in theory, intensification cannot fail 
(by reducing unwanted components) to improve any guess to which it is 
applied; but an argument which neglects inescapable errors cannot yield 
a trustworthy assessment of its practical merits. Much has been written 
which presumes exact solution to be impracticable, yet proceeds as though 
such errors were avoidable. 


Rayleigh’s principle used in conjunction with two intensifications 
8. Advantage having accrued from one intensification, it is natural to 
inquire whether two or more will be justified by the resulting gain in 
information. Two at least should normally be practicable, and their 
cumulative error should not be serious: i.e. they should yield close 
approximations to w,;, w,;; and values for the integrals 





i | [ wy Fwy, dxdy A® S [ARAL] 


I ( [ w dw dxdy = > [Al 

pet | 

i | | wdw,dady =A > [AP/A,| | 
r" [ [ w Iw, dxdy = A Zz | Az rj: | [ [ wy, wy dxdy >, (15) 

wen ** | 

| 

| 


liv | | Wy Pwy, dady = A* > [AZ/AZ] 


of which J, J’, I” have appeared already, in sections 3 and 6. Taking these 
values as reliable, the remainder of this paper draws conclusions in regard 
to the modes and eigenvalues. 
The argument is based throughout on Rayleigh’s principle, occasionally 
extended in the manner of sections 4, 6, and 7. In sections 18—20 it will 
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be applied to a representative set of computed values: here it is illustrated 

by an example that can be solved exactly—a uniform strut with hinges at 

both ends.+ For that system (w denoting the transverse displacement) 
ow d*w/dx*, ow Ww, A PL? B, 

when B denotes the uniform flexural rigidity, Z the length, and P the 

thrust, and when La measures the distance from one end. The eigen- 


values—all positive—are 


Ay, Ag, Ag,---, Cte. = w*X (1, 2, #...., ete.), (16) 
and (9) take the forms 
d?w,/daz?+ Aw = 0, d?w,;/dx?+- Aw, = 0,..., ete., 
so when the guessed mode is w = a(1—2), 
then 12w, = A(a—2a3+-24), 
360wy, A?(3a— 523+ 325—2%), (17) 
and $f = 1, 30/’ = A, 50401” = 17A2, 
144 x 63017” = 31A3, 1386 x 14400/'¥ = 691A4. 





Therefore in this example, according to (7), 
Ap as deduced from w is AI/I’ = 10, 
Ap as deduced from w, is Al”/I” 17x 18/31 9-870 967 74, 
according to (13), 
uz, as deduced from w is A?J/I” 1680/17 (9-941 002 434)?, 
uz, as deduced from w, is A?/”//'¥ = 34x 1980/691 
(9:870 360 84)?. 


Observing that for this system, according to (16), 


A, = 7 = 9-869 604 401 089 359, (18) 
we see that 
Ap > (Ap); > Ay, as predicted in section 6, 
9 9 9 
and that we > (us), > Ai, 
also uz, r},, (uz), < (Ap)Z, as predicted in section 7. 


(Here, of course, the eigenvalues are widely separated.) 


A lemma 
9. Hereafter we shall occasionally need to know the limits imposed on x 
by the double inequality 


) 


2<2+27! < A%1-+e). (19) 


+ The same equation governs the free vibrations of a tensioned string—a case considered 


in great arithmetical detail by Temple and Bickley ((6), section 5.2). 
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This requires that x > 0, so that #2+1 < 2(1+¢)a,—i.e. 
{e—(1+e)}? < 2e+e?; 
therefore it implies that 
ax lies between 1+-€+.,/(2e+e). (20) 
So also does the inequality 
—2 >4+2-! > 21+), (19)A 


which requires that x < 0; while the inequality 


4< (4+a-1)? < 4(1+€)?, (21) 
which permits 2 to have either sign and requires that 
2 < |x|+ |x|“! < 2(1+), the positive root of 4(1+e)?, 
implies that x| lies between 1+-«+.,/(2e+e?). (22) 


‘Optimal synthesis’ of w and w 
I 
10. For a mode defined by 
Ww, = w+ yu, 


(y denoting an arbitrary constant) the ‘Rayleigh estimate’ is 


(Ap),, Say [ ( w, dw, dxdy/ [ [ w, Fw, dady 


2 A ( [ w, dw, deedy | [ [ w, Hwy +ywy,) dady (23) 
A(I+2yI'+21")/(I'+2yl’ +721"), | 
and the alternative estimate 2 is (similarly) | 
(ug),, say = A®%(I4 2yI’+y71")/(I"+2yI”" +I). 


Therefore, as they depend on y, (Ap), has a stationary value when 


> 
t 


I.I”—I["?+y(1.1"—-I'.1")+y(1' .I”—1") = 0, (24) 
(u§), has a stationary value when 
I.I”—1' I’+y(1. Pv¥—1")4+- yr’. Fev— 1" I”) = 0. (25) 


The lower stationary value of (A,),, will normally (always, if the eigenvalues 
are all positive) be nearer A, than either A, or (A,,);: and the lower stationary 
value of (u&), will always be nearer Aj than either p% or (u2);. When the 
curves representing (23) are flat (because w is a close approximation), their 
minima can be found by inserting trial values of y close to the relevant 
root of (24) or (25). 

An earlier paper (8) thus combined two guessed modes w, and wz 
(by ‘optimal synthesis’) to minimize A, and thereby to obtain a closer 
approximation to w, than either w, or wy. (In effect this is the ‘device 
of the variable parameter’—apparently due to Rayleigh, though some- 
times attributed to Ritz. (Cf. (9), section 526 and p. 500.) The y’s derived 
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from equations like (24) are always real and yield two modes which are 
conjugate: for the modes (w+-y,w,) and (w+y,w,) will be conjugate, 
by (4), if 

I+(yity)l + y2l" 0, ’+(yit+ye)l"+y1 721” = 9, 
and both conditions are satisfied when y,, yz are the roots of (24), because 
then (1'. 1” —I"*)(y, +0, ¥1 2) = (U’.L’—I.1"), (1. 1” —I'*). On the other 
hand, when y,, vy. are the roots of (25), then 

(l).P°—L" LD" \yyt+ye v1 ye) = (1"2-—I. 1), (1.1”—I'. 1’), 


and the first but not the second of (4) is satisfied. 


IMPOSITION OF DOUBLE LIMITS ON A, 
First restrictive inequality 
11. When /, J’, J” have been computed in accordance with (15), a value 
can be found for ¢ as defined by 
‘ ” ’ 22) / 2/)/ 
by (15), where BR = Aj(A?2/AZ+1) > 0, (26) 
Ny. = A/Ay,+A,/A. 
This equation has the form of (7), section 3, BZ and dj, replacing AZ and A, 
in that expression for A,: consequently Vaisey’s theorem (section 4) can 
be adduced to prove that 2(1+-e) cannot lie between the smallest positive 
and smallest negative A,, neither of which (since every A is real) can lie 


between +2.+ That is to say, 
either 1+ 1, and then the smallest positive A’ < 2(1+e), 
or 1+e 1, and then the smallest negative A’ > 2(1+e). 


Suppose, then, that /’ and hence (1+) have been found positive, and 
let a stand for that (positive) value of A,./A which yields the least positive d’. 
Then 2<at+a-1 < 2(1+¢), (19) bis 
and the argument of section 9 goes to show that 

x lies between 1-+-e +,/(2e }-¢?), (20) bis 
The same conclusion holds when J’ and (1+) are negative, x then denoting 


that (negative) value of A,./A which yields the least negative 2’. 


A refinement of the first restrictive inequality 
12. The nearer (1+) is to +1, the narrower are the limits imposed on x 
by (20); and when w is replaced by w, as defined in (23), the value of « is 
given by . 
2(1+e) = {74 1"4+2y(1'4-1")4+ PL" + 1) + 2yl’+y1") (27) 
4+ (AJAz—A;/A)?. 
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instead of by (I+ 1”")/I’. From this equation, as in section 10, two values 

of y can be found for which « takes a stationary value: they are the roots of 
MPI+Ll)—-l'('4+-1")4+y" (14+ 1")— ("+ Fev} 

+I" +1")—1"(1" + Try} = 0, (28) 

and inserted in (27) they yield two values for « of which the smaller may 

be used to impose limits on x. 


Second restrictive inequality 

13. Alternatively, a value can be attached to 

4(1 +e)? = (1+21"+]'%)/1" = > (C2/X2), by (15). 
where Ci = AR(A?2/AZ+1)? > 0, (29) 
and (as before) dy = (A/A,+A,/A). 
This equation has the form of (13), section 7, Az being now replaced by A}; 
so Vaisey’s theorem is no longer needed, because Rayleigh’s original 
argument goes to show that the smallest A’? (which + 4: ef. section 1], 
footnote) < 4(1+.)®. Consequently if x stands for that value of ),/A 
which yields the smallest value of A’*, then (section 9) 
x| lies between 1-+-¢-+,/(2e+e?), (22) bis 

(1+-e) now denoting the positive root of (1+). 

The limits thus imposed on 2 will be closer if, as will normally be the 
fact, « is smaller as given by (29) than as given by (26). They could be 
made still narrower by minimizing (J+ J‘*)/J” in the manner of section 12; 
but the resulting equation would involve JY and J¥', which are not known. 
In this respect (26) has an advantage over (29); and it has another in that 
it gives a definite sign to the A in question. 


A test to verify that the restricted eigenvalue is (or is very nearly 

equal to) A, 

14, When zis shown to havea value close to 1, A must approximate closely 
to some positive eigenvalue: when x has a value close to —1, —A must 
approximate to some negative eigenvalue. Thus in every case some eigen- 
value can be ‘bracketed’ by double limits. By methods which a subsequent 
paper will discuss, it can be made almost certain that the gravest positive 
mode predominates in w, and hence that A approximates to A,: then (1+e) 
may be expected to be positive, and in that event (20) will put double 
limits on A,. 

Evidence in support of this expectation is obtainable. Let A, denote 
the eigenvalue to which A is closest, so that 2 = A,/A lies close to 1: then 
in the mode (w’, say) defined by 


w’ = w—axw, = > [ A, w,(1—A,/A,)], by (10), 
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the component w, is nearly eliminated, and the contribution of another 


8 


mode A,. (say) will be 


also nearly eliminated if A, = A 


very little altered if |A,./A,| is large, but 

greatly increased if |A,/A,| is small (and especially if A,/A, < 0). 
[f, then, any |A| exists which is appreciably less than |A,|, the corresponding 
mode will make a larger (relative) contribution to w’ than to w, and as such 
will yield a smaller value for u.2,; consequently, if 12, has a greater value for w’ 
than for w, then no |A| appreciably less than |A,| can exist. (The relevant 
formula is the last of (23) with y replaced by —z as here defined.) 
Exemplification by the uniform strut 

15. In illustration of sections 11-14 we now apply their methods to the 

uniform strut (section 8), for which the J’s are given exactly by (18). We 
suppose that A has been identified with Ap (= 10), so that I' = I exactly. 
On that understanding, 


J=['’=1/3, JI” = 170/504, J” = 775/2268, J'¥ = 17275/49896. 
(30) 
These values inserted in (26) give 
€ 1/168 = 0-005 952 381, (2+?) = 0-109 271 189, 
and so, according to (20), for some A, 
x = A,/10 lies between 1-005 952+-0-109 271. (i) 
Inserted in (29) they give 
(1+)? = 14+247/67320 = (1-0018,)?, 
and so, according to (22), for some A. 
a A,.|/10 lies between 1-0018,+-0-0605,. (ii) 
The limits (ii) are narrower than (i), but do not define the sign of the A,. 
On the other hard, (i) shows that some positive A,./10 has a value close to 
unity, so presumably (ii) defines a positive eigenvalue.t Moreover, the 


limits (i) can be nerrowed in the manner of section 12. 


Thus (27), when the values (30) are inserted, becomes 


€ cia(1 t 2 (9A 5y* 198) (1 + 85y 434 7 " 756), 


} 
t It is known, of course, in this instance that the eigenvalues are all positive ; but that 
knowledge is not used here, because in this respect the problem is not representative. 
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giving « a minimum value 0-000 091 565 898 when y = — 10-305; and (20), 
when this minimum value is given to «, requires that for some ,. 
x = A,/10 lies between 1-000 09,+.0-013 53,. (iii) 


These limits are finer than either (i) or (ii), and the sign of A, is defined. 


16. Now the test of section 14 can be applied: the criterion being 2, as 

given by 

ps = 18°(1 —2a+ 85a2/84)/(1—3102/153+ 345522/3366) 
with 2 restricted in accordance with (iii). This (since x == 1) has a value 
greater than (2), as found in section 8 (viz. 1680/17) in a ratio sensitive 
to x but not less than 13: hence it can be taken as almost certain that ), 
is the restricted eigenvalue (for here we have other reasons for believing 
that A, and A, are not nearly equal). 

On that understanding the limits can be made still closer, because 
(cf. section 7) any pg > A,, and by ‘optimal synthesis’, in the manner of 
section 10, a minimum value 9-869 680 1 is obtained for (14s), as given by 
the last of (23). Taking this as the upper limit, and combining it with 
the lower limit given by (iii), we have finally, 

9°86558 < A, < 9-86968 (31) 
as the restriction imposed by our methods on A, for the uniform strut. 

By adoption of the mean value of the two limits in (31), the correct 
value (7?) is reproduced with an error of about 0-02 per cent. And th 
methods by which the limits were deduced nowhere require an exact solution 
of the problem. 


An example in which the eigenvalues have both signs 

17. This paper had its origin in a recent attempt to improve the accuracy 
of a solution which L. Fox and J. R. Green obtained by relaxation methods 
in 1941 (reference 8). A full description of the problem will be given 
elsewhere: it concerned the elastic stability of a flat-plate cantilever under 
loads which, since they induce both positive and compressive principal 
stress, could also entail instability if their sense were reversed. This means 
that the eigenvalues (critical load-intensities) have both signs,+ so that 
{ayleigh’s method (in its original form) is not applicable: moreover the 
gravest mode is characterized by nodal lines, so cannot be guessed with 
accuracy as it can in simpler problems. Fox and Green, by ‘point relaxa- 
tion’ combined with ‘optimal synthesis’ (section 10), had obtained an 
estimate of A, which they believed to be ‘correct within about 1 per cent.’; 
but had attempted less accuracy in their determination of w,, and had not 


{ This consequence was first pointed out by Dr. D. N. de G. Allen. 
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established that their accepted mode and eigenvalue are in fact the 
gravest. 

Miss Vaisey, in the more recent investigation (not yet published) em- 
ployed, along with methods specially devised, the procedure which has 
been here termed ‘intensification’; but arrived repeatedly (on her finest 
‘net’) at a point where closer approximation seemed impracticable unless 
further indications could be drawn from theory. This difficulty had not 
been resolved in June 1952. She had obtained a mode believed to approxi- 
mate closely to w,, and this mode (w) she had twice intensified to obtain w; 
and w;; but these, though they yielded confirmatory estimates of A,, 
both differed sensibly from the mode obtained by Fox and Green. 

She had deduced from them the integrals defined in (15), identifying A 
with A, as deduced from w and so making J and I’ exactly equal (as in 
section 15). Her other computed values gave 

[”/I = 1-0012,, I” /I = 1-0019,, Tiv/T = 1-0029, (32) 
(a ‘dropped’ figure such as 3 denoting some value between 2 and 4 but 


possibly nearer to 2 or to 4 than to 3). 


18. These values inserted in (7) and (13) give 


Ap/ A as deduced from w Rif’ l, 


Ar \ 9 e9 Wy ag 0-999., 
(p45/A)* o " wow = f/f’ (0-999,)?, 
(us/AP?  ,, » ww, =I'"/I'v = (0-999,)*. 


They thus confirm the prediction of section 6, that Ap > (Ag);, also the 
predictions of section 7, that u2 > (4%), and that 
wy < (Ap), (uS < Ari 
(This is partial verification of their accuracy.) 
Inserted in (23) they give 
(Ap )y/ N= {Lb p24 12+, x 10-4 }/{(1 ey)? (24 yy +1959) x 10-4, (i) 
(u2),,/ A? f(1+-y)?+12-, x 10-4y?}/{(1-4 y)?+ 
| (12+, +38-¢y+29-gy2) x 10-4}, (ii) 
whence (24) becomes 12-,+7y—4-9y? = 0 (iii) 
and (25) becomes To +5 yy- -2+y" 0. (iv) 
According to (i) and (iii), 
(A,),/A has a maximum value — 2-464 15039 when y = —0-99914 and 
a minimum value 0-999 140032, when y = 2-47. 
According to (ii) and (iv), 
(u%),, A? has a maximum value (2-081 307 18)? when y = —0-99907 and 
a minimum value (0-999 118 09)? when y = 3-40. (33) 














270 R. V. SOUTHWELL 


Inserted in (26) the values (32) give 
2(1+e) = 14J"/I = 2-0012;: 
therefore, according to (20), 
x lies between 1-0006,+.0-0347, (34) 
when x denotes that value of A,./A which yields the smallest positive («-+-2-), 
The alternative formula (29) of section 13 requires, when the values (32) 
are inserted, that 
(1+-e)? = 4-0053,/4-0048, = 1-000 12,: 
therefore, according to (22), 
x| lies between 1-0000,+0-0112.. (35) 


These are much closer limits; but they do not define the sign of 2. 


19. Before proceeding to the refinement of section 12 we apply the test 
described in section 14. Its criterion here is the value of 42, as given by 
pe/A? {(1—a)? x 104+-12-, x {(1—ax)? x 104+ 12+, —38-, 7+ 29-, 2?) 
when the range of x is restricted by (34). This value + 1-568; so it is A,, 
or some very nearly equal eigenvalue, to which (34) relates. 
Now introducing the refinement, we insert the values (32) in (27) to 
obtain 
2e (12:, a-sy-4 2+3y") {(1 + y)? x 104+ 24-5 of. 19°, y?}, 
according to which 
2e x 10* takes a minimum value 0-399 023 1488 when y = 2-208 
(and a maximum value — 50 330-1234 when y 0-991). 
When the minimum value is inserted in (20), that inequality requires some 
value of A,./A (= 2) to lie between the limits 
1-000 019 4, + 0-005 217. (36) 
If, then, the A in question is A,, which cannot exceed the minimum (sg), as 
given in (33), we have, finally, 

0-993 80, < A,/A < 0-999 118, (37) 
the upper limit coming from (33), the lower from (36). The limits, now, 
are very close indeed: identification of A, with their mean cannot entail 
an error greater than 0-27 per cent. 

Miss Vaisey’s value of A was 6-002, so the final limits derived from her 
computations are Ron — ..a0n thn PY 
ie : 5-964 805 < A, < 5-996 706. (38) 

If the A in question differs slightly from ),, it lies between 


A(1-000 019,,-+-0-006 217), 
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by (36); i.e. in the range 
5-964 805 < A < 6-039 434. (39) 


Fox and Green ((8), section 32) gave the figure 5-949; so this paper con- 
firms their belief (cf. section 17) that their estimate was correct within 1 per 
cent. (According to (38) it is in error by 0-8 per cent. at most; and according 
to (39), if it is low by more than 1 per cent., then A, and A, are equal to 
within 4 per cent. It may (and of course in theory should) be too high: 
in that event A,, A, will have somewhat wider separation.) 


20. When close limits have been thus imposed on A,, more can be learnt 


about the associated mode w, and about other modes which contribute 


l 
sensibly to w. This aspect will be discussed in a subsequent paper. 

In theory more could be learnt from three or more intensifications and 
from the resulting integrals of type (15); but in practice this further informa- 
tion will be less trustworthy since the cumulative errors will be larger. 
It would seem better, when the optimal approximation to w, has been 
found in the manner of section 10, to start afresh with this as the starting 


assumption. 
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MODES OF VIBRATION OF A SUSPENDED CHAIN?+?+ 
By DAVID 8. SAXON and A. 8. CAHN 


(U.S. National Bureau of Standards, Los Angeles, California) 
[Received 14 February 1952] 
SUMMARY 

In the following paper a method is described for calculating the characteristic 
requencies of a suspended inextensible chain vibrating with small amplitude in the 
lane of the catenary forming the equilibrium configuration. An asymptotic solution 
f the linearized equations of motion is obtained such that the accuracy of the results 
.creases as the mode number increases and/or, as the catenary becomes flatter. 
1. Equations of motion 
We consider an inextensible chain of length L suspended from two points 
it the same level at a distance a (< L) apart. We measure the arc length, s, 
toa given point on the chain from the mid-point of the chain, and denote 








the space coordinates of this point by 2(s,t), y(s,t), where x and y are 
referred to the lowest point of the equilibrium catenary as origin, as 
indicated in Fig. 1. We suppose that the end-points are located at x = +a/2; 


y=b. The equations of motion can now be written as follows (cf. Rohrs (1) 


and Pugsley (2)): a2) ‘ on 
: (7 ‘ ). (1) 
ct Cos cs 


c Cs 


Here 7'(s,¢) is the tension per unit mass per unit length in the string. 


The preparation of this paper was sponsored (in part) by the U.S. Office of Naval 
Research. 
t The authors’ interest in this hundred-year-old problem resulted from a conversation 
th Dr. E. U. Condon, then Director of the National Bureau of Standards. Dr. Condon 
id seen a reference to it in an article entitled ‘Some experiments in dynamics, chiefly on 


brations’ by John Satterly in the American Journal of Physics, 18 (1950), 405. Professor 
Satterly stated that he did not think that the calculation of the period of a swinging catenary 


swaying in its own plane had ever been done. 


[Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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In addition to the equations of motion the chain is subject to an equation 


of constraint ax\2  [ay\? 
(; ) a ( ; ) os . 8 
cs cs 


The equations describing the equilibrium catenary are obtained by 
solving equations (1) to (3) with the time derivatives set equal to zero, 
These static equations are easily solved and, expressed for convenience in 
terms of the angle « which the tangent to the catenary at s makes with the 
x-axis, yield the following well-known equations 


, sinh-!(tan «) 
a 2 tan ap 
L 

y= - (sec a— 1) 

2 tan a» 

, (4) 

| : gL l 

2 tan wp» COS a 

L tana 
s = —- 

2 tan a 





where a, is the angle at the end-point (a/2,b). 


2. Small amplitude approximation 

We obtain a form of the equations of motion for the case of small ampli- 
tude vibrations by assuming a solution differing little from the equilibrium 
catenary. Thus we write 


x = x,(s)+&(s, t) | 
y = y,(s)+n(s,t) +, (9) 
T = Tis)-+7(8,1) 





Cees 


where the subscript e refers to the equilibrium values given in (4). Substi- 
tuting expressions (5) into (1), (2), and (3) and neglecting quadratic terms 
in the small quantities €, , 7, we obtain 


oe _- (7, eg dx, . 
ot* Os\ ~ és ds 
o*n C ay OR dy . 
- = 7 4 = T (6) 
08\ “0a dé 


dx, c& | dy, &n 


— 0 
ds 0s ds és 
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of equations (4); we also assume €, », 7 are proportional to exp(—iwt). In 


this case equat ions (6) become: 





1 
NE cos2a “ | cos — S + (cos x)T]}, (7) 
ax da 
A*n cosa COS a Dy (sin | (8) 
aa q dx 
dé . l 
COS a + SIN ‘: 7 0 (9) 
dx (LX 
Lw* 
where A- = (10) 


29 tan Xo 
This system of ordinary differential equations contains the three unknown 
functions €, n, 7 and is to be solved subject to the boundary conditions 
f, 7 0, a t Xp. 

\s an indication that the problem is completely defined, it can be shown 
that elimination of 7 and » from (7), (8), and (9) leads to a third-order 
equation in €, the solution of which contains three arbitrary constants 
see Appendix I). Determination of 7, once € is known, is obtained by 
integrating (9), which adds a fourth constant to the solution. These four 
constants are just sufficient to satisfy the four boundary conditions, pro- 
vided that A satisfies the appropriate secular equation. 


An important fact should be pointed out, namely, 


constant C,” 
(11) 
& Q, 
NMC 
with 7 : 
COS a 


is an exact solution of equations (7), (8), and (9). This may be verified by 
substitution. Note, however, that 7 = 0, € = constant is not a solution, 
although one might have been tempted to assume this to be so. This 
expresses the fact, which is not immediately obvious, that a hanging chain 
can be moved a small distance vertically without disturbing it from the 
equilibrium configuration, but that it cannot in like manner be moved 
sideways. Use of solution (11) will be made later. 

For a better understanding of the complicated equations with which we 
we dealing, it is instructive to examine the limit of flat equilibrium 
catenaries, i.e. when |a * 1. In this limit, it can be shown that the 
solution is approximately (see Appendix I) 

y= sm (Aa), 
s 


co 


‘ 
=e, 
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The requirement that 7 vanish for « = +a, then gives 
n \7 
A= , : (13) 
e-T 3/ % 


Although this result is only a first approximation, the essential features 
turn out to be correct. The major error is that actually no solution exists 
corresponding to A = 7/2) and the true lowest mode corresponds to 


)) 








A = z/a». Otherwise these rough results turn out to be even quantitatively 
acceptable for ap» 1 (see equations (33) and (34)). 

Other important results follow investigation of this special case. From 
equation (13) it follows that, for small values of a»), A > 1. Further, since 
A is roughly proportional to the mode number, for high enough mode 
numbers, A > 1 for any ap». Also, for small a», the waves are essentially 
transverse since the horizontal component of displacement € is small com- 
pared with the vertical component . However, equation (9) suggests that 
the transversal character is true as long as A is large. For high modes of 
a deep catenary, for example, € <7 near the bottom where a is small, 
while on the steep sides the opposite is true. 

These characteristics make it seem reasonable to seek asymptotic solu- 
tions for A > 1 since such solutions will be applicable over at least some 
interesting range of the configurations and mode numbers. Likewise it 
seems desirable to take advantage of the nearly transverse character of the 
waves by introducing new variables u and v, respectively, locally perpen- 
dicular and tangential to the equilibrium catenary, as indicated in Fig. 2. 
Thus, we substitute 


u —€sin «+7 Cos a, 


(14) 


v €cos «+7 sin a, 
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into equations (7), (8), and (9) to obtain: 


A? , - 
—v vCOSa-—-U COS a- 7’, (15) 
COS" 
A? . — . . 
— U wu" cosSa+wu’ sin a—uUcosa—T+vsiN a, (16) 
COS“ 
vy’ U. (17) 


The primes indicate differentiation with respect to «a. Elimination of u and 


7 from this system yields 


. A ” 
(cos a)u" (2sin a)v (cos x + v 


| navies 1) (2sim ap" a v=0Q0. (18) 
Boy 


cos’ cos? 
Equation (18) is to be solved subject to the boundary conditions: v = v’ = 0, 


when « +- Xp. 


3, Asymptotic solutions 

We next seek the four independent solutions of (18), at least asymptoti- 
cally for large A. The procedure which will be used is roughly equivalent 
to separating the fourth-order equation into two second-order equations, 
one of which leads to a pair of functions varying rapidly compared to the 
slowly varying trigonometric coefficients, the other leading to a pair of 
solutions which themselves vary slowly. That such a scheme is a suitable 
one can be seen quite directly from the following physical considerations. 
For large A (or high mode numbers) we expect the problem to reduce 
substantially to one of transverse waves on a string with slowly varying 
tension. The asymptotic (WKB) solution for such a problem can be written 
down almost by inspection, of course, since it involves only a second-order 
differential equation for the transverse amplitude. Consequently, for large 
A, we expect that it is possible to separate out of the rigorous equations 
which govern the motion a second-order equation leading to such a 
solution as a first approximation, and this indeed turns out to be the case. 
Use of these rapidly varying solutions iteratively furnishes higher order 
corrections which essentially take into account such factors as the curvature 
in the chain and the dynamic changes in the tension. 

In addition to these corrections, and more important, the rigorous dif- 
ferential equation must, of course, have two other solutions of a different 
character which physically have their origin in the fact that the waves are 
not exactly transverse and the small, but not negligible, tangential displace- 
ments must also satisfy the boundary conditions. These then are the two 
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slowly varying solutions which are rather like two integration ‘constants’ 
required to fit the two additional boundary conditions. These slowly 
varying solutions cannot be found exactly, but they can be obtained 
asymptotically to any required degree of accuracy using iterative methods 
as we shall see. 

We now follow the procedure indicated above, considering first the 
slowly varying solutions since they are simpler to treat. Because these 
functions presumably vary only trigonometrically with «, their derivatives 
are of the same order in A as the functions themselves, and hence equation 
(18) is seen by inspection to consist simply of terms of order A? and terms 
of order unity. As a first approximation, we retain only the terms of order 
A* so that the differential equation reduces to the second-order equation 


” 


vy’ +(2tan a«)v’—v 0. (19) 
However, we already know one exact solution of (18), namely 
v = Csina, u = Coosa (20) 


corresponding to the solution € = 0, » = C [see (11)]. This is, of course, 
a slowly varying solution, independent of A; and hence satisfies (19), as well 
as (18), exactly. The construction of the second solution of (19) is thus 
trivial and we find at once 


v = D(cosa+ asin «), uU Da cos x. (21) 


This then is the first approximation to the other slowly varying solution. 
To improve this approximation, one can establish a simple iterative 
procedure as follows. Denote by vy the solution (21) and expand v in the 


series o <a Se. 
— 


n 
Substituting this series into (18) and equating coefficients of like powers 
of A to zero, we obtain the system of equations 


” 


vy, + 2(tan «)v,—v 


n 


” 


— CQS3ySaiv Mts hid Lia Dts ’ \ 
= —COS*aj,U__ 1 — 2(tan a), _1+U,_1— 2(tan a)v,_} 


t; 
so that each v, can be determined in order from the preceding one. This 
procedure is easily carried out in principle since both solutions of the 
homogeneous equation, and hence the Green’s function, are known. Thus, 
for example, we find in this way, 


sin4 
v, = —2(1 ies : 


Actually, since the slowly varying solutions are essentially correction terms, 
we shall not require more than the first approximation (21). 
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Next we obtain asymptotic solutions for the two rapidly varying func- 
tions. In order to indicate the nature of the results we first show how to 
do this in a very rough way directly from equation (18). Since the rapidly 
varying functions are expected to be of the form exp[+iAf(a)], the nth 
derivative of such a function is to be regarded as of order A” times the 
function. Equation (18) then contains terms of order A* through A and 
retaining only the fourth-order terms, it reduces to 
ae 


— WV 0. 
cos? X 


io exp! =,\ | S| 
J costa 


and thus establishes the self-consistency of the assumptions. 


which yields at once, 


Using this kind of argument, we could then analyse (18) more carefully 
to get more accurate results. Actually, however, it seems to be easier to 


accomplish this using equations (15), (16), and (17). The main reason is 


vm ol, u), 
A 


and thus v is not required as precisely as u. This can be regarded as an 


simply that, as seen from (17), 


explicit statement of the fact that the waves are nearly transverse. Using 
the above relation, we see from (15) that 

t= O(u) 
and hence, dropping terms of order u from (16), we obtain as the first 
approximate equation for u 


9 


wu” —w' tana+——_ u = 0 (22) 

cos*a 
and, as expected, this is precisely the equation for transverse vibrations on 
a string with tension 7' proportional to 1/cosa. To order 1/A, the solution 


of (22) is easily found to be ‘ 


so Et Oe ; 
u ~ costaexp| +7A | : ' (23) 

J costa 

0 

We now use this result iteratively to obtain a solution to order 1/A?. We 
can immediately get v to the required order using (17); thus 
] . 
utOj—u). (24) 
re 

Using (15), we can now express 7 in terms of u and after a little manipulation 
we find 


COS? a 
~ 


T 1A 


] 
7 =~ 2(cosa)ju+Ol[—wu}). 
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Substitution of these relations for 7 and v into (16) then gives 


wu” —(tan «)w’ + a +340 . w= 0. 
cos’a A 


u exp( 4 iA | fo(a) da — { F(a) datija | fo(ax) da ~} 


0 0 0 


Writing 


#eR 


we easily find the adiabatic solutions 


| 
Cu costa ettAs+a/A) 4 o( =U}, (25) 
A? 
. ‘ Costa yr py l ‘ 
and from (24) _ — ettAs+9 4 Of — w), (26) 
t iA A? 
a 
' ra 
where f (x) | sa (27) 
; J Costa 
0 
: 
and q(x) L | (14-4, tan*a)cos!a da. (28) 


Using (20), (21), (25), and (26), we can then write the complete solutions, 
to order 1/A?, as 


u = A costae’AS+9% + Bcostae-iAS+9% + C cos a+ Da cos x, 
Z ee ae eee ie ate ; 
v — cosia et4S+9/4)_ —__ eggig e-*A+9/4N + C sin a+ D(cos «a+ asin «). 
u/ U/ 
nde e rec 4 and v vanish at «a ta», the solutions 
Under the requirement that uw and nish at »: the solution 


decompose into odd and even pairs which lead immediately to the following 
transcendental equations for A: 


odd modes (u odd about « 0, v even) 


l x COs! l 
tan(Af-+-q/A) eee he (29) 
+9 ~ \ ecosa+asina A 
even modes (u even about a 0, v odd) 
] cosia ] - 
cot(Af+g/A) : = —k(a). (30) 
I+9 A sin « A 


For convenience, the subscript 0 has been dropped from ap. 
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By observing that forA > 1, we obtainAf ~ nz and Af ~ (n+-4)z respec- 
tively, we can easily obtain explicit solutions correct to order 1/A? by 


expanding. In this way we finally find 
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We have thus obtained comparatively simple explicit expressions for 4, 
and hence, by (10), for the characteristic frequencies of the chain. The 
functions f, g, h, k are tabulated in Table I. The quantities A,,/z for the first 
five modes, as determined from these functions, are tabulated in Table II 
and are plotted against « in Fig. 3. Included in this figure are the following 
experimental results: 

(1) Data for the first five modes as determined with the kind assistance 

of Professors I. Rudnick and R. Leonard of the Department of 
Physics, University of California, Los Angeles. These points were 
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TABLE I 


g(a) 


1 2 } 
44 tan“ax)cos?a¢ 
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*992C 


235 
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———__————— 


TABLE II 








As can be seen from the figure, agreement is 
lowest mode for angles as large as 70°. 
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(3) Data for the first three modes as determined by Pugsley (2). 
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for « <1, f(a) ~ « and since g, A are of order unity, for even modes we 
obtain = 

A~— (a <1), (33) 
x 


in agreement with the previous simple result. For odd modes, on the other 


hand. since k(a) > 1/a as a 0. we obtain 


(n 4) a ] 


a 34 
; | fea = 





so that a small correction is required for these modes. It is interesting to 
note that A = nz/a, A = (n+ })/a extrapolate exactly to the modes on a 
stretched string starting with the second mode. As might be expected, no 


analogue exists for the lowest string mode. 
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APPENDIX I 


We have 
" _ ° on 
AzE cos*a [(cos a)é’ + (cos a)r], (7) 
, dx 
\ at ‘+ (sin a)r] (8 
A- -0S* COS 4)7) =(SIN Q)T), ‘ 
n co Xa )n ( ) 
é’cosa+7’sina = 0. (9) 
Expanding (7) and (8) 
AE , rot , : 5 
. €” cosa—€’ sina+7’ cosa—Tsina, (35) 
COs*a 
A*1 : , . 
“I n” cosa—7’ sina+7’ sin a+7 cosa. (36) 
COS" 


Eliminating 7’ from (35) and (36) 


A(sina)E  A*y si ‘s » - oo ver . , - 
(sin a cos a)£” — (cos*a)n” —(sin?a)€’+-(sinacosa)n’—7. (37) 
COS“a COS & 


Eliminating n’ and 7” by use of (9) and rearranging, (37) becomes 


l gr, Manta, Meine 2: (38) 


7 sin x (cos xe” — - aang 
7 Ccos*a cos a 


sin « 
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Substituting (38) into (8) and again using (9) we obtain 


9 a 2) ce 

= 2 n 2cosa ol , sing 

(cos a)&” ——— € +( a OOS Op ———— )e+a2= ——— £ Q, (39) 
sin « sin?a COS*ax/ cos?a 


which is the desired third-order equation in €, 


APPENDIX II 


For the case when |q| Xo < 1 we modify equations (7), (8), and (9) 
, a dé 
WE cos" [eos x — +7 cos |, (7) 
’ dx dx 
| _ dy 
dn cos" COS a rTSM a], (8) 
dx dx 
dé ; dn 
cos a—-+sin a 0 (9) 
da dx 


by replacing sina by « and cosa by 1 and thus obtain: 


Pe dé dr 


2 ~ —- — |+-——-— aT, 
vS ~ da? “dada (40) 
, d*n dy dr 
> Le, ( 
A°y da? = "in sa 
dé " dy _—" (42) 


da * da ~ 


In the limit of a flat equilibrium catenary, one might expect that the chain would 
behave similarly to a stretched string. If such is the case, then it is not hard to see 
that A O(1/a»), and we here make that assumption. The approximate solution then 
obtained proves to be consistent with this. 

In order to make evident the order of magnitude of the quantities involved, we 
introduce a new independent variable 


w= ra (43) 
so that equations (40), (41), and (42) become 
: dé dé dr w 
A2E A? — =. —w — +) —~T, (44) 
> ” dur dw dw AX 
d*r dy dr “3 
d?n 2 } —w l w Te (45) 
dw dw dw 
dé dy : 
/ x w } 0. (46) 
dw dw 
Since we are assuming A 1, let us consider solutions in the form of power series 
in A“! ae 
é > on A af" 
0 
7 = 270A", (47) 
0 
T 3 T,X". 
0 


Then equations (44), (45), and (46) yield sets of equations relating the terms of (47). 
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‘ dE, 
(30) f, = <8 
7 dw 
d?no 
10 = dw? 7 
dé 
“0 _g 
dw J 
é dé, dt, ) 
de dw: ' dw 
d*n 
1 ¢ 
3) dp dw? (49) 
dé dy 
= tt a 0 
) lu dw 
For the limit 1/A ~ ay 1, we are principally interested in the leading terms of 
(47), and from equations (48) and (49) these are seen to be 
4() & 0 
sin sin 
No (w) = (Ax) | 
(41) cos cos >. 
sin 
T 2 (Aa) 270 
(49 cos J 
Equations (50) give the results stated in equations (12) of the paper; equation (13) 
uid verifies the assumption that A O(1/ap). 
It might be mentioned that this procedure can be continued to obtain an improved 
a pproximatior However, in this case, it should be noted that more terms are 


expansion of su 











¥ and cos a in (7), (8), and (9). 














VIRTUAL MOMENTUM AND SLENDER BODY THEORY 


By JOHN W. MILESt+ 
(Department of Engineering, University of California, Los Angeles) 


[Received 18 September 1952] 


SUMMARY 
Munk’s virtual momentum approach to the calculation of the transverse forces 
acting on slender pointed bodies is reconciled with the more rigorous analysis of 
Ward (as extended to unsteady flow) by comparing the dipole fields at infinity, using 
a method of solution due originally to Rayleigh. The virtual momentum of a given 
cross-section is exhibited in terms of the parameters associated with the mapping 
of this section on a circle. 


1. Introduction 

THE problem of the flow of an ideal fluid past a very slender body was first 
attacked by Munk (1), who used the concept of virtual momentum to 
determine the forces acting on airship hulls. More recently this same con- 
cept was applied by Jones (2) to very low aspect ratio pointed wings. These 
analyses were, in large measure, based on physical intuition, but a rigorous 
analysis of the supersonic steady flow problem has been given by Ward (3) 
and extended to unsteady flow by the writer (4). 

[t is not immediately obvious that Ward’s results for a general cross- 
section are the same as those obtained on the basis of the virtual momentum 
concept.{ The following development, which establishes the required 
identity, is based on Rayleigh’s solution of a remarkably similar problem 
in diffraction (5). In the case of a body of revolution, we find an even earlier 
antecedent in Rayleigh’s analysis of diffraction through a circular aper- 
ture (6). 


2. Ward’s formulation 

We consider the supersonic flow of an ideal compressible fluid past a 
slender pointed body of length / located in the immediate neighbourhood 
of the s-axis, the undisturbed flow, U, being directed parallel to this axis. 
The coordinates (s,2,y) are a fixed Cartesian set, with s measured down- 
stream. Then, consistent with the definition of a slender body (see ref. (3) 

+ This paper was written for the Aerodynamic Research group of Douglas Aircraft, Inc. 
(Santa Monica, Calif.), where the author is a consultant. 


+ 


t Originally the writer was of the opinion that the (two-dimensional) virtual momentum 
approach was not generally valid. However, upon being informed by Mr. J. R. Spreiter 
that the approach did give Ward’s result for a winged body, the present analysis suggested 
itself. It is understood that the matter is also being investigated by Mr. Alvin Sacks. 


[Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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for details), the linearized equation for the velocity potential, ¢, reduces to 


d,.+¢,, = 0. (2.1) 


LX vy 
All lengths, times, and velocities are dimensionless and referred to the 
characteristic quantities 7, U, and 1/U, respectively. 
Ward’s solution, as generalized in ref. (4), is developed in terms of the 
complex variable z and the complex potential w, defined by 


x+ty, (2.2) 
w(s, 2, y, t) did (2.3a) 
ao(s,t) Inz+b,(s,t)+ > a,,(s,t)z-™. (2.3 b) 

m=1 


The terms a, In z+-b, are due essentially to the axial flow, and it is found 


3, 4) that 1 DS 


a 


So ogee 2.4 
°— > py’ (2.4) 


where D/Dt indicates substantial time differentiation, viz. 


D()_ &) , A) 
Dt és at 


, (2.5) 
and S (which may depend on ¢ as well as s) is the cross-sectional (transverse 
to s) area of the body. The term b, does not enter the calculation of the 
transverse forces, while the remaining terms in (2.3) can be determined by 
the solution of (2.1) for the cross-flow past the cross-section S at any 
particular station. 

The extension (4) of Ward’s calculation (3) of the transverse force 
components (X, Y) leads to the result 


~ F(s,t) = —(X+iY) (2.6a) 
~ - Pa —_ 
po! il pe) ; nie 


where p, is the mass density in the undisturbed flow, z, the centroid of S, 


and a, the coefficient of z—! in (2.3). 


3. Solution in body coordinates 
We now seek to relate the expression (2.6) to the virtual momentum at 
any local cross-section associated with the passage of a cylinder of this 


cross-section with velocity components (wu, v) through a fluid at rest, corre- 


sponding to the translation 


Dz, 


(3.1) 
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(We emphasize that w, v represent the prescribed velocity of the body, not 
the fluid velocity.) In establishing this relation it is expedient to introduce 
the complex body coordinate 


, 


z = z—2Z,(s, t) (3.2) 
and rewrite the solution (2.3) in the form 
w = ayinz’+b,+w’, (3.3) 
x 
w= > a; - Sea (3.4) 
1 


where w’ is the complex potential due to the two-dimensional cross-flow, 
The coefficient a, is still given by (2.4), as may be shown by repeating the 
derivation thereof in the body coordinates. 

To relate the coefficient a to the desired coefficient a,, we expand In’ 
in inverse powers of z, viz. 


Inz’ = In{z(l1—z,z-1)] = Inz—z,2-1+.... (3.5) 
Then, since the leading term in the expansion of z’—! in inverse powers of z 


is simply z~!, we have ’ 
is simply we have a, = a’, —Ay2,. (3.6) 
Substituting (3.1), (3.6), and (2.4) in (2.6b) yields the force in the form 


oF D , ? ~ 
— Po ll "pe a "- (3.7) 


It remains to determine the coefficient aj, which is the strength of the 
dipole field associated with the cross-flow w’. 


4. Virtual momentum from dipole field 
Let the (dimensionless) coefficients of virtual mass (A, H, B) of the cross- 
section S be defined by the expression 
T = hp, U2l?(Au?+ 2Huv-+ Br?) (4.1) 
for the kinetic energy of the flow field associated with the motion of transla- 
tion prescribed by (3.1). The corresponding virtual momentum (dimension- 
less with py Ul? as reference quantity) is given by 
M = M,+iM, = (Au-+Hv)+i(Hu+ Bo). (4.2) 
In terms of these coefficients, Lamb (7), following Rayleigh’s analysis of the 
analogous diffraction problem (5), has shown that the distant field due to 


the motion is given by 


w’ ~ {[(A+S8)u+ Av]+i[ Hu+ (B+ S)v]}}(272’)+1. (4.3) 


j 


Equating the coefficient of z’ to the dipole strength and substituting V 
from (4.2), we obtain 


27a} = M-4 Sq. (4.4) 
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Finally we substitute a, in (3.7) to obtain the desired result 


oF .» DM 


l 
ds PO” De 


; (4.5) 
which exhibits the transverse force on any local section as the rate of change 


of its virtual momentum. 


5. Mapping solution 
We conclude by constructing an expression for the virtual momentum 
in terms of the mapping of S on a circle. Let 
C= f(z’) (5.1) 
be the conformal transformation that maps S on |¢| = c(s,t). In order to 
preserve the flow at infinity, we impose the restriction 


im 2 f(s’) = 1. (5.2) 


Then, using the well-known solution for the flow (of complex velocity q at 
infinity) past a circle of radius c, viz. 

w’ g6+qet- (5.3) 
where the bar indicates the complex conjugate, we find the total cross-flow 
in the z’-plane by substituting ¢ from (5.1) and subtracting out the flow at 
infinity, namely gz’, whence 

, -[ fy; _ 9/ om - 
w’ = qf (2')—2"]+4c?/f(z’). (5.4) 
The coefficient of z’~! in the expansion of f(z’) will be designated as Res(f ), 
while, by virtue of (5.2), the coefficient of z’-! in the expansion of f- is 
simply unity. Hence we obtain 
1 = (Res f+ qe’. (5.5) 


a 
Substituting this in (4.4) and solving for M, we obtain 


M = (2nc?—S)q+27q Resf. (5.6) 
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THE STEADY FLOW OF A VISCOUS FLUID PAST AN 
ELLIPTIC CYLINDER AND A FLAT PLATE AT SMALL 
REYNOLDS NUMBERS 


By 8. TOMOTIKA and T. AOI (University of Kyoto, Japan) 


[Received 24 April 1952 


SUMMARY 

In Part I the exact solution of Oseen’s equations of motion for the steady flow 
of an incompressible viscous fluid past an elliptic cylinder is obtained, in terms of 
elliptic coordinates. The solution gives as a limiting case the flow past a flat plate. 

The drag experienced by the elliptic cylinder is discussed in Part II, with special 
reference to both the pressure drag and the frictional drag. It is thus found that, 
as long as the calculations are based upon Oseen’s equations, the total drag can be 
analysed into pressure and frictional drags proportional to the axes of the cylinder 
respectively for any value of the Reynolds number. Further, expansion formulae 
correct to the fourth power of the Reynolds number are derived for the drag on thi 
elliptic cylinder and on a flat plate, and numerical values of the drag coefficient ar 
computed from these formulae. 

Part III deals with the general expression for the stream function representing 
the flow past the elliptic cylinder. An approximate formula for the stream function 
for the flow in the neighbourhood of the flat plate is then obtained. The velocity 
profiles as well as the stream-lines for a flat plate edgewise to the stream are drawn, 
while for a flat plate at right angles to the stream, the stream-lines show that there 
exists a pair of standing eddies behind the plate. 

Lastly, the pressure distribution on the surface of the elliptic cylinder is calculated 
in Part IV. 


1. Introduction 

As is well known, the steady flow of an incompressible viscous fluid past a 
solid obstacle can be dealt with successfully on the basis of Oseen’s linearized 
equations of motion when the Reynolds number is small. In previous 
papers (1, 2) we have investigated in detail the steady flow past a sphere 
and a circular cylinder by making use of the exact analytical solution of 
Oseen’s equations, and it has thus been found that the actual flow patterns 
observed around these obstacles can be deduced theoretically from this 
solution. 

Similar investigations on the steady flow past an elliptic cylinder or a flat 
plate based on Oseen’s equations have been made by several writers, but 
they were mainly concerned with discussions on the drag experienced by 
the obstacle. For instance, an approximate expression for the drag on an 
elliptic cylinder or a flat plate was obtained by Harrison (3), by Bairstow, 
Cave, and Lang (4), and by Piercy and Winny (5). Davies (6) deduced an 


[Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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approximate formula for the drag on a flat plate edgewise to the uniform 
stream, but, as will be shown in this paper, his results are erroneous. Further, 
Sidrak (7) has also recently derived an approximate formula for the drag 
in an elliptic cylinder and on a flat plate, both placed parallel to the stream. 
3ut his analysis is vitiated by several errors and his results are unreliable. 

As far as we are aware, no other detailed theoretical discussion based on 
Oseen’s equations seems to have been made in the case of the steady flow 
,round an elliptic cylinder or a flat plate placed parallel or perpendicularly 
to the uniform stream. The main object of the present paper is to discuss 
this problem and in particular to obtain the drag experienced by, and the 
flow patterns around, the obstacle by making use of the exact solution of 
Oseen’s equations. The pressure distribution on the surface of the obstacle 
is also computed. 

Detailed analytical calculations are given only for the case of an elliptic 
cylinder placed parallel to the stream, since the analysis for the case of the 
cylinder placed perpendicularly to the stream is similar to this case. How- 
ever, essential parts of the analysis for the latter case are given in the 
Appendix. 

Part I. SoLtuTIon OF OSEEN’S EQUATIONS 
2. The general solution 

We consider a two-dimensional cylindrical obstacle fixed in a steadily run- 
ning uniform stream of an incompressible viscous fluid of infinite extent. Let 
O(x, y) be the rectangular coordinates having the origin O within the obstacle 
ind the axis of xin the direction of the velocity U of the uniform stream. We 
denote by U’ +-u, v the components of the fluid velocity at any point, so that 
1,vare the velocity of perturbation which become vanishingly small every- 
where at a great distance from the obstacle. If squares of u, v are omitted, 
we get, from the Navier-Stokes equations of motion, the well-known 
linearized equations of Oseen 

U ~ (u,v) “{- = )p+oV%(u,0), (1) 

Cx p\ox cy; 
where p is the pressure, p the density, v the kinematic viscosity, and V? 
stands for ¢?/éx?-+-e?/ey*. The fluid being assumed to be incompressible, 


the equation of continuity is then 


ou , w = (2) 
ox oy 
If we put k = U/2v, these equations are satisfied by 
l / C l ¢ 
a. Sah. ee ee (3) 


’ ‘ ~ i 
Ox 2k da dy 2k ey 
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and p= pu, (4 
provided that V*o = 0, (5) 
and Vix—2k—* = 0. (6 


Consider an elliptic cylinder whose major and minor axes are of length 
2a and 2b respectively, and let this cylinder be placed parallel to the 
uniform stream with its centre at the origin of coordinates and with the 
major and minor axes along the axes of x and y respectively. Let é and » 
be elliptic coordinates defined by 


a- iy = c cosh(é- in), ; 


there c = 2__p2 
where c = ,/(a?— 6). 

Since the flow past the cylinder is symmetrical about the axis of x and 
perturbations due to the presence of the cylinder must vanish at infinity, 


the appropriate solution of equation (5) in terms of € and 7 is given by 


x 
‘ . l 
d= cUa,é—cl > x, e-"— COs nn, (8) 
n 
1 


n 


where the «,,’s are constants of integration. 
To solve equation (6) put x = e*"S(£)H(n); the equation becomes 
separable in the variables and we are led to the modified Mathieu equa- 


tion for &(€) and the Mathieu equation for H(7), namely 


a - (A+ > cosh 2) == 0, (9) 
= 4+ {A+ ~*~ cos 2n|/H = 0, (19) 
dy? : 


where « stands for kc. For a discrete set of characteristic values of \ the solu- 
tions of equation (10), which are periodic and even in », are expressed as 
0 
ce,,,()) = > (—1)™*"AL™ cos 2ryn 
r=0 
(11) 
x 
CComii(n) = > (—1)™* BEn{P cos(2r-++- 1) 
r=0 
where the coefficients A” and B’™ are functions of «. The solutions of the 
second kind of equation (9), which are associated with ce,,(y), can be 
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expressed in terms of the modified Bessel functions J, and K, in the 


n 


forms:t 


FEK,,,(£) { (2m) S Ay" K,(4xe§) I (Axe ) 
— £4 
|. (12) 
wRK $ »(2m+1).- 
FEKansi(8) = remrD > Ba" 





(K...(4uce€)I,.(4uce-€) — K,(4uce£)I,, 4(frce-)} | 


These series involving Bessel function products are known to be rapidly 
convergent for small values of x. 


Thus the appropriate general expression for x is given by 


TK cosh & cos _ nh Pe _ € 
X Ue 7 , z BI EK,,,(€)ce,,(7), (13) 
m=0 
where the 8,,’s are constants of integration. 

The components u, v of the perturbation velocity and the pressure p at 
ny point in the flow can be obtained by inserting the above general 


expression for ¢ and x into (3) and (4) respectively. 


3. Determination of the constants of integration a, and £,, 

The constants of integration «, and £,, can be determined from the 
boundary conditions at the surface of the elliptic cylinder. To do this it is 
convenient to introduce the components of the perturbation velocity v¢ 
and v, normal to an ellipse € = constant and a hyperbola 7 = constant 
respectively, which are orthogonal. These components are given, in terms 


of d and X; by 


| Cdh | cy ] ° 
v # xX. x sinh € cos 
' ch c& © 2xh C& Oh 
(14) 
| Cob | C q ] — 
v. X x cosh € sin 
ch O71 2h ¢ n h 
] 
9 l < v, y) 9 9 
where h? ee cosh?é— cos? 7. 
C* « (¢ 7) 
The tations used here for the Mathieu function and the modified Mathieu function 


similar to those in N. W. McLachlan’s Theory and Application of Mathieu Functions 
Oxford, 1947), except that our FEK,,(é) differs from McLachlan’s function Fek,,(€) by a 


Onstant multiplier. 
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Substituting the expressions (8) and (13) for ¢ and x into (14) we have 


1 a A 1 y x 
v= — — x, e-"€ cos nn + eX cosh £ cos 7 > x 
é h > n 7-1 2h Bin 
m=0 


n=0 





x {FEK(, (é)ce,,(7)—« sinh € FEK,,, (€)ce,,,(7)cos 7} 


U x UV x 
Sy pela “ty . nt o* | « cosh € cos 7 es Wa 
ye as > ne" sin ny + —— ¢ 5one Bn & 
m=0 


n=1 


‘e) (15) 





< {FEK,,,(€)ce,,(4)-+« cosh € FEK,, (€)ce,,,(7)sin 7} 


where dashes denote differentiation with respect to € or 7. 

We now rewrite these expressions in forms which are suitable for fitting 
the boundary conditions at the surface of the cylinder and for calculating 
the stream function. For this purpose we expand the carters sides in 
(15) into Fourier series. Let us firstly introduce functions F,, ,,(€), @,, (8), 
such that 


] 


aK 


gramgeny K{FEK,,,(€)ce Cn(n)— K sinh & F K K,,(€)ce,,(1)eos n} 


™ de Fn, n(€)COs nn 
: (16) 


~ ex cosh Fos 1 KEK (E)ce;,,(7)-+« cosh  FEK,, (€)ce,,(7)sin 7} 
mK 


= > G,,,,(é)sin ny 


ro 
n=1 
Then F,, 9 is clearly given by 
F | PEK, (£)Ce’, (¢)FEK’, (é)Ce., (¢)} 
m,0 = a Dic UK, (€) CAS) Cinlé ) CnlS)5, 
mF 
here Ce (¢ 1 [ x cosh £ COs 7) ae ] 
where e,,(§) — ss , ce,,(7) dy, 
d 


the modified Mathieu equation (9) associated with ce,,(7). Therefore F,,o 
is a constant independent of € but dependent upon m. This constant value 
can easily be evaluated by assuming a large value for € and replacing 


FEK,,(€) and Ce,,(€) by their respective asymptotic expressions, thus: 


FEK,, (€) ~ Jz) z Ce,,,(€) ~ ce,,(0) l e=. 
22 4 (2772) 


where z = ixef. Hence we have 


Pe aes ce,,(0). (17) 


m,0 ¢ 
2K 


and represents the definite integral form of the solution of the first kind of 
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Straightforward calculations give the expressions for F,, ,(€) and G,,, ,() 
(n = 1, 2, 3,...) as 


(£) LF K Kn (E)Dinn (¢) —-FE K,,.(€) Ln, n(€)$ 


K 


2 Pes , 
G(E) FEK,,(€){« cosh € J, »-1(€)— Ny, n(€)—K Cosh EL nay 


m,m 
- 
where 


1 (€) } S ( l ym tr Ave) a _op(K cosh €)+ I, orl cosh €)} ; 
J (19) 


2m,n\> 


n—2r-1(k cosh €)+- J), .5,.,(« cosh €)} 


Im s1n(8) = 4 > (—1ym Bees PU 


With these expressions for F,, 9, F,,,(€), and G,, ,(€) (w = 1, 2, 3,...) the 


m,0° m,n 


expressions for vg and v, become 


‘_ 2 
v, — S on ng > Bin Fon n(é)}eos nn 
n=0 * m=0 (20) 
U = 
" r h ‘3 jan en" -> Brn Gaal n as sin nn 


The constants of integration a, and 3 contained in (20) are determined 
by the boundary conditions. Since, however, the conditions at infinity are 
satisfied automatically, we have only to consider the conditions at the 
surface of the elliptic cylinder. 

Let the equation of the surface of the elliptic cylinder under consideration 
be € = &), then the boundary conditions become 


VU: sinh € cos 7, 
i 


; 
+ cosh € sin », 
? 7 : 


at € = €. Using the expressions (20) for vg and v,, these conditions yield 


immediately 


t 


I | 
Xo Se (0)£,, | 


2k on m | 
m=0 | 
2 sinh € n l | ‘ 

Xp > Bin Fnn(€o) °0 \ ) i (21) 
m=0 . 10 (n ye 
| 
od > BG (é) { cosh &, (n 1) 
moh ZmnlSo) =o (mw = 2, 3,...) | 


Eliminating the «,,’s from the second and third of these equations, we have 


— £o 
S BSF (é) (é,)) (—eee (n 1), 
hom nk mn\>0 n,n\> € 
0 ; ‘ms 0 (n = 2, Oe 
and the 8,,’s can be determined by solving this system of simultaneous 


# : 
linear algebraic equations. 
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Part II. THe DRAG ON THE ELLIPTIC CYLINDER AND ON THE 
FLAT PLATE 


4. The general expression for the drag 
We shall next calculate the drag experienced by an elliptic cylinder and 
also by a flat plate. As was mentioned in (1), two different methods are 
commonly used; either the drag is obtained by summing up the viscous 
stresses exerted by the fluid unon the surface of the obstacle, or it is caleu- 
lated by applying the theorem of momentum to an infinite mass of fluid 
surrounding the obstacle. 
If the analysis is based upon Oseen’s equations, the first method gives 
the drag D on a cylindrical obstacle in the form 
D= —U | ts, (23) 
én 
where integration is taken round the circumferential curve s of the cylinder 
and é/én denotes differentiation along the outward normal n to s. The 
second method gives r ad 
D = —pU | “? do, (24) 
on 
where integration is taken round a large closed contour o everywhere at a 
great distance from the cylinder and @/én’ denotes differentiation along the 
outward normal n’ to o. 


For the elliptic cylinder, € = &,, (23) becomes 


whilst (24) becomes D = —pU | (2) dn, 
© aes 


o being taken to be a large ellipse confocal with the cross-section of the 
cylinder. 

Making use of the expression (8) for d the above two formulae both give 
for the drag on the cylinder the value 


D=-— 2mpU*cag, (25) 
which, by the aid of the first of equations (21), can also be written in the 
form : 

D = 2mpU Y ce,,(0)8,,; (26) 
m=0 


where denotes the coefficient of viscosity of the fluid. This gives the 
general expression for the drag on the elliptic cylinder placed parallel to 
the uniform stream. 
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Defining the drag coefficient Cp, by Cp = D/(pU?.2a)-1, we have 


= ae 
Y — vs 97 
“p R ce,(0) Bn (27) 
m=0 
where R 2aU |v is the Reynolds number. 


The general expression for the drag on an elliptic cylinder placed perpen- 
dicularly to the uniform stream can be obtained likewise, as will be shown 
in the Appendix. 


5. The pressure drag and the frictional drag 

As was shown in (1), the total drag D on a solid can be analysed into the 
pressure drag J, and the frictional drag D,. In the case of a cylindrical body, 
adopting Oseen’s approximation, the general formulae for D, and D, are 
given respectively by 

D, pU 1? ds, D, = —pU | m a ds, (28) 
oy 

where (J, m) are the direction-cosines of the outward normal to the cross- 
sectional curve s of the cylinder. 


In the particular case of an elliptic cylinder these become 


27 
es * [¢ 
D, pl csinh €, | a cos 7 dyn, 
J \ex}e-¢, 
0 i 


on 


D, pl ‘c cosh €, [ (=) sin n dn, 
J oy é=aé 
0 i 


=Co 


and substituting the expression (8) for ¢ and carrying out integrations, we 


1 
nave 


D, 2rpU cay € fosinh &,, 
D, 2rpU*c Xp ¢ £0 cosh &,. 
\dding these we immediately get again the expression (25) for the total 
drag D. 

Further, remembering that a = ccosh€, and b = csinh &€, we arrive at 
an interesting result that 


b a 
D, D, D, D. (29) 

, a+b a+b 
This result shows that, if caleulations are based upon Oseen’s equations of 
motion, the total drag on the elliptic cylinder placed parallel to the uniform 
stream is divided into the pressure and frictional drags in the ratio b:a, 
whatever the value of the Reynolds number may be. 
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On the other hand, in a case when the elliptic cylinder is placed perpen. 
dicularly to the uniform stream we have 


; a D. b 


— - at 
a+b (30) 


) ), = -D. 
2 IS a+b 


Thus in this case the contribution of the pressure and frictional drags to the 
total drag are in the ratio a:b. 
We can now deduce the results for the special cases of a flat plate and 
a circular cylinder. Thus, for a flat plate edgewise to the uniform stream 
we have, by putting b = 0 in (29), 
D,=0, D=D, (31) 


while for a flat plate at right angles to the uniform stream we get, by putting 
= 0 in (30), 


D=D, D=0. (32) 

Further, in the case of a circular cylinder (a = 6) we have, from either 
9 ¢ > 

of (29) and (30), D, jie LD, D, — 1D. (33) 


Thus, as we already pointed out in (1), the total drag on the circular cylinder 
is equally divided into the pressure and frictional drags. 


6. Expansion formulae for the drag 

In order to calculate numerical values of the drag coefficient directly from 
the general expression (27) we must first determine the 8 
solving the system of simultaneous algebraic equations (22) for each given 


m 8 numerically by 
value of the Reynolds number; the work is thus rather cumbersome. It is 
therefore desirable to derive an expansion formula in powers of the Reynolds 
number which may be conveniently used for the purpose of computing 
the drag. 

To this end, we have to express 8,,, as well as ce,,,(0) in powers of «, or in 
powers of the Reynolds number R. To do this, we must expand the quanti- 
ties F,, »(€o)— Gn.n(€o) in (22) in powers of « by using the series expansions 
for the coefficients A’” and B’” in ce,,(7) as well as the series expansions 
for the modified Bessel functions J, and K, which occur in F,, ,,(&) and 
Gy»(€). Assuming R to be small it is thus found that By, B,, By,... are 
of order 1, R?, R4,..., respectively,+ while it is easily seen that the quantity 
ce,,(0) for any number of m is of order unity. 

Using the above estimated orders we have derived an expansion formula 
for the drag coefficient correct to the order of R*. The expansions for 8, 
or Fi, n(€o)—Gn.n(€o) necessary for our derivation need not be given here 
and only the final result will be stated. 


+ In order to obtain the expansions for Bo, B;, and B, correct to the order of R%, it is 
necessary to solve the first five equations (n 1-5) in (22). 
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For convenience we introduce a parameter defined as 
1—(b/a) 
C= rr eert, 
1+ (b/a) 
which depends on the shape of the elliptic cylinder. Then the expansion 
for the drag on the elliptic cylinder whose major axis is parallel to the stream 


is given by 


dor R? 
Cy 
RS 32(1+oc)?8 
{S2?—4(1+ 20—o0?)S+ 4,(15—4o0— 180?— 120%—o*)}— 
R! ; — P 5. or 
—{S4—4(4— 30" 203) S3— 3.(31la0+ 100?— 603+ 1004+ 0°) S?+ 


327(1+0)4S 


1. (280-+ 11850— 98702 — 206507 — 9004+ 52505+ 850%+ 307) S— 


1 


sapq(225— 1200—52407— 21603-39004 + 4400°+- 18008 + 2407 +-o*)} 
(34) 
R 
where S 1(1+oc)—y—log , 
a Y os 8(1 T a) 


and y = 0-57721... is Euler’s constant. 

Asimilar expansion formula for the drag coefficient has also been obtained 
for the case when the elliptic cylinder is situated in such a manner that its 
major axis is perpendicular to the direction of the uniform stream. The 


result will be given here for comparison. Thus, 


Aor | k? 
Cp \1 7 
RS| 32(1+0)28 
{S?—4(1—20—o0?)S+ (15-4 4a— 1807+ 1203—o*)!— 
R4 
} (S41 (4302+ 203) S34 4(31o0— 100?—603— 1004+ 05) S?2+ 
322(1+-c)4S?' os 7 
1 (280— 11850— 98707-+ 206503 - 9004 — 5250°-4 85a°— 307)S 
»3p4(225-+ 1200 —5240?+- 21603 + 39004 — 4400° + 1800°— 2407 + 08)} 
(35) 
where S }(1—o)—y—log- 


8(1-+-o) 

It may be noted that the first approximation Cp) = 47/(RS), obtained 
by taking only the first term in the square bracket in (34) or (35), is identical 
with that given by Harrison (3) and also by Bairstow, Cave, and Lang (4). 
Apart from this, our formulae (34) and (35) for the drag on the elliptic 
cylinder of any thickness ratio appear to be new. 
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In a recent paper (7) Sidrak has obtained similar expansions, as far as 
terms of order of R?, for the drag on an elliptic cylinder and for the skin 
friction of a flat plate, each placed parallel to the direction of the uniform 
flow. A careful examination, however, has revealed several errors? in his 
analysis, so that his work is unreliable. 

Expansion formulae for the drag on a flat plate and a circular cylinder 
can be readily obtained as limiting forms of the above formulae (34) and 
(35). Thus, in the first place, when a flat plate of length L (= 2a) is placed 
edgewise along the uniform stream, we get, by putting o = 1 in (34), an 
expansion formula for the skin-friction coefficient of the plate correct as far 
as terms of order FR‘, in the form 


4a l aa. 4 5, ol 
_ a: = §44.1,93_ 23 92__ 133 9__ 25 
( Dp = RS P sea gl S i3) 128 2 (S*+ i3S 34 3605 14) 1282|° 
(36) 
, ee 
where S = 1—y—log +t 3:1954—log R, 
) 


and the Reynolds number, R, is defined as R = UL/v. 

This formula may be compared with the similar expansion formulae ob- 
tained by other writers. In 1933 Piercy and Winny (5), following Bairstow, 
Cave, and Eang’s method of analysis (4), obtained a second approximation 
(correct as far as terms of order R?) for the skin-friction coefficient of the 
flat plate. When expanded fully in powers of R their formula becomes, in 


our notation . 
: An Bice  «. 

Cp = =a |1——(S*—S—f) |. 

RSS| S 2128, 


Thus it is seen that up to the order of R? our formula (36) is identical with 
Piercy and Winny’s formula. One of us (S.T.) has extended Piercy and 
Winny’s analysis to a third approximation, and it has been confirmed that 
our formula (36) is correct up to the term in R‘, 

Davies (6) has also derived an expansion formula up to the order of R? 
for the skin friction of the plate placed edgewise to the uniform stream. 
But, as is noted by himself, no agreement is found between the expressions 
by himself and by Piercy and Winny. Recently Yosinobu and one of us 
(S. T.) have examined carefully Davies’s analysis and found that his result 
is incorrect. It seems certain that this is due to bad convergence of the 
series-form of the modified Mathieu function used by him. Indeed, it has 
been found that if use is made of the function FEK,,,(&), defined by (12), 


m 


for the modified Mathieu function, we can derive, by pursuing Davies s 


+ One of the errors occurs at the beginning of his analysis, the expressions for the rectangu- 
lar components of velocity of perturbation u and v being, in fact, erroneous. 
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analysis, the correct second approximation for the skin friction of the flat 
plate, which is identical with Piercy and Winny’s formula.t 

Next, if we put o | in (35), we get an expansion formula, correct as 
far as terms of order R‘, for the drag on the flat plate placed at right angles 


across the uniform stream. We find 


; Arr l (ge yi] R° 4.010931 192_19 Re 
Cp = pal — gi tT 4) Tog ge a+ 85 — Sis) Tog3 ; 
(37) 
with S Y log i R — 2-1954—log R, 


where, as before, the Reynolds number, R, is defined as R = UL/v. This 
formula does not seem to have been obtained before. 

Lastly, an expansion formula for the drag on the circular cylinder of 
diameter d = 2a can be obtained by putting o = 0 in either of our preceding 
formulae (34) or (35). The result is 


bcr - ic eee Bom oo. 3s. eee 
Up | L—(S?—-38 +98) 55 — aq (S43 + eS — ase) aaa] (38) 
tS S ‘ 32 S* ; i ~ ——e 
where S 4—y—log4R = 2-0022—log R, 


and R = Ud/visthe Reynolds number. As should be expected, this formula 
is in complete agreement with that obtained by us previously (8). 


7. Numerical values of the drag coefficient 

3y making use of our expansion formulae (34)—(38) we shall now calculate 
numerical values of the drag coefficient C,, for the elliptic cylinder (including 
the flat plate and the circular cylinder as limiting cases) at small Reynolds 
numbers. The range of values of the Reynolds number in which these 
formulae can legitimately be used was first estimated as for the case 
of the expansion formula for the drag on the circular cylinder of our previous 
paper (8), and it has been found that all our formulae can be used with 
sufficient accuracy for values of the Reynolds numbers less than about 4. 

Numerical values of C;, thus calculated are given in Table I and are shown 
graphically in Fig. 1. It is of some interest to note, as will readily be seen 
from Table I or Fig. 1, that the value of the drag coefficient of the circular 
cylinder at any small Reynolds number is always larger than the corre- 
sponding value for the elliptic cylinder or the flat plate, each having the 
same length as the diameter of the circular cylinder. This ‘rather curious’ 
result has been already pointed out by Bairstow, Cave, and Lang (4, p. 429) 
by using the first approximation C,) = 47/(RS). 


+ Our analysis will be given elsewhere. 
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1. Curves showing the variation of the total drag coefficient Cp with 








Reynolds number R. Curves I, II, III, IV, and V correspond respec- 
tively to the cases D, E, G, C, and A in Table I. 
TABLE I. Values of Cp 
Case A. Flat plate 
Case B. Elliptic cylinder (b t) } parallel to the stream. 
Case C. Elliptic cylinder (6 4) 
Case D. Circular cylinder. 
Case E. Elliptic cylinder (6 oe Sener to the stream. 
Case F. Elliptic cylinder (b 1) 
Case G. Flat plate 
Case Case B Case C Case D Case E Case I Case G 
4 7°61 5°47 9°24 10°63 10°25 10°10 10°04 
oe 5°61 6°31 6°95 8°13 7°82 7°70 7°66 
8 4°54 5706 5°73 6-78 6°51 6-41 6-38 
3°87 4°43 4°95 5°93 5-69 5:60 5°57 
2°39 2°05 3°24 4°07 3°88 3°84 3°83 
1°84 2°25 2°62 3°47 3°22 3720 06] |= 318 
155 2°04 2°35 2-92T — 2°89 | 2:84 





+ This value 2-92 has been calculated by the general expression (given by (48) in 
our previous paper (1)) for the drag on the circular cylinder. 


In the case of the flat plate placed edgewise to the uniform stream, the 
value of the drag coefficient for R = 4 has been computed by several writers 
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by using different methods. Thus Bairstow, Cave, and Lang (4) have 
obtained graphically the result C,; = 1-527, while Southwell and Squire (9), 
using the extended Oseen’s equations due to them, obtained numerically 
the value of C, = 1-563. These values may be compared with the value 
( = 1-551 calculated from our formula (36). Although there appear to be 
slight discrepancies among these values, the agreement between them is 
nevertheless rather good, especially if we take account of probable errors 
which can arise in graphical computations. Lastly, it may be added here 
that Pierey and Winny (5) obtained the value Cp = 1-499 for the skin- 
friction coefficient of the flat plate by using an asymptotic solution of 


Oseen’s equations. 


Part III. THe FLow PATTERNS AROUND THE ELLIPTiIc CYLINDER 
AND THE FLAT PLATE 


8. The general expression for the stream function 


We shall now calculate the general expression of the stream function for 
the flow past the elliptic cylinder by using the exact solution of Oseen’s 
equations. If we introduce the stream function y’ for the velocity of 
perturbation so that wu és’ /ey and v oy’ /ex, then vz is expressed in 
terms of x’ as », (ch)-* eb’ /@n. 

If we denote the stream function for the flow past the fixed cylinder 


in the uniform stream by %, we have 


7 
ub Uy+ | chvz dn, 


0 


provided that the constant of integration is so chosen that the axis of x 
forms a part of the stream-line & 0. 
Inserting in this equation the expression (20) for vz, and taking account 


of the first and second of equations (21), we get 


% = cU (sinh €—eé-£ sinh €,)sin n-4 


UY  Y BulEnnlE)—e" OF y 9 (Eo)} ——-- (39) 
n=1 m=0 n 

This is the required expression for the stream function for the flow past the 

elliptic cylinder whose major axis is parallel to the uniform stream. As will 

readily be seen, the first term on the right-hand side is the stream function 

in the case of a perfect fluid and consequently the second term represents 

the effect of viscosity of the fluid. 
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If we introduce the non-dimensional stream function %, defined as 


%, = %/(Ua), we have 


#, = sech &,(sinh ¢—ee-€ sinh €,)sin 7+ 


» << sin m1 
= f (5-8) Ff ’ ] 
+sech £5 > > Bini Finn(€)—e" bo-€ Fin(€o)} - (40) 
n=1 m=0 n 
A similar general expression for the stream function in the case of flow 
past an elliptic cylinder placed at right angles to the uniform stream will 
be given in the Appendix. 
g 


9. Approximate formulae for the stream function in the case of 
flow around the flat plate 


We shall now deduce from (40) an approximate formula for the stream 
function for small Reynolds numbers R, which will be used to compute 
the flow patterns. 

To this end, the general expression for the stream function will be 
expanded in powers of R as far as terms of order R?. We first derive 
expansions in powers of R of the functions F,, ,(€) given by (40), and of 
the quantities 8,,. Expansions of the £,,,’s have already been obtained. 
Expansions for the functions F,, ,,(€) have been obtained by using the series 
definitions of the modified Bessel functions J, and K,, (which occur in 
the expressions for F,, ,(€)) on the assumption that their arguments are 
small. Accordingly such expansions are justified when the variable € is 
also small; our approximate formula for the stream function may be 
adequately used for the flow near the obstacle. 

In this way, we have obtained an approximate formula for the stream 
function for the flow past the elliptic cylinder of any thickness ratio and 
we have deduced the corresponding formula for the case of a flat plate by 
making the parameter €, tend to zero. Since it is of special interest to 
discuss the flow patterns around the flat plate, only the results for this case 
will be given below. 

For the flow past a flat plate placed edgewise to the uniform stream an 
approximate formula for the stream function, correct to the order R?, is 


f= {((4y,3+5,,3 €)e%* + (Ay) +5,, E)eE+ (ay, +b, _, é)e#+ 
+ (a, 3+, -3&)e 3é}sin q+ ((42,2 +b. 2€)e* + de9+ (ae, a+b, »€)e-*F}sin 29+ 
(433+. Ee% + (ag, +b Eye + (43 _,+3_, Eje$+ 


+ (dg _3+b, _s )e-**}sin 3n, (41) 





where t 
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where the coefficients a,,,, and b,, , are given by 





i ee ree | 
“8 = ~~ 1024” aE eae | 
} 
Re(, 13 6 ee ie oe 
“us si isa) “aa = —~"L-1 = 98° 1024\8* 383)’ | 
?? 13 
1-1 10 Al a | 
e213 
al, 3 51 > 128 
2 | | 
— As (1 , Ring -b._. = , = | 
- . 64 8 - 3 648 | 
R I l | 
ly 9 at 3} | 
| 
R? 2 R? | 
— " b. — —f = — » 
133 cant - "38 3,-3 = 3072 8” 
R? |] R2 1 
fle a b. == a \ — ad 
‘31 = 3048 S ™ 8-1 = Toad 8? 
ai ee (1 3) 
, 1024 38 
?? 17 
ala » I (1 ). 
ie 1536 12S 


(42) 

the quantity S being defined in (36). 
On the other hand, for the case of a flat plate placed at right angles to 
the uniform stream, the corresponding expansion formula, correct to the 
order of R2. is of the same form as (41), but in this case the coefficients 


a,, and b,,.. are given by 


R? l R2 |] 
"ee Pak : b,. ij P ‘ 
; vagal 3) 1,3 “1-3 = 1094 8 
l i 5 | ] R? [1 1\ | 
fl, i (1 a _- bi, b, 1 _— ( +t a 
IS ° 256 SS 4,87} . , 28 1024\S — S82)" | 
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R(,, 1 
72,0 = 35 +59) 


— | : R2 l 
433 = — sora? - as} bss = 5.3 = 3079 8” 
my R? 1 
a3 = ark +5) sa = Psa = — Tog 
3.2? 
a34 = 1024 


, R2 1 7 
ade == — re 
3-3 768 248 


where the quantity S is given in (37). 


10. The flow patterns around the flat plate 

Using the above approximate formulae for the stream function, we shall 
now discuss the flow patterns around the flat plate. In the first place, we 
shall compute the velocity distributions along the lines perpendicular to the 
plate placed edgewise along the uniform stream, in a case when the 
Xeynolds number is equal to 4. The results are given in Table II, where 
the values of the velocity components v, = U-+u and v, = v along the 
lines x = —1-5, —1-0, —0-5, 0, 0-5, 1-0, 1-5 are given, taking, for con- 
venience, the uniform velocity U to be unity and the length L of the plate 
to be 2. The velocity profiles are shown in Fig. 2, in which the stream-lines 
around the plate are also shown. 

It will be seen immediately from this figure that when the Reynolds 
number is small, the flow near the flat plate resembles closely the shear 
flow with uniform velocity-gradient, except in the vicinity of its edges. 


TABLE II. The velocity distributions around the flat plate 


edgewise along the stream in the case when R = 4 








x I°5 x I°o x ors x ° x o's x I*o I 

y zr v / r / vy I “Vv r v 
fe) 054 oO |0 ° ° ° ° ° fe) ° ° ° 0°53 fe) 
orl oo 0°90 |0°29 0°09 0°13 O'OI |O°"IO O00 | O'I! 0°00 | 0°22 o°0o7 | O°5 ° 
o*2 06 o'r; O°39 O'13)/0°24 O02 O'IQ 0'00/0'2I O°o!l | 0°33 oll | o°6 fe) 
o°3 O47 O15 0°34 0°04/0°29 O°’O! | 0°30 0°02 | O°42 or1l3 
o"4 0°53 O10 0°43 0°00/0°37 Oro! 0°39 0°04 | 0°50 o°160 

5 0°60 o'18/) 0°50 0°08/ 0°45 O°0l | 0°47 0°05 | 0°55 O'l7 





The stream-lines around a flat plate placed at right angles to the uniform 
stream have been drawn in case when R= 4. The flow patterns thus 
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Fig. 2. Caleulated stream-lines and velocity profiles for R = 4 around a flat plate 


placed edgewise along the uniform stream of velocity U. In the half-field of flow 


above the central stream-line %, = 0, the stream-lines correspond respectively to 
%, = 0-01, 0-05, 0-1, and 0-2, 











Fic. 3. Calculated stream-lines past a flat plate at right angles 
across the uniform stream for R = 4. In the half-field of flow 


above the central stream-line y, 0, the stream-lines corre- 
spond respectively to x, 0-01, 0-05, 0-1, 0-2, and 0-3, while the 
stream-lines in the region of eddies are at y, 0-0005 and 


0-001 respectively. 


obtained are shown in Fig. 3: it will be seen that there exists a pair of 
standing eddies behind the plate. It is worth noticing that when the flat 
plate is placed at right angles to the stream, a pair of standing eddies will 
always appear behind the body even when the Reynolds number is very 
small. 
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Part IV. THE PRESSURE DISTRIBUTIONS ON THE SURFACE OF 
THE ELLiptic CYLINDER 


11. The pressure distributions 


Lastly we shall compute the pressure distributions on the surface of the 
elliptic cylinder. For this purpose we first derive from the exact solution 
of Oseen’s equations the general expression for the pressure on the surface 
of the elliptic cylinder. Now, the pressure p at any point in the fluid is 
given by (4), together with the expression for ¢, given by (8). Since, how- 
ever, 6¢/€x > 0 as € > 00, p represents the pressure relative to that at 
infinity. Thus, if we denote the absolute pressure at any point in the fluid 
by p also and the corresponding pressure at infinity by p,,, we can write 
, Od 

ex 


P—Px = pl 


In particular, if we denote the pressure on the surface of the cylinder by p, 


we have De—Pn 2 [0 
spU2 — U\ ex} ¢-z, 


Substituting the expression (8) for ¢, we find 


- » 
es te > > x, {eM DSo cos(n+ 1)n—e-—™ +9 cos(n—1)n}, (44 
pU® ho 


n 
where h? = cosh?&,—cos*y, and this gives the general expression for the 
pressure on the surface of the elliptic cylinder whose major axis is parallel 
to the stream. 

We have obtained, likewise, the general expression for the pressure 
distribution on the surface of the elliptic cylinder placed perpendicularly 
to the stream, and this will be given in the Appendix. 

It is of interest to compute the pressure distributions on the surface of 
an elliptic cylinder whose major axis is parallel to the stream. Using the 
above expression (44), computations of the values of the pressure ¢0- 
efficient (p,—p..)/(4pU?) have been carried out for the two cases when 
R = 4:8, b/a = }and R = 3-2, b/a = } respectively. The results are given 
in Table III and are shown in Fig. 4.+ From this figure it is seen that when 
the Reynolds number is small the pressure gradient is inappreciable for the 
most part of the surface of the cylinder except in the neighbourhood of the 
leading edge where a very sharp pressure gradient appears. 


+ In Table III the values of the pressure coefficient at 2/a 1-0 have been omitted, 
since the series on the right-hand side of (44) is slowly convergent in the neighbourhood of 
the trailing edge of the cylinder. In Fig. 4 the curves near the trailing edge have been 
drawn by extrapolation. 
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TapLe III. Values of (p,—p..)/($pU2) 





R=48 | R=3° 


b/a 4 | bla ] 


ze) 3°90 7°67 
o"9 1°72 1°03 
05 0°90 0°33 
06 0°09 or’l4 
o°4 0°37 0°37 
o*2 0°69 0°53 
°) 0°57 0°64 
ool 0°05 
o°d4 0°69 
o'75 0°65 
o7I o"74 
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Pressure distributions on an elliptic cylinder placed parallel 
to a uniform stream. 
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APPENDIX 

For completeness, we shall give here, as an addendum, various results for the cas 
of an elliptic cylinder whose major axis, of length 2a, is at right angles to the uniform 
stream of velocity U. For brevity, only relevant parts of the analysis will be given 
here. 


1. The general solution of Oseen’s equations 
The fundamental solution of Oseen’s equations is given by ¢ and y which satisfy 
equations (5) and (6) respectively. In this case, however, it is convenient to introduce 
the elliptic coordinates (€, 7) defined by 
xtity esinh(é-4 n)s 
where c /(a?—b?). 


The appropriate expression for the function ¢ is given by 
, _ a 
$= cUaré—cU Y 


/ 
Lag 
n=1 


v* e—"§ cosny, 


where the a* are constants of integration.t The appropriate expression for the 
function y is given by 
x 


+ U ekcsinh € cos 9 p3 B* FEK*#(€)ce*(7), 


m=0 


where k U/2v and the B* are also constants of integration. Here, the function 


ce*(7) is the Mathieu function, namely 


y 
: 2 —_ 
ce om(7)) = a cos 2rn, 
r=0 
x 
CeIn +1(7) > AQt si" cos(2r+ 1)n, 
r=0 


m 


while the function FEK*(&) is the modified Mathieu function associated with ce*(y) 
and is defined as 


l S , 
FEK#,(é) = — 2 ALQ™K, (Ae®)I,(4ne-£), 


I * eet 
FEKS, . (é) SS Ag DK, (hucef) Ident) + K,( de®), (4uce-€)}.- 
cat) 


r=0 


The constants of integration a* and £* in the expressions for ¢ and x can be 
determined by the boundary conditions at the surface of the cylinder which is defined 
by € = &. Thus we have 


£ 


1 —~ 
wf = —5 S ce*(0)B%, 
2e Ly 
m=0 
x , . 
x* e-"ho— YF B* F* (€) (cosh £5 (n A) 
n a Pm* mn\S me 
m=0 \0 (n = 2, 3.,...), 


+ The various quantities in this case are distinguished from those previously given by 
adding an asterisk. 
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{—sinh€, (n 1), 
10 (n a oe 


uf ~ Qe ye (¢ 
* n&o = p= GF ee) 


m=0 


l 
Ft Ae {(FEK*’(é)7* ,(€)—FEK#(€) 1%’, (€)}, 
= 
l 
G* ,(€) FEK*(&){« sinh €1* ,_,(€)—nI* ,(€)—« sinh €1% ,,..,(€)}, 
mn 


[zg AP) LS APT .(« sinh €)+T,,,(« sinh €)}, 


2 hk ar =6(*R 


n 


[* (€) L > AQMIYT, op_1(« sinh £)+ J, 0741(« sinh €)}; 


as before, x stands for ke. 
Further, the 8* are determined by solving the following system of simultaneous 
lgebraic linear equations: 
y £0 7 
Y Bx F* (é,)—G* (Eo)} f—¢ (n L), 
— scsi \0 (n = 2, 3,...). 
2. The general expression for the drag coefficient 


By introducing the drag coefficient Cp D/(pU*.2a), we get 


x 
I~ <<“ 
27 aia = 
( D R ce m9) m? 
m=0 
where R = 2aU /v is the Reynolds number. 


3. The general expression for the stream function 
The stream function, expressed in the non-dimensional form, is given by yf, = %/(Ua), 
wher 
sech &€)(cosh € — e6-€ cosh & )sin 7 4 
i Pears ¢ a sin ny 
sech £4 > > Brit Ser f m.n($o)} —— . 


mt m,n 
? 1 m=0 n 


4, The general expression for the pressure coefficient 
The pressure distribution on the surface of the cylinder can be computed from 


the equation 


: S *{e—("-1)€o cos(n + 1)n+¢ "+Nfo cos(n- 1)y}; 


wher h** = sinh?€)+cos?y. 
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THE TRANSIENT MOTION OF SOUND WAVES 
IN TUBES 


By J. D. PEARSON 
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SUMMARY 


This paper falls into two parts; sections 2-5 deal with the transient motion of 
sound waves in an infinite tube of constant cross-section ; while in sections 6—9 the 
effect of slight distortion of the cross-section is considered. In the first part a formula 
s obtained for the velocity potential by use of Heaviside’s operators, and an asymp- 

‘tic expression for this is then found by the method of steepest descent. In the 
second part of the paper a series solution is obtained for the velocity potential in 
a tube of slightly variable cross-section. The particular case of a tube of near 
rectangular cross-section is investigated, and for such a tube a longitudinal harmonic 


listortion is considered in detail. 


1, Introduction 

THE motion of sound waves along an infinitely long straight tube of constant 
cross-section, filled with a perfect gas, when the motion at every point is 
assumed to have a harmonic time dependence, is well known and was con- 
sidered by Lord Rayleigh (1). 

A slightly different problem is considered below. The gas is taken to be 
initially at rest in a state of equilibrium, and a known disturbance is assumed 
to commence at a particular cross-section of the tube. Linear approxima- 
tions are used, the velocity potential is assumed to satisfy the wave equation 
and is taken to be continuous and have continuous derivatives everywhere; 
thus discontinuities such as shock waves are not considered. The velocity 
potential may be expanded in a series of eigenfunctions connected with the 
cross-sectional shape. Due to the linearity of the wave equation each 
such term in the expansion of the velocity potential may be considered 
separately. 

An expression for the velocity potential valid at all points in the tube is 
obtained by the use of Heaviside’s operational method. It contains an 
integral, however, which is difficult to evaluate. An asymptotic form of the 
velocity potential valid for large values of the time ¢ is obtained from the 
operational form of the solution by use of the method of steepest descent. 
A harmonic time dependence of the forcing motion is considered in some 
detail. 


[Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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For this case the resulting field can be split into two parts, namely, the 
transient motion, and the steady state solution or main signal which js 
propagated or attenuated according as the frequency of the forcing field 
is greater or less than the critical frequency of the mode considered, 
The transient disturbance travels with the free space velocity of sound, 
while the main signal rravels with the group velocity of the mode con- 
sidered; this group velocity is less than the free space velocity. 

For a propagated mode the method of steepest descent due to Debye 
‘annot be used at sections of the tube near the front of the main signal, 
since at such points a pole of the integrand is near the saddle point. The 
method developed by Clemmow (2) can be used, however, and a solution 
obtained in terms of Fresnel integrals. 

In the second half of the paper the effect of small variations in cross- 
section of the tube is considered. It is assumed that this distortion of the 
tube modifies only slightly the results obtained from the corresponding 
tube of constant cross-section. If terms of the second order of the variation 
in cross-section are neglected, the problem becomes that of determining a 
solution of the wave equation with given normal derivatives on a tube of 
constant cross-section. This is solved by means of a Green’s function. 

The prescribed disturbance at a particular cross-section may be split into 
a sum of modes. The resulting sound wave will also consist of a sum of 
modes, but energy will be continually transformed between the modes with 
passage along the tube. The modes are no longer independent but are 
coupled to each other by the distortion of the tube. The transfer of energy 
to any particular mode is assumed to be small. If the prescribed disturbance 
consists of a single mode, the effect of the variation in the section is to 
transfer some of the energy from that mode to other modes. Energy in any 
one mode, on passing through a distorted section of the tube, is partially 
transformed to other modes, which, if they have a cut-off frequency below 
the frequency of the disturbance, will be propagated with different group 
velocities from the original mode, or if their cut-off frequency is higher than 
the frequency of the disturbance, they will be attenuated. 

If the variation in cross-section is assumed to be harmonic, resonance 
phenomena occur for certain values of the frequency of the original dis- 
turbance. For such frequencies a large amount of energy is transferred 
from the original mode to another mode, thus violating an assumption on 
which the theory is based. This difficulty does not arise, however, if signals 
of only small duration are considered, in which case the energy-transfer to 
any single mode does not become large. 

I wish to thank Mr. D. 8. Jones and Mr. E. Wild for their very helpful 
advice and suggestions. 
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2. Integral expression for the velocity potential 

Consider an infinitely long straight tube of constant cross-section with 
rigid walls. The tube is filled with a perfect gas with constant specific heats. 
Any motion of the particles of the gas will be taken to be small, and the 
accompanying thermal changes to be adiabatic. 

Let x, y, z be a right-handed system of cartesian coordinates such that 
the x-axis is parallel to the wall of the tube. The wall of the tube is denoted 
by c, anormal cross-section of the tube by A, the bounding curve of A by C, 
and the unit outward-normal to C by v. 

The gas is initially in a uniform state of equilibrium, physical properties 
in this state being denoted by the suffix 0. At timet = 0 the gas is disturbed 
over the plane face of a normal cross-section of the tube, which will be taken 
to be the section x = 0, and it is required to determine the subsequent 
motion. 

If terms of the second order of small quantities are neglected, the velocity 


potential ¢ satisfies the equation 
V2p—— 0, (1) 


where a is the free space velocity of sound in the gas. Having regard to the 
above remarks on the origin of the time scale, ¢ and @¢/ét are taken to be 
zero for t 0. 
The boundary condition at a rigid wall for a perfect gas requires that the 
normal velocity shall be zero; hence 
Cd 


CV 


bo 
— 


0 one. (: 


A disturbance at a normal cross-section is completely prescribed if ¢ is 
given at points on that cross-section. Thus a solution of (1) is required 
with given values of ¢ at x = 0 subject to condition (2) and a radiation 
condition at infinity. 

The operator @/ét will be denoted by p, the Heaviside operational repre- 
sentation of a function of f(t) will be denoted by f(p), and the symbol 
will be taken to read as ‘the operational representation of’. 

If ¢ is prescribed at x = 0 for all positive values of the time, being zero 
fort < 0, subject to ¢ at x = 0 satisfying the condition on C stated above, 
then for points inside and on the curve C it is possible to expand ¢ in the 


().,~0 > a, (t)bn(y, 2), (3) 
where a,,(t) [ (d)..~0 b*(y, z) ads. (4) 


A 
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The set of eigenfunctions y,,(y,z) are solutions of the equation 
no c 2 
é é 
“bn + “bn k24, = 90 over A, 
dy? © ez? 
with boundary condition 


Opn 0 


on C, 


and %(y, z) is the complex conjugate of y,,(y,z). It can be shown that 


( Pn wn dS = 0 (n = m), 


A 
and we choose y, such that 
» cs 
| ly,,|? dS == i, 
A 
Equation (1) becomes, in operational form, 
p* 
V°*h wp 3P 2 0, (5) 
a? 
a solution of which may be written 


$ = > A,(x, p)b,(y, 2), (6) 


and the A,’s are to be determined. 


On multiplying (5) by %*(y,z) and integrating over the area A it is seen 
that 29 7 9 
c2A « , Tk = 
= — — (i + I ) A n — 0. 
Ox" wy 


The general solution of the last equation is 


A,,(x, p) = a,,(p)e Ane 4 B,,(p)ern®, (7) 
where 2 = a + 2, 
a2 


Considered as a function of p, A,, exists on a Riemann surface of two sheets 
with branch points at p = +iak, (points B and B’ of Fig. 1) and branch 
line joining these two points along the imaginary axis. We take A,(p) in 
(7) to lie on the sheet of the Riemann surface on which X,, is real and positive 
for p = 0, this sheet being defined as the first sheet of the Riemann surface. 

The first term on the right-hand side of (7) represents a wave travelling 
in the positive x-direction and the second a wave travelling in the negative 
x-direction. Thus since waves must diverge from their source, we must have 


X,(p) = 90 forx < 0, 


B,(p) =0Q forz > 0, 


while, at x = 0, A,, = G,. 
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It follows that, for x > 0, 


$ = > G,(p)e**,(y, 2), (8) 
and, he 2 < 0, $ = LA, (pem*yfn(y, 2). 


n 
Now it is known that 
x x J, {k,,./(a*t? —2x?)} 
pe Ant — §{t |—H t——jzk,. nvi ——~ My (9) 
a a a (t°—x*/a*) 
where J,(x) is the Bessel function of order 1, 5(t—2/a) the Dirac delta 
function, and H(t—2x/a) the Heaviside unit function. 
Hence, by the operational representation of a Faltung integral, it is seen 
that 
f 
dy 
v(y?—2?/a?) j 


(10) 


t 
$=) H(t a) (y,2)\a (t- a} } An(t—y)ky Stn V(y?a*—2*}} gy| 
aw iat Bias J 


a 


ria 


for z > 0, and a corresponding expression for ¢ exists for x < 0. 
The integrals in the above expansion for ¢ are difficult to evaluate for 
most forms of a,,(¢), but an asymptotic form for ¢ may be obtained by the 


method of steepest descent. 


3. Asymptotic form of the velocity potential 

Since the modes are propagated independently of each other, only the 
propagation of a single mode need be considered, and thus at x = 0, ¢ is 
taken to be equal to a, (t)y,,(y, z). 

In order to obtain an expression for ¢ valid for large t, the following 


integral will be considered, 
l a, (p) 
I exp(pt—aA, 2) dp, (11) 
a Pp 


where d > 0,x > 0, and @,,(p)/p is regular on the Riemann surface of d,,(p), 
and the path of integration lies on the first sheet of this surface. 

By putting 2 = cta, where 0 <c < 1, and regarding c as constant, (11) may 
be written 


d+iza 
l 4 a,(p) 
awn! 1-2 9 
I — expitg(p)} dp, (12) 
aTl . P 
d 
where q( ~P) p ca, (p). 

Cases for which ¢ > | are not considered, since they correspond to points 


ahead of the disturbance. The integral (12) is in the form for which an 
asymptotic expression can be found by the method of steepest descent. 
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It is easily shown that g(p) has saddle points at 
p = +tak,(1—c?)-* 
(points S and S’ of Fig. 1) on the first sheet of the Riemann surface of r,(p). 


The path of integration of J is deformed to pass through the saddle point 
and to lie along the line 


im(p—cad,,) = constant. 


y 


RI 


—> 
P 





IMp 











Sheet 1 Sheet 2 
Fia. 1. 


A diagram of the path on the two sheets of the Riemann surface is given 
in Fig. 1. 

Applying the method of steepest descent, contributions are considered 
only from the portion of the integrand near to the saddle point, and Debye’s 
formula gives for large t 


> Bos 27 + (k,, a)te .” iak, eitk, av yi ag —tak,, e—itk, av ot) | 
" t tak, (1 ~ c?)é 7 (1—c?)4 . . (1—c?)# 


(13) 
The expression (13) is valid only for large values of \tg”{iak,, (1—c?)-*}| and 
thus cannot be used at a zero of g"{iak,(1—c?)-4}, where accents denote 
derivatives. The asymptotic form becomes invalid, therefore, for values 
of c near c = 1, i.e. for points near the wave front. 
[f the original disturbance is taken to be such that at x = 0 


$=a,(t)p,(y,2) for t > 0, 
then d = d,,(p)e»*p,, (y, z) (x > 0), 
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which gives on inversion 
d+i« 
a, 
271 L 


d—i« 


) 
d ws, (y, z) f ) otovw) dp, (14) 
which is of the form of the integral (11). 
In order to apply the formula (13) it is necessary that a,(p)/p is regular 
on the Riemann surface of A,,(p). If a,(p)/p has poles, then contributions 


from these must be considered. 


4, Harmonic forcing motion 
Now consider the forcing term 


a,,(t) H(t)cos wt, 


then a,(p) = — _ 
er 
Substitution of this expression for a,(p)/p in (14) gives rise to an integral 
which is formally similar to that of expression (11). The integrand, however, 
has poles at p Liw. On deforming the path of integration to lie along 
the line of steepest descent, contributions to the integral may arise from 
these poles. As is seen in Fig. 1, this is the case if 
w k, a(1—c?)-#, 
which implies w > k,, a and t > xw/a(w*—a?k?)* and a contribution 
x W x ‘ P P 
Hit ———=5 }Pn(y,2)e0s wt ——(w?—k? a?)}) 
\ @4/(w*—a*k;) a 
must be included in our asymptotic form of (14) as well as the term given 
by (13). It is easily shown that the path of steepest descent crosses the 


branch line of the Riemann surface at the points given by 


p +-itak, (1—ce?)# 
(the points P and P’ of Fig. 1). Hence if w < k, a(1—c?)?, which implies 
that w < k,a and t > 2k, (k2 a®—w?), a contribution 
vk. F x 2.9 2\h 
Hit . —— |b (y, z)exp (k° a*—w?*)*) cos wt 
a(k= w*/a~)? a 


arises from the poles at p |-iw, for x > 0. 
The asymptotic solution for ¢ such that 
d H(t)cos wtb, (y,z) atx=0 
is for n 0 given by 
» 


‘ l r »je2 2 eosltk __p2)t 17} i 
as (== Pal x fe cos{tk, a(1—c*)*—} 4 lt (15) 


tat a} (1—c?)k? a2—w?(1—c?)} 














320 J. D. PEARSON 


for x > 0, or 


Pre l—-( 2 }Pa(t—2) atak? cos{k, (t?a? —x?)! — 47} mi $0 


mk, a a) (t2—a?/a?)#t?(w? —k2 a?) —2?w?/a?} | 
x PT oa 
where ¢@ = H\t——_—~ ne —— |cos wt —~ (w2—k? a2)! 
a(w?—a?k?)! a 
if w > k,a and 
xk Fe 
’ o Bt . exp| ——(k? a2—w?)!) cos wt 
( a(k2—w? a) i _ ' 


ifw <k,a. For the case w > k,,a, ¢’ will be called the main signal term, 
If n = 0 contributions to ¢ arise only from the poles at p = +iw to give 


$= al — “to cos( — =I. 
a a 


Thus, a disturbance which is constant over the cross-section is propagated 
unchanged. 

The velocity of the disturbance represented by ¢’ is seen to be the group 
velocity of the nth mode. 

The transient disturbance represented by the first term on the right-hand 
side of (15) takes infinite values when c = 1 and when c = (1—a®k2/w)t, 
which correspond to cross-sections at the front of the disturbance and at 
the front of the main signal, respectively. Thus in the neighbourhood of 
these sections the asymptotic form (15) becomes invalid. For values of ¢ 
near the front of the main signal, the pole at p = iw and the saddle point 
at p = iak,(1—c?)-! are close to one another and (15) tends to infinity. 

At a point near the front of the main signal, Debye’s formula may be 
applied if ¢ is taken sufficiently large, but the formula breaks down at the 
actual main signal front. An asymptotic form which is finite and con- 
tinuous at the main signal front is obtained in section 5. The transient is 
seen to oscillate about zero with decreasing amplitude and increasing wave- 
length, for fixed t, as 2 tends to zero, or for fixed x, as t tends to infinity. 

For general prescribed value of ¢ across a normal cross-section containing 
all the modes, with a time dependence H(t)cos wt, it is seen that all modes 
are attenuated except the modes with critical frequency smaller than w. 


5. Approximation near the front of the main signal 
In order to obtain an asymptotic expression for ¢ which is continuous at 
the front of the main signal, the following integral will be considered, 


I = Ls — P expitg(p)} dp, (16) 
2 pe 
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shere w > k,a. This integral can be written 
d d 
1 l i | exprig( p)} dp | explg(P); 
47ra\ . p—w J p+tw 
d d 
* expilt , xpit , 
eXPUMP)S dy | € PUMPS gy ; (17) 
J P+ tw P. P tw 
d ix 


The third and fourth integrals in (17) have integrals which are regular on 
he Riemann surface of A,,(p) in the range of integration. Debye’s method 
f steepest descent can be applied to these integrals and their sum is 
symptotically equal to 

“| cos{tk,, a(1—c?)?— 47}. (18) 
nt} (1—c?)Nk, a+w(1—c?)?*} ™ 

The first and second integrals in equation (17) will be evaluated by an 
mended form of the method of steepest descent. 

Consider the integral . 
* expitg(p)} 


P iW 


dp. (19) 
The path of integration is deformed to lie along the line of steepest descent, 
the contribut ion from the pole at P iw being dealt with as in section 4, 
It can be shown by methods similar to those used in the proof of Watson’s 
lemma, that, neglecting contributions from the pole, J, defined by equation 
(19), is asy mptotically equal to 
. hf? lz 
( . de 
I, — et] — = a (20) 
: [—¥ (Po)}*(Po iw) 


where py ik, a(1—c?) 


The function D is defined to be equal to 


(| py| —w)(1—c?)!/e(k,, a)? 
and thus (py—iw){—g"(po)}3 i Detiz 
Equation (20) may be rewritten 
P —"* 
L~ - pa CXPUI (Po) +- Laar— diz?) dz. 
Ni 2 27)) 
NOW = pe | exp{—a(z?+ D*1)} dx 


U0 






ft 9 
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and thus ao 
I~ | | iD exp{tg( py) —4tz?—a(z2+ D®i) + lim} dadz. 
—« 0 


In the analysis below the order of integration is reversed and some justifica- 
tion of this is required. The order of integration may be reversed if the 
double integral of the modulus of the integrand is convergent. It can easily 
be shown that the double integral 


ies 
| | |Dexp{tg(p)—}tz?—a(z?+ D%i)+ hiz}| dadz 
—x 0 


is convergent for finite X. The integral 
x x 

| exp{— 4}tz?—a(z?+ D%i)} dadz 
x xX 


R(X) 





will now be considered. Carrying out the integration with respect to 2 in 
R(X), it is easily seen that 


v4 


R(X) = ( exp{ — }tz?— X(z?+ D%)}— =; dz 
} 2 


yy2 1 D2 
€ € x " X 
. : i - lexps tz2?— X (22+ D7 
or R(X) = | | | =— ——_ J Bde 
a e e 2+ D% 
t. € € 
where « is positive. 
Now 
a ¢ if N(2?4D?i) l a 
| dz | e—t+ X)z* dz 
J 24D% D2 | 
€ € 
_ 2e 
- pe 
and 
*6 x 
ha ¢ his X(22+.D?7) l fe 
es dz “ € (4+ X)z? dz 
P 22 -_ DP )2 
€ € 
x 
e l X 2 ¥ 4422 
= € € e dt dz 
D? ; 
0 


We choose e€ such that Xe? is large, say « = X-+, and thus R(X) = o(1) as 
X—>o. Thus 








o « a « 
e~He*—a(2*+D*) dadz = | | e-te*-ae*+D%) dzdxz-o(1) 
-x 0 0 
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and taking the limit as X — oo we have 
IL~ | | iDexp{tg(po)+ }ia—}tz?—2(2?+ D*i)} dzdz. 


0 x 


Carrying out the integration with respect to z in the last integral, we have 
‘ - ) 


; it D 
im 
xD fs . 
L~ | exp{tg(py)+ jim—ixD?} Gan! dx. 
0 a 
; y [z\4 
Putting (x+ it)? ” 
i a) 


the last integral becomes 


ir D 24 


j , sae es 
L~ | exp{tg(po)+ jia— himy?— }uD*} ——__ dy. (21) 
; . a S | D| 
D\ J (t/7 
The Fresnel integral &(x) is defined to be 
é (x) | exp(}y?ai) dy (22) 


I 


and &*(x) is the complex conjugate of &(2). Equation (21) becomes, on 


) 


using (22) and including the contribution from the pole, which has been 


omitted in the above analysis, 


[, ~ 227i exp{iwt —(ix/a)(w?—k? a?) 


1 


i] 


2'ri& *{| D|(t/7)expitg(po)+ }ixa—3uD?} (23) 


—— 


nw 


for w Po| and 
I, ~ 2'mié*{| D| (t/7)}expttg( po) + fia— 3tiD?} 
for WwW Po}- 

The transformation from equation (20) to the form (23) may also be 
carried out by the use of the definition of the confluent hypergeometric 
function (2). 

By a method similar to the above it can be shown that the integral J, 


d 
, ms expilg(P)} dp, 
Pp tT bW 
has the asymptotic form 
I,~ 27 exp{ —iwt + (ia/a)(w?—k? a?)#}— 
for w Po and 
I, ~ 2mi&{| D\ (t/7)}exp{—tg( po) — fia $itD*} 


for W « Po . 
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A few properties of the Fresnel integral defined by (22) will now be 


given (4). 


[It can be shown that 


6 (2) { B(x) +-i1A(ax)redriz®, 


a i » 
where B(x) | e-inaty_Y dy 
72? | l+y? ° 
0 
t £ y-* 
and A(x) | e—tra*y _: dy 
»s 9 J* 
aD? : ] J o 


For large values of x 


We define the angle y such that 
x(u) = arg! B(w)+7A(u)'. 
Tables of A(w) and B(u) are given by Rankin (4). 


ii at z= 0 d H(t)cos wt b,(y,z) (w > k, a) 


then for x > 0 db wb (y,2)L, 


where J is the integral (16). Thus 
d~ In(ys2),2 1! {| D| (t/z)!\|cos| tk, a(1—c?)! x{|D|(t/7)}+ 4a)4 


1 (k,,a)'e cos{tk, a(1—c?)} 
¢ " 2\4 9\] 
V2nt {k, a+w(1—c?)!(1—c?)t | 
for w < k,a(1—c?)-+, and 


db~ daly 2)) 2-*|&{| D| (t/z)*\\cos|tk,, a(1—c?)! —x{|D|(t/7)#}+-42]+ 


1 (k,,a)'c cos{tk, a(1—c?)!—17} 


. . cos! wt — (a/a)(w?—k? a2) (25) 
\2nt {k,a+a(1—c?)!(1—c?) ; ee 


for w > k, a(1—c?)-4. 

It should be noted that at the front of the main signal expressions (24) and 
(25) are equal. Thus the asymptotic expression obtained above is con- 
tinuous at the front of the main signal: it is therefore of greater value than 
the expression (15) obtained previously, for values of c corresponding to 
sections near the front of the main signal. 


For large t \&{| D\(t/z)8}| ~ 
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and thus the expression (24) and the first two terms of expression (25) be- 


come for large ¢ equal to 


2k, a\4 chia ; 
Pri¥s 2)| ; }° 2 » = 9 9  cos{th, a(1 c?)t— jo}. 
sis mt ‘(k a)?—w?(1—c?)}(1—c?)! 


? 


Thus the expressions (24), (25) reduce to the values of ¢ obtained by Debye’s 
method of steepest descent if ¢ is taken sufficiently large. 


P/V 


of 


Enve/ope 








A graph of the amplitude of ¢/%,,(y,z) is given in Fig. 2 for particular 
values of ¢, k,, and w. For points far from the main signal front the ampli- 
tude of d/yus(y,z) was calculated from (15), and for points near the main 
signal front expressions (24), (25) were used. The values of t, k,, and » 


used were 


t 5007 /a see. k 0-1 em.~! w 4,000 cycles/sec. 


6. Tube of variable cross-section 
The propagation of sound waves along a tube of slightly variable cross- 
section will now be considered. As before, ¢ is prescribed at x 0 for all 
0, and the form of the resulting sound wave is to be found. 
The tube of constant cross-section, the propagation in which has pre- 
viously been considered, will be described as tube 1, and a tube of variable 


cross-section whose cross-section for any 2 differs but slightly from that of 
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tube 1 will be denoted as tube 2. It is assumed that the cross-sections of 
tubes 1 and 2 are the same at x = 0. The velocity potential in tube | 
will be denoted by ¢,, and in tube 2 by ¢,; ¢, and ¢, are assumed to differ 
but slightly from one another. If the mathematical solutions obtained 
do not agree with this assumption, the method will automatically be 
invalid. 

The bounding curve of tube 1 is assumed to be continuous and possess 
piece-wise continuous first derivatives. It is assumed to consist of N con- 
tinuous curves whose equations may be written 


fiy.z) = 0 ( = 1, 2....,.¥). 


It is assumed that the variation in section is such that the bounding 
curve of cross-sections of tube 2 consists of NV continuous curves whose 
equation may be written 


fy, z)+eR (x,y,z) = 0 (r oy See (26) 
or Sty t+e¥ (x,y,z), 2+¢Z,(x,y,z)} = 90 (r= 1, 2,..., N), 27 


where ¢ is small. It is assumed that F,, Y, 


, 4, and their derivatives are 
bounded. Thus equations (26), (27) are equivalent when terms of O(e) 
are neglected, if 


* 
y 02 


Fi(z,9,2) = Y, oh Z Ch (r = I, 2,..., NW). (28 
Cc 


Terms of O(c?) will be neglected throughout and the value of ¢ accurate to 
O(e) will be determined; ¢, and ¢, are taken to be equal at 2 = 0 for all 
positive time, and it is assumed that 


d, = +9, (29) 


where e® is O(e) and terms of O(e?) are neglected. 
The boundary condition on the wall of tube 2 is that the normal derivative 
of ¢, shall be zero. Ifl,, m 


to the rth section of the bounding curve of a normal cross-section, then from 


n, are the direction cosines of the outward normal 


= 
equation (26) it is seen that 


oF 
«e—, m 


oF. éf. oF. af. 
= €- ra) , — €—+—. 
Cx cy cy Cz Oz 


_ 


r 


(30) 


The boundary condition on the wall of tube 2 may be written 


Obs | Od, Od. 
l, —=+m, —+n, — = 


¢ Ys 
OX oy Cz 
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Substituting from equation (30) in the last equation, we have 


Of, Obs Of, Obs . ate ate a = (31) 


Oy Cy Oz Gz x Ox Oy OY G2 @& 


Thus from equations (29), (31) the boundary condition on the curve whose 
equation is f.(y,z)+«RM(y,2,z) = Ois 


of. Od, Of, Od; f F’ 0d, OF, é¢, oF, é¢4, of, oD ; aad 0 
am -€ + -+ i + at = 


Cy Cy cz Cz Ox Ox Gy Cy Gz Gz ody OY Gz & 


(r = 1,2,...,N). (32) 


Expanding @f,/éy, Of,/0z, 0b,/ey, 6d, /éz in Taylor series and neglecting terms 
f O(c?) the boundary condition (32) may be expressed as a boundary 
condition on f.(y,z) = 0. Thus on f(y,z) = 0, having regard to the fact 


that 
Of. 0d, , Of, Ody 


CY CY Cz C2 


on this curve, it is seen that 


eM (: L Y, oy, Z,| Od, | (: Sr Zz. | : 4 r) od, 
CV, 


(bay | Oh 2) Ala, Pay) 
2 : dz?" ° Gydz * 
oF. ed, oF, cd, _ OF, 64, 


Ox Ox cy cy C2 Oz 


, (33) 


where v, is the outward normal to the curve f(y, z) = 0. 

The problem of finding a solution of the wave equation with zero normal 
derivatives on a tube of variable cross-section has thus been transformed 
into the problem of determining a solution with given normal derivatives 
on a tube of constant cross-section. 

It is seen from equation (29) that since d, = ¢,atx = 0,® = Oatx = 0; 
then ® is a solution of the equation 

yo... & (34) 
a> ct? 
neglecting terms of O(e?). 

A solution of equation (34) is therefore required with initial conditions 
that ®, @@/et are zero at t = 0, and boundary conditions ® zero at x = 0 
and ¢@/év given on the walls of the tube. There is also the condition that 


® is zero everywhere for t < 0. 
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The application of Heaviside’s operational methods to equation (34) gives 


vo? & = 0. (35) 


a~ 


At points inside the curve C, a solution of (35) may be expanded in the form 


D = DE,(X, Palys 2), (36) 
where Cy (2%, Pp) : @ b*(y, z) ds. 
A 
The boundary condition on ® at x = 0 implies that ¢,(a, p) is zero at 


x 0, 
Multiplying equation (35) by *(y,z) and integrating the resulting 
equation over the area A, it can be shown, by the use of Green’s theorem, 


that “ae 2 ‘ame 

oC es oP o® , a 

—_ (is 4 = Cn : or dS = U, (37) 

Cx a~ a Cv 

A solution of this equation will be found by the use of a Green’s function. 
The Green’s function G,,(x, x’, p) is defined such that 

eg 

? 


$A? G,, d(a—z2’), (38) 


» 
Cx* 

where 8(a—z’) is the Dirac delta function, 
G.—>0 as zxr—>- 00, 


and G, = 0 at x = 0, It is easily shown that 
G le A,|e—2"| __ pA, a+ | forx > 0,2’ > 0. (39) 


Consider the function 


C_(x’, a n (x’,p)G, dx’ forzx>0 (40) 
Be (a’)? n\® > P) ¢ | O(a’)? vs PIGr 
0 0 


a] 977 . 9- 
Cc , oc“ 


which is seen to be zero on integrating by parts. Thus substituting from 
equation (37) in the expression (40) we have 


"7 e@ 
€,(x, p) | ) ' be ds-G,(a,x',p) da’ for x > 0, (41) 
CV 
0c 


and a corresponding solution exists for x < 0. 

Thus a solution of equation (35) has been obtained as a sum of eigen- 
functions, and by inversion of the operational form an expression for ® 
an be obtained. 


It is seen from equation (41) that the change in the velocity potential 
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caused by a distortion of part of the boundary curve C is independent of 
the distortion of the rest of C if terms of O(c?) are neglected. 


7. The effect of distortions in a rectangular tube 

Little advantage is to be gained by proceeding with the analysis for a 
tube of general cross-section any further, and a particular example will 
now be given. 

A tube which at « = 0 has a rectangular cross-section, with sides of 
length b and d, will be considered. Since the contributions to ® from each 
of the four walls have been shown to be independent, a distortion of one 
wall only need be considered. 


Let the equations of the walls of the tube be 
y= eF(z,z), y=6, 2-=0, and <=—d, (42) 


where F and its derivatives are bounded. 


From equation (33) the boundary conditions on the tube walls are 


e® 2, oF e®, @éF e® 
1 F(z, z) ann wf — ony= 0, (43) 
cy oy~ COX CX Cz C2 
e®/ev 0 every where else on C. 


A further simplification is now introduced, that the cross-section is 
assumed to remain rectangular everywhere. Thus F becomes a function 
of x only, and the third term on the right-hand side of equation (43) becomes 
zero. The function F will be taken to be equal to sin ax, where «x is such 
that ae is O(c). A more general form of F(x) is considered briefly in section 9. 

If, at x = 0, d, consists of a sum of modes, energy will be continually 
transferred between the modes as the sound waves move along the tube. 
The case when, at x = 0, ¢, consists of a single mode is given below. If this 
is regarded as the consideration of a single term of a more general form of 
dat « = 0, then the analysis given can be regarded as the energy lost in 
one mode neglecting the energy gained by that mode from the other modes. 
\n expression for this neglected energy can easily be obtained, however, by 
a summation over all the modes of the energy converted into the mode 


under consideration. 


Let db a (p)eb, (Y¥,z) atz 0), (44) 
Z 7 8712 
where us cos - J cos Pe Pas 
(hd) b d 
where P, | (7 Q), Po 2 + 
>» >» g2 » 
2 fa #2" 
and hrs " » 
b? d? 


For a rectangular tube the summation parameter n of the tube of general 
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section is replaced by the two parameters r and s and a double summation 
is required, and thus for a rectangular tube equation (36) becomes 


D = er (x, p)bs,.(Y, 2). 


It has been shown that for x > 0 with the value of dat x= 0 given by 
equation (44) we have 


-, = P)br (Ys z)eAne, (45) 


9 
° 

where 2, = i, a=... 
a 


a“ 


Substituting from equation (45) in equation (43), the boundary condition 
on y = 0 can be written 


(: 7} = [¥,.(Y, 2) ],-0e 'wedn(p)| (FF) sin XX — x COS aX re - (46) 
cy y=0 ; b 


It is easily shown that 


—hC/oaxe 
. e® 2 = ‘[rzx\?.. ) 
— yr, dz = Pr Pn Sem 4 6A, P) sin ax— a COS a A, | , 
cy b b 

:" y=0 

0 

(47) 
where 6,,, is the Kronecker delta function. 


Substituting from (47) in the general expression (41) and noting that 





e® o® 
ey ov 
at y = 0, it is seen that 
- 9 7 rrr 2 
_ ie -A,,2’ | 8 ms (ple 
Cum “,p) = —_ ——_ Pr Pn 7 pars rT 1 SIN wa ( — XCOS ax ‘r, .|@ nm dx 
0 . . (48) 
for x > 0, where 
] 
ad -/ 7 a : — Te ‘ 
Gam(%, x’, p) = rv -[e Nnm|C—2'|__ g—Aym|ar tar ]. 


nm 


In order to carry out the integration in equation (48) the values of two 
integrals will now be given. It can be shown that 


y Ant’ os _/ -! 
| Game sin aa’ da 


l A, a x24 f2 —k2 a2 . 
—_— (e Ave? e App cos x2) ak ( — == ra) sin a2 ors 2 ; 
x( p? + 2 mn) és x*( p* T Ynm) 


(49) 
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and 
( asia AX’ eos ax’ da’ 
0 
l r..a? : 24 (2 [2 \q2 { 
hs . sin az é re (OF ‘ __ ra) e (¢ Anm2__ eos x € Ars) . 
2|0(p"+Ynm) 2a°(p*+Ynm) 
(50) 
where y,,,, is defined by the equation 
2 a2 x" | k*. T Kum (ke n- =) 
fnm 4 i 2 | 4q2 
It should be noted that 
Yom 2 l 9 9 9 z - 
o Kom ) (a* +-k2,— Kim) = 0 (51) 
a* | OX 
ae ree we - 
and a= — ke, | (a?+ kh? —k?,) > 0. (52) 
a* 2a 


When square roots of the expressions (51), (52) are taken it is understood 
that the positive root is to be taken. 
Substituting from equation (49) and (50) in equation (48) it is seen that 








fe) ». p,, a? { ram\? o? +k om 
é (x, p) sm Pr Pn a,.(p) ~_ TT *nmT8__ 2 I gin ar e—Anet4 
nm } > 9 ij rs t 9.2 rs 
p° T Ynm , | 6b aa" 7 
A.,[ra\? Xr, . , , oe 
(e-em —cos ar e-'ve)| eT) 4 (a+ Ky —K8)| (58) 





The operational form of the solution with a general time dependence at 
« = 0 has thus been obtained. A particular form of this dependence will 
now be considered. 

Let (4),-9 = cos wt %,,(y,z) for t > 0, where k,,,, > w/a > k,,. 

in the analysis given below, only the steady state or main signal terms 
of ® will be considered, the transient terms being neglected to simplify 
somewhat the form of the results. The region in which the transient terms 
may be neglected in comparison with the main signal terms will be dis- 
cussed later. For values of © outside this region the transient terms must 
be included and can be calculated by methods similar to those used in 
sections 3—5. 


The functions 7,,,,, J, and e,,, are defined by the equations 


9 9 € 9 
1 km 5 ] ark; $ 
Pnm 2 ? Tn m 2 
w Yn m 


and € 0 ifn - € os J ifn <r. 
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Neglecting the transient terms it is easily shown from equation (53) that 


a a> 1 v\ 
Cum, (X) Onn Pr Pn aba v2 alr H(t Pang xX 
= fnm 


er 2 1 , 
P ‘ . ro. W ‘ é e a Z ‘ P 
K jot+k2, —k2.4+-2 k?.)* sin{ wt —— (w2—k?2,, a2)! 
“N nm rs 9 rs nm 
b a? a 
w 9 I 
=—Kr,]" COS ax 
aa 


© Ho 9» 9 9 4 
 sin| wt —— (w?—k?, a*)t) 
a 





2 2 2 [yae\ 2 9:2 
9 ].2 x? + nm hk rs[ 17 WW 
2KS, 0 - ———a 


x 6 az 


f | = Jo) p2 219 ‘r7\") 
Alt Ynm ‘)| \* Kom ae 2 h 1 x 


9 9 9 \ 
xT ki, mm ire . x 2 2 2\4 
( sin Ynmt os (Yn . oun a“)? T 





sin 0% cos (w? k,a*)!} +. 
a 


2a a 


H(t—dn2) [>t +n He, +2(°2) 
a ; 
x 


; sin(y,.,t— (Ynm ks )” we) es 
\ a J 
In this expression for c,,,,(7) the transient terms have been neglected. 
These terms are small compared with the given expression, for values of f 
such that 


(a?+ 


9 9 
“nm kes) wv 


2a 





t > p,.%/a, t> pr t/a forn <r, t > 9,.%/a, t > dnm/4. 
Now p,,, < p,, ifm < r; hence the second condition is implied by the first. 
Thus if a value of the time 7'(«,a,n) can be found such that 

T(x, x,n) = maximum value of [p,., Zs Inm |t/4, 
then for ¢ T'(«,«,n) the transient terms may be neglected. 

The function q,,,,, becomes infinite when 


+h? —k2 = 0 (« 40), 


"nm 
and the function q,, becomes infinite when 


x22, —k2, = 0. 
Thus for a?-+-(k?,—k?2,,) = 0, T’ becomes infinite and the third and fourth 
terms in expression (54) become zero. The transient terms in this case 
cannot be neglected for any value of the time t. 

In the work that follows ¢ will always be taken greater than 7'(«,7) 
unless otherwise stated. 


From equation (29) it is seen that 4, has been regarded as the sum of ¢, 
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which has been calculated in sections 2—5 and the term e® which has been 


calculated above; «® represents a disturbance which will be called the 
secondary wave and which may be regarded as having its origin in the 
wave ¢, acting on the distortions of the tube walls. 

It is seen from the form of equation (54) that c,,,,(a) is composed of four 


nm 
distinct terms. Each term represents a wave with a velocity which is 
different in general from the velocities of the waves represented by the 
other three terms. 

The wave represented by the first term vanishes except for modes below 
cut off and travels with the group velocity of the nsth mode. It can be 
regarded as representing the accumulation of the secondary waves below 
cut off, on passing down the tube. 

The coefficient of the term in the rsth mode of , i.e. the coefficient of 


the term in the original mode of ¢,, is seen to be c,,(x) from equation (54), 


and equals 
( a? Sa?4+-2(ra/b)**(w? » \t F w o \4 
(p 2) K?.}* (1 COS «x )sIN | wt | = -k. at 
| 2ba Y- w* a~ a“ 





a® {k? l(ra/b)?—w?/a*} . { w* o \} 
: : SIN ax COS| wt | ed 


a® {x*-+-2(rz/b)*} , ’ 
2 2 (SIN(Y,,.¢ I $ WW) sin(y,.¢ Loa)} |. 
1b (Y7s—w*) | 


If a pulse of finite duration is considered such that 


(oi. cos wt! H(t)— H(t T) hb es 
then for points such that ¢ > g,,a%/a-+-7, i.e. for points behind the main 
signal, 
a2 fa2+ 2(ra/b)?}(2 cos(y,.t—4y,.7-+ Fax = 
‘an (p,)? ‘ : if Yrs'— BY re 71 2 d\n Yr 
th (y5,—w*) | 2 COS(Y,44—$Y rT Sax) 2 


Thus behind the main signal there are two waves of equal amplitude 
travelling in opposite directions. Alternatively, the disturbance behind 


the main signal can be regarded as a system of standing waves. 


8. Resonance 


It is seen from equation (54) that a resonance phenomenon occurs in the 


nth mode when w? = y%,,, i.e. when 
2 4 2 2 a2 .2 \2 
w x | pe T a (k nim ks)? 
a 4 2 42 


For this value of w the expression (54) for c,,,,(x) becomes infinite and is 


meaningless. 
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Resonance occurs in the working mode, i.e. the rth, when 


a 
2 
a sl thre 
A graph of the resonant wave-length of corrugation (27/«) of the tube 
against the wave-length of the initial disturbance (27a/w) is given in Fig. 3, 
The initial disturbance is taken to be in the 1, 0 mode and graphs of the 








2b 
2 ys ode 
£77 0.0 
W 
JO mogo 
b 24 2b 3b 4b 
a 
Fic. 3. 


resonant wave-lengths for the first four modes are given. It is seen that 
resonance cannot occur when the wave-length of the disturbance is greater 
than the critical wave-length of the 1, 0 mode. 


It can be shown that, neglecting transient terms, 


y? wa Ls oa a 
—— —_— Aygt — (/ —.—— sin{ wt —— (w?—k?,a?)*}, 
(p?+w?) .  A(w?—k?. a?)t}2w a 

4 . , 

y WwW Wa ° a 9 — ‘ 
=. —eAnt = —=|t— —._—__.y ; SIn{ wt — — (w* k2 a*)t) + 
(p?+w?)? 2 a(w? — k?, a?)! a 

1 x 9 9 2\h 
+ cos| wt — — (w? ~Ke,a)), 
a 
= P A e—Arsh 


] 3 9 I 0 ¢ 9 2 9 
2w\ a? a(w? —k?, a?)! a 


1 fw® ,,\-4. a ‘ad 
+53 — #8 sin (ot = (w*- Kat), (59) 
2a? 


a2 a 


On substituting these values in equation (53) it is seen that the disturbance 
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will remain small in the rth mode for values of ¢ such that 
t me (wheat), 
a 
which is satisfied at the front of the main signal. 

To obtain an expression for ® valid at the front of the main signal in the 
case of resonance, equation (53) must be interpreted by use of equation (55), 
to the right-hand side of which must be added the transient terms which 
can be calculated as in section 5. 

The solution for ® so obtained will become invalid at some distance 
behind the main signal front and will then be comparable with the term 4,. 
For a short signal, however, ® will never be comparable with ¢, and thus 
the solution indicated above will be valid everywhere. 


9, More general variation in section of a rectangular tube 
If the function F(x, z) of equation (42) is taken to be a function of x only 
and of the form : 
F(x) = .(2/m) [ sinawf(a) do (56) 
fr 
then f(x) 4 (2/77) ( sin ax F(x) da. (57) 
0 


Now F(x) and ¢F'/éx must be bounded, and thus the following must be true: 


} sin aa f(a) da = O(1) 
0 
and ( vsin ax f(a) da = O(1). 
0 
These conditions are satisfied if f(a) is of bounded variation and |af(«)| 
integrable from 0 to oo (3, p. 425). 

In order that F(a) shall be expressible in the form (56) it must be odd 
about « = 0. The results of sections 6-8 may be extended to this more 
general form of variation in section by multiplying them by ,/(2/7) f(«) and 
integrating with respect to «a from 0 to oo. The form (56) could have been 
taken as the variation in section throughout the analysis of sections 7 and 8, 
all differentiations, integrations, summations, etc., being then carried on 
under the integral sign of «. The form of f(«) is assumed to be such that the 
change in the order of integration and differentiation, etc., is valid. 
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A NOTE ON GRAVITY WAVES OF FINITE AMPLITUDE 
3y J. N. HUNT (Imperial College, London) 
[Received 7 February 1952] 


SUMMARY 

Some corrections are found necessary in Struik’s paper on travelling waves. Thes 
do not invalidate the proof of the existence of finite travelling waves on water of 
finite depth, but lead to corrected forms of some of the approximate expressions for 
the wave profile and the speed of propagation. 

It is pointed out that the amplitude of highest waves must decrease with decreasing 
depth. The definition of the velocity of propagation is considered, and an expression 
for the velocity to the third order is obtained, which indicates some necessary 
corrections to certain technical reports. 


1. Introduction 

THE purpose of this note is to draw attention to, and to correct, some 
inaccuracies in Struik’s paper (1) on travelling waves in water of finite 
depth, and also in some technical reports (2, 3) which use the results 
obtained in Struik’s paper. The corrections are mainly in the third-order 
terms of the wave profile equation and in the expression for the velocity 
of propagation. A small correction is also required in the expression given 
by Stokes (4) for the velocity of propagation of these waves. A number 
of mistakes in sign also occur in that part of Struik’s paper which establishes 
the actual existence of waves of finite height in water of finite depth. For 
the sake of completeness the necessary corrections are indicated and are 
found not to invalidate the proof. 

Levi-Civita finally established in 1925 the existence of gravity waves of 
finite height in water of infinite depth (5). Struik’s investigation uses an 
extension of Levi-Civita’s method, but it was recently pointed out by 
Bowden (6) that Struik’s solution does not reduce to that of Levi-Civita, 
nor exactly to the solution obtained by Stokes, when the depth of water 
tends to infinity. When, however, the above corrections are made, it is 
found that all the solutions are consistent. 


2. Corrections to Struik’s proof of the existence of travelling 
waves in water of finite depth 
In Struik’s induction proof of the general form of the solution, several 
errors in sign occur, principally in equations (42), (54), (60)-(63) of Part I 
of his paper and in section 10, which deals with the convergence of the 


[Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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method. It will be convenient to summarize briefly those parts of Struik’s 
paper which need altering; the equations in sections 2—4 will be numbered 
as in his paper, those equations which are different from his being indicated 
by a dash. 

Following Levi-Civita and Struik, cartesian axes are taken with the axis 
of y vertically upwards with ¢ = % = 0 at the origin and % = qg on the 
free surface, 6 and y& being the velocity potential and stream function 
respectively. The complex velocity is w = u—iv. The region 0 < % <q 
is transformed into the annulus 1 < |{| < R of the ¢-plane, where 
(= erifed, R = e®r7acd — e% and f = 6+ ip; A is the wave-length and c the 
wave velocity. 

Levi-Civita (7) has shown that the free surface condition leads to the 
mixed equation 


d Ac ] ] + 
C —.{w(e-*C)w(e%Z)} _g cf _ | a (E) 
dl 2m |w(e-*f)  w(e%L)} 
to be solved on either boundary |f R or {¢ R- of the annulus [, 
and where w is a holomorphic function of ¢ in I’, which is real on |f l. 


Struik’s transformation of this equation is obtained by writing 


uW re—m iw w A ir. 
and on | R (i.e. Z Re'?) equation (E) becomes 
dr ae . 
ke-37 sin 0 = 0, (F) 
do 
ga 
where k e3m, 
arc2 


His equation (42) should then read 


“d k,@ = kP(0,r)+26, (42’) 

“ago 

where z= k—k, s ” (43) 
n=1 

and P(6,7) = e-7* sin 6—8@. (44) 


Equation (42’) differs from Struik’s equation in the sign of the last term. 
To solve equation (42’) the following expansions are required: 


S n - liar 
Ww 2enFe > ao, = 6, tT p> 
us 1 


P(6,7) => P.n", kP(0,7) => wp". 
n=1 


n=1 
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The first-order solution is simply w, = —i({+¢-1), ky = cotha. Equating 
powers of u in equation (42’), Struik’s equation (54) takes the form 
| 
oT kn = Xnthn19, (nm = 2; 3,...), (54’) 
da 
n—2 
where Xn = w, + > £0, (55’) 


v=1 
Considering the form of the function P(@,7) and the assumed form of the 
functions w, (r = 1, 2,..., n—1) and k, (r = 1, 2,..., n—2), Struik obtains 
the polynomial for y,,, 


Xn = Inn Sin No+q, n-28in(n—2)o+...+-¢,,8in ea, (56) 


where e = 1 for n odd and e = 2 for n even, the q;;'s being constants. The 
unknown function w, must possess an expansion of the form 


a ox 
. , 
= 2 Yn» | aaa 2 Yny 8,5 
v=2 v 


and equation (57) should read 


x 
. ' , , 
6, = — > {—D,Ynvsin vo+ D, yn, cos vo}— D, vy, Cos a, 
v=2 
x 
- y aS a ae nm) 
tT = — D{S,y,,cosvo+S, y,, sin vo}—S, y,, sino, (57’) 
v=2 


where 7, = (’+¢-", 8, = —(”’, 8, = R°+ R”, and D, = R’— RR”. 

Substituting equations (57’) in (54’) and equating harmonics in o, we 
have the following results which contain differences in sign as compared 
with equations (60)—(63): 


tun = Yny = 9 forv > n, and for vy < n when n—v is odd; 
(—v8,-4 ky D,)yny — ia for 2 <v<cin; (60’) 
rat Dy kn»+ = 9 for n odd, k,-1 = 9 for n even. 


The required solution is thus given by 





n 
‘ ap q r 21! 
= ; Nv - 61 ) 
Wn Zz —vfS, ns he D, Ny ( 
1 k (c)sin o d ; , (c)sino d 
P a _ — ~ Ss ac= — @D,\a)}S o ao, 
ANC n-1 7D, ; Xn oO Co 7D, | n 
7 —7 (62’), (63’) 


where >’ denotes summation over 7,,, ),-25-++) Ne: 

A number of changes of sign are necessary in Struik’s section 10 following 
the changes in equations (42) and (61), but they do not affect the proof 
of convergence. 
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3. The amended solution to the third order 

To obtain the velocity of propagation and the form of the wave profile 
to the third order, it is necessary to calculate yo, y33, and k,. It is found 
that a correction is needed to the value of k, obtained by Struik, and some 
further corrections are needed in the subsequent calculations. The corrected 





value of k, is ve ’ 
= S- + € S. — 1 LS ’ 
Bit i eae F 11(9’) 
2 DR 
Struik obtains to the third order the correct results 
e-™ 1+p? 11(10) 
and 
° ° 2a j 47 S,-+- ° 
2=x+1y f — 26 sin J 4." 92 =2 f : sin oof _ 
, Cc CA “A D} cA 
27\? ,, 38,+458,14 . 6 i 
wifinge SP cee nf 1I(17) 
A D} cA 
where, however, the correct form of 6 is given by 
2a S, + 4 
| 2 3 Re 
X ‘—¢te 3 11(15’) 
The velocity of propagation is found from the equation 
A ' , , 
= emk = (14+ 3p") kotha”), 11(18’) 


which, after some simplification, gives 


2 _ VA D, - S+2%+4+12 , 
. D? 


2a S, 


The last result differs considerably from Struik’s result. It agrees with that 
obtained originally by Stokes (equation (33) of his paper), provided one 
obvious small correction to Stokes’s work is made, namely the insertion of 
the factor (277/A)? in the second-order term. 

On the free surface /, f = ¢+-7q, and separating II(17) into real and 
imaginary parts, Struik obtains Stokes’s parametric equations for the wave 
profile. When, however, ¢ was eliminated between these equations by 
Struik, he omitted some third-order terms and his explicit wave profile 
equation I1(30) should read 





7 27x 2Qr S,(S.4 _ 
y= qi 6? D, —bD, cos a | 2 1(S2+ 4) a dara _ 
ae | ig A A 2D, r 
_ (=) 43 9S 10S,- S,— 36 on aoe 3 So+ 65,41 5S8,+ 28 . ~_ 
A | 8D} r 8 D3 ws 


11(30’) 
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which differs from equation II(30) in the coefficient of the second harmonic 
and in both third-order terms. Struik’s expressions for the elevation above 
and depression below mean level and total-wave amplitude remain un- 
changed. 


If we now write the coefficient of cos(272/X) in I1(30’) in the form 


b, D bp, —?(27)" paps 

i) = —O —_ 2 5 +... 

1 1 1 :( x 1 

and assume b, D, >a as the depth of water tends to infinity, where a is 
finite, then the limiting form of I1(30’), after a suitable change of axes, is 


2rx 1 


a — 2 sia 4nx  3/27\? Pe 67x 
= acos “ COS es ‘ Ss , 
y » 7 2\) »  B\ar r 


which is essentially the same result as that obtained by Stokes for an 
infinite depth of water. 

The corresponding wave amplitudes for a particular depth of water may 
be determined from the relations 


A {, D,+-2D;+2D,—5D,)\ 
= — 2uD,+3y3 — ott Rece Mint 
” 2n | = Di | 
and h = , 4. = uD, 


where H is the height of the waves from crest to trough, and h is the mean 
depth of water. 


4. Note on the limiting cases of infinite depth and zero depth 

By an extension of Levi-Civita’s use of ‘majoring functions’ to the present 
case, Struik showed that a solution exists for a sufficiently small, but non- 
zero, range of values of ». Examining the equations satisfied by these 
functions, an estimate may be obtained for that range, and it will be seen 
that it tends to zero both in the case of very shallow water and also for an 
infinite depth. 

Using the following equations of Struik’s paper 


KuD, = 2KG(H) (106) 
and H = (14+ R*)Ru+ JKG(A)+ J KH, (108) 
where K = K—k, and G(H) = H?Y(H), (95), (111) 


it follows that K < J-1 since H > 0, and G(H) > 3H2. Since 
H > (1+ R?)Rp, 
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then G(H) 3(1+ R?)? R2u?, and using (106), 
Ku D, > 6(1+ R*)2K Rep? 
or K{D,—6R*y(1+- R?)?} > 6ky R2y(1 + R*)?, 
whence we obtain the sufficient condition p < py, say, where 
R?—1 
Mo 6 R3(1+ R2)2(1+dky) 

It follows that the solution is certainly convergent for values of u satisfy- 
ing this inequality. As the depth increases, R — oo and py > 0, and the 
corresponding range for » decreases. This appears at first sight to make the 
wave height vanish in the limiting case of infinite depth of water, but this 
is not, however, a necessary condition, and more refined inequalities can be 
found which lead to an upper bound for » which tends to zero less rapidly 
than py. This indicates the possible existence of the non-zero limit a 
referred to in section 3; it is known, of course, that a non-zero value of a 
does exist from Levi-Civita’s original work. 

It is perhaps worth noting that Levi-Civita’s investigation establishes the 
existence of permanent travelling waves only for quite small (though finite) 
amplitudes. In fact an examination of the inequalities used in his demon- 
stration restricts his existence proof to waves of height about 1/200 of the 
wave-length. Some refinements in his inequalities can be made which 
bring this ratio up to about 1/98. This falls far short of the height (about 1/7) 
estimated by Michell (8) and others for the sharp-crested waves, but there 
appears to be no physical criterion which would indicate a limit to the height 
below this value. 

The case when the depth of water decreases to zero is essentially different. 
In this case R > | and again p, > 0 and this now requires that the wave 


height shall tend to zero. 


5. The velocity of propagation: corrections to certain reports 

It was shown by Stokes that it is possible to define the velocity of 
propagation of permanent waves in a number of different ways, depending 
on the method adopted for solving the problem. It is not at once obvious 
that the various methods so far used in problems of this type do employ 
the same definition of velocity, but it will now be shown that, to the third 
order, these results are in agreement. 

The first method of solution is that used in the supplement to Stokes’s 
paper (4), where the expansion 

x ¢ =. 3 A, 


f 
_ ul 


elnm/exyp + ky e—(nm cds sin nm (1 ) 
C n 1 c 
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is employed, where m = 27/A and ¢ = —k on the canal bed. This method 
is essentially the same as that of Struik who employs the complex variable 
notation. The resulting expression for the velocity of propagation is given 
to the third order by 

y 99d | 9 9 2 
2 ga D, F , Sy +284 12 (2) ) 


e +." — . 2 
on St DS r) J “) 


the first harmonic in the wave-profile equation. 
The second method is an extension of that of Rayleigh, and uses the 
expansion d ‘ 
— as fonmy +h) > —RMY +A] - § 
et Dank! Uth) ¢ v+™\sin nm, (3) 
where / is the mean depth of the water. Again proceeding to the third 
order the velocity is given by 


aw Di{y_, Sat dl 4 (4) 
27 S;\ Di! A} J 
: , 2nah 
where D', = 2 sinh ‘ = +] iS — 2 cosh ~ ae 


To show that equations (2) and (4) are identical, we use the relation 
obtained by Struik to the third order, namely, 
k 2rk 
h ‘ 
; cA 
using the above definition of b. Substituting in (4), and using the notation 
of equation (2) we have 


2 ga D, 1. l ee »\° 14. 1 (27b\? hy 4 Sat 16 /27b\?) 

27 S, 2\AD, J |L | 2\A | Dt \Ass 
gr Dy{,_, Se+2S+12(27b)?) 
2a S,\ D} A | 
Exactly the same expression has thus been obtained by each method, 


and so to the third order, the same definition of velocity applies to both 
methods. Referring equation (3) to fixed axes, we have 


x 
¢ _— > a, {ermy +h) ¢ nmy sin nm(x—ct), 
Cc n=1 


h= (5) 


7 
+ — b? cot 
A 


and thus at each point of space occupied by the fluid referred to fixed axes, 
the mean horizontal velocity is zero, where ‘mean’ denotes the average with 
respect to time. This satisfies Stokes’s first definition of velocity and it may 
be measured from an ordinate through a fixed point on the canal bed. 
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hethod It should be noted that the expressions given in two reports of the U.S. 
wiable | Beach Erosion Board (2, 3) are incorrect, since equation (2) is used with 
'siven | the values of D, and S,, appropriate to equation (4). This is the same as 
omitting the second-order term in equation (5). 
Clearly, equation (4) is simpler and more convenient than (2) since it 
9 . 
\° involves the known h instead of k/c. 
— In conclusion I wish to thank Professor A. T. Price for help and advice 
lent ¢ . ° 
in the course of this work. 
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SUMMARY 
In this paper we obtain solutions in terms of complex potentials to static problems 
of generalized plane stress for infinite strips of isotropic elastic material. We solve 
the following types of problems: 
I. Specified stresses (including isolated forces) along both straight boundaries, 
(Sections (1)—(11).) 
II. Specified displacements along both straight boundaries. (Section (12).) 
III. Specified displacements along one straight boundary and specified stresses 
along the other. (Sections (13)—(15).) 
[V. Isolated forces in the interior of infinite strips with 
(a) both boundaries stress-free ; 
(b) both boundaries free from displacements ; 
(c) one boundary free from displacements and the other stress-free. (Sections 
(16)—(21).) 
PROBLEMS OF TYPE I 
1. Introduction 
FouRIER integral solutions to particular problems of infinite strips with 
specified boundary stresses have been given by Filon (1, 2) and Howland (3) 
in terms of Airy stress functions. Filon’s method has been modified by 
Hopkins (4) to solve problems of a strip with one boundary free from 
displacements. The present investigation gives rise to fewer difficulties 
with regard to convergence of integrals and applies to specified boundary 
displacements as well as to specified stresses. Further, it includes cases in 
which the boundary stress or displacement functions are not expressible as 
Fourier integrals and deals with isolated boundary forces without recourse 
to improper integrals. 


2. Notation and fundamental equations 

The following investigations are confined to elastic material in equilibrium 
and in a state of generalized plane stress. If there are no body forces, then 
it is known (see Stevenson, 5, 6) that the stresses and displacements may be 
expressed in terms of complex functions Q(z), w(z) of a single complex 


~ 


variable z (= 2+7iy) by means of the equations 


20 = 2(7x+97y) = O'(z)+0'(2), (1) 
—20 = —2(27—yy— — = 2Q"(z)+w’"(z), (2) 
8uD = 8u(ut+iv) = «Q(z)—2O'(Z)—a' (2). (3) 


[Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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The notation is that introduced by Stevenson in the papers cited, primes 
denoting derivatives, D being the complex displacement, a Lamé elastic 
constant, and x = 3—4o, where o is a modified Poisson’s ratio which is 
defined in terms of the Poisson ratio n by the equation (l1—o)(1+-) = 1. 
The complex potentials Q(z), w(z) are analytic functions of z and, therefore, 
possess derivatives of all finite orders at all interior points of the material. 
From (1), (2) we see that 

2(0—®) = 4(77+izy) = OQ'(z)+-O'(Z)+-2Q"(z)+o"(z). (4) 
We use the notation 


WY y 0 yy”, Wy Co yy", etc., 


where ¢, is a constant, and 
fr(u) | f(t)e-™ dt x,(u)+-ta,(u), (5) 


where «, and a, are real functions of the real variable u. By the Fourier 
integral theorem 
iy : ' 
f (x) (x, COS 7U—ay sin ru) du. (6) 
0 
Similarly, for functions 4(x), ys(a), x(x), we write 
by (u) €, +1€, Ur (u) 0, +10», X7(u) T+17¢. (7) 
We assume that the functions f (x), 4(x), o(x), x(#) satisfy conditions suffi- 
cient for the validity of the Fourier integral theorem and that a, etc., are 
bounded for all non-negative values of the parameter u. For brevity we 


shall write 


S = sinhyu, C= coshyu, A= uc, s=—sinhA, c=coshdA. (8) 


3. Specified distributions of stress along both straight boundaries 
We solve the problem of the strip subjected to the following boundary 
stresses: 


(yy +ixy)? = f (x)+i4(2) (9) 


(yy+tyx)} yb(x)+-2x (2x) 


ei ’ , : 
| (o, +77, )cos ru—(o_+-17,)sin xu] du, (10) 


where f, d, ws, y are real functions of x. In a previous paper (7) it has been 
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shown that the conditions in (9) may be satisfied by the sum of the two pairs 
of potentials 


Q,(2) — { T(z) dz, wo(2) — —2zQ,(z)+2 | Q(z) dz, (11) 


Q,(z) = | J(z)dz, a, (z) = —20,(2), (12) 


Sr(uje™” du, J(z) = — - [ hue" du. (13) 
* J 7 
For any real integrable functions 8,(w), By(w), y,(#), yo(w) of the real variable 
u we define the integrals 


where J(z) = 





x 


H(z) = - | [h(u)e*"+h(w)e-*"] du, where h(w) = B,(w)+i8,(u), (14) 
0 


x 


K(z) = : [ Lecupet-+-E(uye ‘eu) du, where k(u) = y,(w)+ityo(u). (15) 
"5 

From the results of paper (7) quoted above we see that, since H(z), K(z) 

have zero real and imaginary parts respectively when y = 0, the following 

pairs of potentials leave the line y = 0 stress-free: 


O(z) = | H(z) dz, w,(z) = —20,(z)+2 | Q(z) dz, (16) 


Q,(z) = [ K(z)dz, — wa(z) = —2,(2). (17) 


The suffixes 0, 1, 2, 3 will be attached to all quantities derived from Q,, Q,, 
Q,, Q, respectively. The solutions obtained below are, in the first instance, 
purely formal, differentiation with respect to z under integral signs being 
assumed valid whenever necessary. The validity of the solutions is con- 
sidered after their completion. Problems are solved for elastic material 
occupying the region R defined by —o <4 <0, 0 <y < ¢q, by using 
(11), (12) to obtain specified stresses along the boundary y = 0 and choosing 
h(u), k(u) so that (16), (17) give specified stresses along y = cy. From 
(4), (11), (12) we obtain 


x 


——— 1 f¢ - 
(yy+r12y)) = _ | e-¥"fl (1+-yu)a,+iyuas|eos 2u+ 
0 
+[iyua,—(1+yu)ag|sinazu} du, (18) 
anne ae 1 ~yus | , a 
(yy+izy), = 7h hate {[ywe,+(1—yu)ie, |eos xu 


0 


-+-[yue, —(1—yu)ie,|sin xu} du, (19) 
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From (1), (4), (16), (17) we obtain 


— — d : y */ Y Y 
(yy+i2Yy)o | {| B, Syu+iB,(¢ yu—S8)]cos xu + 

+-[B, Syu+iB,(Cyu—S)|sinxu} du, (20) 
ae l ' : ’ : 
(Fy +ixYy)s {| y, Syu+ ty,(Cyu S)]cos xu . 


+|y_ Syu+iy,(Cyu+S8)|sin vu} du, (21) 


>) . 
(7x+7Y)o = — | S(B, cos xu-+-B, sin xu) du, (22) 
0 
» fe 
(7x +YY)s | C(y, cos ru—+y, sin xu) du. (23) 


0 

Comparison of the integrand in (10) with the sum of the integrands given 
by (18)-(21) when y = ¢y yields the following equations for B,, By, y;, 2: 

By(s—Ac)—y, Ase ay (1-+A)+€,A]—9, = 9, 

B, As -Y1(s- Ac) ré | x, A—e,(1 = A)] rs = 0, 


YI (24) 
B,(s—Ac)+y2As+e | — x9(1-+-A)+, A]4 o, = 0, 
B, As+-yo(8+Ac)+efagA+¢,(1—A)]—7, = 9, 
whence 
B,(s?—A?) ¥9(A?— s?+-A-+-sc) —o,(Ac+-8) —A(Ae, +-87}), 
B,(s*—A?) x, (8s? —A?—A—sc)- a,(Ac+-s)— A(Ae,+- 873), (25) 


y,(s2—A?) €,(A2— s?-++-cs —A) +-7,(Ac— 8) +A(Aa,—80}), 
yo(82—A?) = e€,(s?—A?+A—sc)+7,(8—Ac)+A(Aa,—809). 

The above solution is valid provided that the J, J, H, K integrals exist 
and the formal expressions for I’, J’, H’, K’ are uniformly convergent with 
respect to z in the region R. In a previous paper (7) the author has shown 
that J and J integrals are, in general, satisfactory. The integrals for H’ 
and K’ obtained by formal differentiation are uniformly convergent in R 
if 


x 


[ eAu(|B,|, |Bel, ly1!, lye|) de < 00. (26) 
0 
From (25) we see that, at u 20, 
B,,B, = O(a, re, 7, Ae) + O(e-), (27) 
Vas Yq = O(a, de~>, 7, Ae) + Ole. (28) 


Assuming for the moment that all integrands are finite for all non-negative 
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values of uw, we see from (26)—(28) that sufficient conditions for validity of 
the solution are that, at u = 00, 
4) Og, €1,€, = O(1), 01, Fg, 71, T2 = O(u-%). (29) 
We will now illustrate the above theory by considering three particular 
cases. 
4. The infinite strip with specified, symmetrical, distributions of 
normal pressure along the straight boundaries 
In this case 
ay? _— yy) f (x), ft ( —2x) f (x), xy? ry} - 0. 
where f (x) is any given real function. In this case 
xy 01; Xs = Oy €1 €. = 7, = 72 0. 
If «,(0) is finite, then the H and K integrands are bounded at u = 0. (Com- 
pare Filon’s treatment (2).) The potentials for this problem are 
. ] x,(u) dr 
Q(z) = ¥ Q(z) = — [ elzu(eA— ]) 4 e-2u(e-A_])] M1 , (30) 
\ ) 2, ; 77. | Ne T A) 
0 
3 tf. 
w(z) = > w,(z) = —2zQ(z)— | [eu(eA+A—1)4 
0 TY 
; 0 
x,(U)e, dA 
}-e-teu(A+ 1—e-A)—4)] 0 (311) 
Provided that (29) is satisfied, the above potentials give a complete solution 
to the problem, boundedness of the integrand in (31) at « = 0 having been 
obtained by the introduction of a term independent of z. For Gaussian 
distributions of normal pressure given by 
” : ] r: 62 ~§2 42 4 39) 
f(x) = ~s5" Dries x,(u) = —e (v2 


the parameter 5 may have any non-zero value. As 5 —> 0 the boundary 
stress distributions tend to symmetrical, normal, isolated forces, but the 
solution cannot then be justified without a revision of the notion of integra- 
tion. 


5. The infinite strip with specified, symmetrical, distributions of 
shear along the straight boundaries 
In this case 
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The potentials in this case are 


1 ¢ a €.(u) dA a 

Q(z eizu(eA4 ])—¢e-izu(e-A4 1)] © 3: 
(2) = = | [ete+0) e+ cay’ (33) 
w(2) 2Q(z) 4 =Co | (eizu—¢ on) (34) 


0 
Ife,(u) = O(u) at u 0, the integrands in (33), (34) are bounded at u = 0 
and, provided that (29) is satisfied, the solution is valid. When this con- 
dition is not satisfied the solution requires modifications which are described 


later in this paper. 


6. Symmetrical, normal, isolated boundary forces 
In this section we assume potentials corresponding to isolated boundary 
forces +7Y at (0,0), (0,¢9) respectively and use the general method given 
above to obtain additional potentials which make the remainder of the 
boundaries stress-free. The potentials 
) » 
Q(z) — iY log z, w(z) = —-—iYzlogz (35) 
77 7 
correspond (5, 6) to an isolated force Y in the direction of the positive y-axis 
at a point O on the boundary of elastic material, O being the origin of 
rectangular cartesian axes O(x, y). These potentials correspond to a simple 
radial distribution of stress and, hence, leave all points of the real axis, other 
than the origin, stress-free. Along the line y = cy, however, they give the 
stress distributions 
9V 3 
2¥¢ f (x), ty 


r(x?-+-c?)? T(x 


(x). (36) 


IV c2 x 
- 2¥ ca 
UY os — 


Similarly, the potentials 
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2iY : 
Q(z) log(z—icy), 
‘ (37) 
20Y ; : 2c,Y . 
w(z) (z—1c,)log(z—ic,) — —*— log(z— icy) 


7 7 


correspond to an isolated force Y in the direction of the negative y-axis at 
the point « = 0 of the boundary y = ¢,, and along the line y = 0 give the 


stresses yy f (x), xy 7 d(x). (38) 
Thus, for the solution of the problem of an infinite strip subjected to 
symmetrical, normal, isolated boundary forces, we require the potentials 
given in (35), (37), together with potentials which have no singularities 
inside or on the boundary of R and which satisfy the boundary conditions 
WG =f), 2 = —3y' = dia). (39) 


yy 
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The required additional potentials are, clearly, the sum of those obtained 


in sections 4, 5. In the present case 

















m4 = Y(A+1)e*, eg = —Yde 
and 


a a | 
Q(z) = "| (cisu_e-tew-dygh_y_1) © 


A(s-+-A)’ 


77 
0 


w(z) = —2zQ(z)+ 


at } {etzul (A+-1)(e-A—1)—A(2A+ 1)eA] + (40) 
7 . rer nr ) 


0 


ir 
4 e—izu-Al () 4. 1)(e-A— 1)—)]+- 4X(A-+-1 ry __ 

e [( + 1)(é ) ] (A+ l)e (64 » 
The potentials in (40) give satisfactory solutions for stresses and displace- 
ments. 


7. Failure of the above general solution 

The formal solution obtained in section 3 to the problem of the elastic 
strip with specified stresses along the boundaries is unsatisfactory because 
the expressions for 8,, B., y,, y2 in (25) are, in general, infinite at w = 0 and 
the integrals defined in (14), (15) are divergent at their lower limits. 
Sneddon (8), who has recently given a formal solution to this problem in 
terms of Airy stress functions, has derived integrals of a similar type. We 
discuss only the cases for which 


B, = O(u-*), B, = O(u-), y, = O(u-"), ye = O(u) asu>O0. (41) 


Equations (41) are suggested by the following considerations. Consider the 
inequality 


XD 


| [tf (t)|, |t°(t)|, |t?H(2)|, |t?x(t)|] dt < oo. (42) 


We shall refer to the above as (42) when p = 0, as (42)* when p = 0, 1, 
and as (42)** when p = 0, 1, 2. If (42)* is satisfied, 


ys €y, Oy, T, = O(1), QM, €9, Og, T2 = O(u) as u> 0. (43) 


« 


From (25), (43) we can easily show that f,, etc., satisfy (41). In order to 
make the general solution satisfactory it is essential that we modify the 


integrands in (14), (15) without violating the boundary conditions along 
y = 9, Co, or introducing divergence of the integrals at w = oo. 


Wy =: 


H,(: 


Hy, 


The 
(45 
ten 
to | 
wil 


W 





btained 


(40 


Lr)” 


place- 


‘lastic 
cause 
0 and 
mits, 
mM in 


. We 


(41) 


r the 


(42) 








r to 
the 


yng 











GENERALIZED PLANE STRESS PROBLEMS 


§, Modification of the complex potentials 
Formal operations on the following pairs of ‘complex potentials’ give 
7/7 = ty = 0, for all functions g(u), m(u), n(u): 


H,(z) - [g(u)+g(u)] du, Q(z) = [ (2) dz, 


w(z) = —2Q(z)+ 2 | Q(z) dz, (44) 


H,(z) = - : zulm(u)—m(u)| du, Q(z) = [ Hate) dz, 


. 
0 


w(z) = —2Q(z)+2 [ 2¢2) ae, (45) 


. 


K, (2) : [ m(w)+- (wu) | du, Q(z) = [ Ke) dz, 


. 
0 


w(z) = —2Q(z). (46) 


The potentials in (44) correspond to rigid body displacements, those in 
(45) to tensions which are constant across the strip, and those in (46) to 
tensions which increase linearly from zero as we cross the strip from y = 0 
toy = ¢y. They may be subtracted from the potentials given by (14), (15) 
without affecting the stresses along the straight boundaries of the strip. 
We therefore consider the potentials 


H*(2) = H(z)—H,(2)—H,(z), — K*(z) = K(2)—K,(@), (47) 


where g(u) = m(u) = h(u), n(u) = k(u). We define potentials QF, w3, 
0*, w* in terms of H*, K* by equations similar to (16), (17). From (47) we 
obtain 


\(e="— 1 —izu)+(B, —iB,)(e-*"—1+izu)|du, (48) 


) 


H*(z) =~ | [(B,+%B, 


Kk *(z) (v1 |-Gyy)(e%*"— 1)+(y,—ty)(€ tzu 1)| du. (49) 


Assuming that (41) is satisfied, we can easily see that the integrands in 
(48), (49) are bounded at uw = 0, whilst, if (27), (28) are satisfied, the 
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additional terms cannot cause divergence of the integrals at uw = 00. Using 
the modified integrals, the four potentials Q), Q,, QF, QF become 


28 





9 a 
Oss) we — | (on én, er — 1), (50) 
mi | u 
0 
oP 
Ot x = | te, rier, (51) 
v u 
0 
QO*(z) = : (B, +72B.)(e%" — 1 —izu+22u?/2)— 
Pe 1 
‘ , ’ : ae du . 
— (B, —iB,)(e-*"— 1+-izu+-27u?/2)| “— (52) 
Or ] r ‘ “re : . sa ae ai 
Q3(z) = — | [(ya ttye)(e*"— 1—izu)—(y,—ty,)(e-#"— 1+ izu)] A (53) 


iy 


0 

We have added further rigid body displacement terms to ensure con- 
vergence at wu = 0. Determination of wy, wy} requires second integrations 
and the introduction of additional terms which are independent of z to 
retain convergence at u = 0. The process presents no difficulty but is not 
described here in detail for the sake of brevity, w(z) being required neither 
for the determination of stresses nor displacements. We have now obtained 
a complete solution to the problem of the infinite strip with specified stresses 
along the straight boundaries. This solution is, however, not in a form con- 
venient for the determination of conditions at infinity and we therefore 
consider alternative ways of modifying the potentials. 


9. Alternative method of modification of the potentials 
In order to proceed further we assume, in addition to the conditions 
imposed by (41), that the functions B,, B., y,, yz have the following forms: 


B, = = + e = (54) 
B, = _ r = + Boo; (55) 
n= + ¥1,05 (56) 


where A,, Ay, B,, B;, M,, N, are constants, and f, 9, etc., are bounded at 


u = 0; note that 8, contains no term B,/u? and y, contains no term M,/u. 
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Sufficient conditions for the validity of (54)-(57) are that the boundary 
stress functions satisfy (42)**. However, whereas (54)-(57) cover all 


physical cases, (42)** are not always satisfied. If we let 


g(u) = B,(u), m(u) = 0B,/u, n(u) = M,/u?, 
n (44)-(46), then the integrals thus obtained are convergent at their upper 
imits and we may use the divergence at their lower limits to cancel that 
f H(z), K(z). With the above interpretations of g(w), etc., we again con- 


sider the modified integrals H*, K* defined in (47). We obtain 


. 


t 


H*(z) [ (8, +-7B.)e*"+-(B,—tBy)e i2u_ 98 (u)+2B,z u?| du, (58) 
ts if : . , a . 
K*(z) [(yittye)e (yi: — tye" — 2M, /u?] du. (59) 


\ssuming (54)-(57), we can show that the integrands in (58), (59) are 
bounded at wu = 0. The potentials determined from the above modified 
integrals are more convenient than those given in (52), ete., for the deter- 
mination of conditions at infinity. If (55) had contained a term B,/u? or 
56) a term M,/u, then the required modifying integrals would have been 
divergent at w ©, and reversion to the previous method of modification 
given in (50), (51) are required in all 


: 


isnecessary. The potentials Q), Q, 
cases to give specified conditions along y = 0, whilst, using H*, K* in 


16), (17), the required conditions along y = cy, are given by 


, if ' , i : iB, 2*|du 
(5 (2) G +-7B,)(e' 1)—(B, —iB,)(e-*"— 1) — 28, tzu + — | : 
7. u |u 
(60) 


(61) 


u u 


2M, ae 
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The modified complex potentials of the general problem are 


Q(z) = Q(z) +(e) + OF) + Q3@) 


; ‘ (64) 
w(z) = —2Q(z)+2 | [Qo(z)+ QF (z)] dz} 
10. Fundamental integral formulae 
The following formulae are required in the subsequent analysis: 
(i) If (42) is satisfied, then 
lim wr = Q, (65) 
porto J" sin xu} 


(ii) If, in addition to satisfying (42), f(w) is of bounded variation in the 
neighbourhood of u = 0+, then 


).. 
sin are du +42 f (0+). (66) 


— 1—cos xu 
(111) du = Ana}. (67) 
‘s u ~ 
0 
11. Stresses at infinity 
Using (65) we can show that, at x +00, 
Ty = YYo = Wo = FX, Th, = XY, 0. (68) 


Referring to the potentials QF, QF derived from H*, K* respectively, the 
equations corresponding to (22), (23) are 


x 


») fe 
(%#+-9H). = = | [| S(8, cos 2u+B, sin xu)—yB,/u?] du, (69) 
0 
ee 
(7x-+7Y)s | [C(y, cos ru—y, sin xu) —M,/u?| du. (70) 
0 


Equations (20), (21) remain unchanged by the modifications in the integrals. 
Using (20), (21), (69), (70) and the formulae in (65)-(67), we obtain for the 
stresses at x = +00 respectively 


— - 


Vy 0, ry = +y(B,y+2mM,)/2, (71) 


La lim [—|2|(B,y+M,)+(A,y—M)]. 
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In order to proceed further in the general case we assume that (42)** is 
satisfied. With this assumption we obtain from (5) that 


va 


fir(u) v5, (u) +ia,(w) i | tf (t)e-™ dt, (72) 


where the primes denote differentiation with respect to uv. Similar expres- 
sions may be obtained for 47,(u), ¥7(u), x'7(v). Further, it follows that, 
it u , , , , , 
X15 €45 Fj, T1, Xa, €9, Fa, Toe O(1), a 
ian: (73) 
Xj; €1, Tj, T1, XQ, €g, FQ, Te O(u). 
From (25), (73) we can show that the conditions in the neighbourhood of 
u = 0 are as given in (54)—(57) and that the constants A,, B,, M,, N, are 
given by ‘i — : ? 
c3 A, 6[ a5,(0) —o4(0) |— 3eg/ €,(0)+-7,(0)], 
ce) B, 6| 0, (0) x, (9) ], 


c5 MU, 3[ x,(0)—o,(0)], 


aN; 3{ «,(0) —o,(0)]—e,[ 2€,(0)+-7,(0)]. 


( 


The symbol {X, Y, NV} will be used to denote a force wrench at the origin 
which has foree resultant X+-iY and anti-clockwise moment ‘, and the 
suffixes 0, 1 denote the force wrenches at the origin which are statically 
equivalent to the external stresses applied along the boundaries y = 0, ¢,y 


respectively. From (5), (7), (72), ete., we obtain expressions for these 


wrenches in terms of «,(0), ete. For example, 


x ( ry dx €,(0). 
We thus obtain 
x 7a f— (0), —a,(0), a,(0)}, (75) 
oR AF. 4 fr, (0), 0,(0), —oy(0)—¢y7,(0)}. (76) 


From (71), the force wrenches at the origin which are equivalent to the 


stresses at a % are given by 


lim (X 


¢ *) F 4 ») . r d T 7 »? 7 7 
f o+X,)/2, —(Yo+h))/2, [a(Yot+%,)—(No+,)]/2}. (77) 
Thus equilibrium of the strip is maintained by the wrenches shown in Fig. 1. 
It is evident that infinite stresses occur at infinity if and only if (Y,+Y,) + 0. 
Their existence is due to the fact that, if transverse balancing forces occur 
at infinity, then their moment about the origin is infinite, and, subject to 


. 


(42)*, the infinite couples required for equilibrium of the half-strips x > 0 
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« <0 can be obtained only by infinite values of 7%. In the diagram the 


strip is assumed to occupy 


’ the region between x 











finite. 
Y, 
A 
mopar > : 
U(¥+%) N, a 31(¥)+ MN) 
1(X,+X,) < <— }(X, 
ss Y Yy Y . és 
(N+M) 4(%4+Y¥) A MY+Y¥,) 3X4 
+N - X, 
Fic. | 


PROBLEMS OF TYPE II 


= +1, where / is large but 


12. The infinite strip with specified displacements along the 


straight boundaries 


The author has shown previously (7) that the condition 


may be satisfied by the su 


xr 








Q,(2) 4 e&"f(u) du, 
TK 
0 
Q, (2) ee e="d,(u) du, 
0 


Also, the potentials 
Q,(z) = H(z), 
K 


» 
Q,(z) = =— K(z), 


K 


leave the line y 


D® = f (x) +id(x) 


m of the two pairs of potentials 








w (2) 20,(2)+(1—«) | Q 
w,(2) 20, (z)+-(1+k) Q 
w,(z) = —20,(2)+(1—n) f Q,(z) dz, 
wa(2) = —2Q4(2)+(1+) | Q(z) d: 


(78) 


(82) 


- 0 free from displacements. The functions H, K are 
defined in (14), (15). As in the previous case, the suffixes 0, 1, 2, 3 will be 





attack 
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ittached to all quantities derived from Q5, Q,, Q,, Q; respectively. Using 


3) we obtain for the displacements along the line y = Cp, 


a7kD} | e-A{[ (x —A)a, +7Aa, |cos xu+[tAa,—(«—A) xy |sin xu} du, 

0 
xD} eM i(«+A)e, +Ae, ]eos eu+[Ae,—t(«+A)eg|sin xu} du, 

; (83) 
akD} {{ —iB, As+-B,(«s+-Ac) ]eos eu+|[B,(«s+-Ac)+- 7B, As|sin xu} du, 
7KD} | {| y, As+-ty.(Ac— xs) |cos ru [ iy, (Ac KS) —y2 As|sin xu} du. 


If we require 
D} ws(v)-ty(x) |(o, iT, )COS xU (02-4 it,)sin vu| du (84) 


we have, by comparison of the integrand in (84) with the sum of the inte- 


grands in (83) 3,(ks+Ac)+y,As+p = 0, 
3, As+-yo(Ac—Ks)+ 4 0, 2 
; ’ (85) 
3, (Ac-+- Ks) yoAs+r 0, 


> AS y,(Ac Ks)-+t 0. 


where » ole A)a,-+A€es] KO, 
4 whe hie, Mia (86) 
} ( (A—k)a,4 re, | -KO», 
t e* (K A)e€s i Axx, | KT, 
and, therefor 3 (K2s2—)2) r(Ac—xs)-+-qAs, 
Yo(K*8"— A*) q(Ac +-Ks)-++-rAs, - 
(37) 
3y(K 78" —A*) p(Ac— ks) tas, 
y4(K28? A?) t(Ac+-«s) pas. 
The integrands in (14), (15) are not always finite at wu = 0. They are 
satisfactory if 
yy O(1 Bo, Y2 O(u-!) as u> 0. (338) 
We assume that » 1 and consider the following cases: 
1) a4, Qo, € Eo, 01, Oo, Ti, To O(1), 
p,q,7,t = O(1), By, Bos yy Ye = O(u-!) as u> 0. (89) 
Hi) 04, €,, 0, 7 = O(1), lo, €9, Fg, T2 = Olu) as u> 0. (90) 
ili X45 , €y, Eg, Oy, Oo, Ti, To O(1), 


O(u), O(u) asu>d. (91) 
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When either (ii) or (iii) are satisfied, 


O(1), r,t = O(u), Bi. y 


P. 4 v1 = O(1), 


Be, Ye 
The conditions given in (89) are the most general of the three, those in 
(90) occur when the boundary functions satisfy (42), whilst those in (91) 


occur in important physical cases which we shall consider later in this paper, 


O(u-1), 


[t can be seen that the conditions in (88) are satisfied in cases (ii), (iii), but 
that in case (i) the H, K integrals are divergent at the lower limits. We 
therefore seek potentials which do not violate the conditions along the lines 
y 0, cy but restore convergence, when necessary, at w 0. From (3) it 
can be seen that, formally, the following pairs of potentials give zero 
displacements everywhere: 


x 








Aui . r 
Q(z) sa B,(u) du, w4(2) 2Q),(z)+(1—«) | 0),(z) dz, (92) 
0 : 
4 . 
(2-(z) P y,(u) du, w;(z) 2Q.(z)+-(1+k) | Q).(z) dz. (93) 
TK 
v0 


Using the above potentials we obtain modified potentials Q*, w*, O*, w*, 


> 


corresponding to Qs, ws, Qs, ws, respectively, as follows: 


E Qui i ‘ . : 
O*(z) = Q,(z)—Q,(z) | (8, iB.)(e"="—1) 
TK 
0 - 
(8, —1tB,)(e-""—1)] du, 
a (94) 
* ~()* (- <-- Q . iru 4 
Ws (Zz) 2425 (Z) (1 K) [(B, UDs)(¢ ] zu) 
TK 
0 
’ ; du 
-(B, —2B,)(e-="— 1 izu)| 
3 u 
2u F ; 
Q*(z) = Q,(z)—Q,(z) | Ls iyy)(e!="—1) 
TK 
0 
(v¥i—tya)(e“*"— 1) dau, 
; ' ni ; : (95) 
ws (Z) 203 (z)—(1-+k) | [iv ’yg)(e'="— 1—2zu) 
TTK 
0 
. du 
(v1 —ty2)(e-*" 1+-7zu) | , 
u 


The above integrals are convergent at their lower limits if (89) is satisfied. 
Evaluation of stresses as well as displacements from the above potentials 
requires the determination of Q”, w”; the integrals obtained by two formal 
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differentiations of (81), (82) or (94), (95) are uniformly convergent in R 


and the solution is satisfactory if 


e*u*( |B, 30|, l¥l> lve )du < @, (96) 
0 
which requires, from (86), (87), that 
B., Bes ¥1> Ve = O(e“ue-*), p,q, r,t = o(u-*) as wu> @. (97) 


Hence the potential integrals and their derivatives are satisfactory at their 


upper limits if 


bas Gay Qa; Hq = OLA), 01, Fg, T}, Te = O(u-*) as U> OO. (98) 


PROBLEMS OF TyPE III 


13. Specified displacements along one straight boundary and 
specified stresses along the other 
The potentials given in (79), (80) may be used to obtain specified displace- 
ments along the boundary y 0 and the potentials in (81), (82) or (94), (95) 


may then be chosen to give desired conditions of stress along the boundary 


y= Cy. From (79)—(82), (4) we obtain the following expressions for the 
stresses given by the above four pairs of potentials along y = Cy: 
~—\4 pu ; \ » “ ») 
jy--ixy)s ef xg u(K— 1—2A)+-ta, u(2A—k 1) |ecos xu 
TK 


x, U(K— 1— 2A)+ ia, u(K+ 1—2A)|sin ru} du, 


(yy+ixry)} eS e, U l—«—2A)+ te, u(1—xK 2A) |cos wu 
TK 
0 
e, u(1+«-+2A)+ ie, u(1—x«—2A)|sin xu} du, 
— A ) . 
(yy --ixy)s f {1 uB,(s—«s—2Ac)+21uB,(c+K«e- 2As) |cos zu- 
™ 7K . 
0 
18(«s—s+ 2Ac)+-iuB,(Ke+e-+ 2As)|sin ru} du, 
— 1 LL F . ) 
(yy CY )a | Uyo( 2ZAS—C—kC) iuy,(2Ac 1 K8) |cos vu ! 
TTK 


Uy (K8 —s—2Ac) + uy,(2As—c—xc)|sin xu} du. 


(99) 
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If we require the stresses along y = ¢y which are given in (10), we compare 

the sum of the integrands in (99) with that in (10) and obtain the equations 
B, u(s—Ks—2Ac)+-y, u(2As—c—xc)+p = 0, 


B,u(e+«e+ 2As)+-y, u(2Ac+s—xs)+q = 0, 


B, u(2Ac+-Ks—s)+-y, u(2As—c—x«c)+r = 0, me 
B, u(e+x«c+ 2As)+-y, u(ks—s—2dAc)+t = 0, 
where eAp = Ua («—1—2A)—ue, («+14 2A)—Kero, |p, 
eAg wx, (2A—K—1)+- we,(1—«— 2A) — Kez, |p, 
er = wa, (e—1—2A)+ weg(«+ 14 2d)+ Kero, /p, — 
et = ua,(«-+ 1—2A)+-we,(1—«—2A)+ Kez, /p, 
and therefore gB,u = p(2rc-+-8—xs)+-t(24s—c—xe), 
Jy2U p(e+x«e-+ 2As)+#(2Ac+x«s—s), 
(102) 


gB,u = q(2As—c—xc)+r(Ks—s—2Ac), 
Jy, u = q(s—Ks—2Ac)+7r(e+x«c-+ 2ds), 
where 
g 402+ («+ 1)2e?— (x«—1)2s? 4r24+ «2+ 1+ 2n(c?+-8?). (103) 
Using the above expressions for f,, etc., the required potentials are the sum 
of those given in (79)—(82). The H, K integrals are satisfactory at u = 0 
if (88) is satisfied, and this is certainly so if 
Nas €4, Oy, Ty O(1), Ne, €9, Og, Te O(u) as u— 0, (104) 
for in such cases 
PY O(1), rt O(u), B71 O(u), Bass ye O(1) as w+ 0. (105) 
However, we can show that (89) applies to this problem as well as to the last 
so that boundedness of «a, (0), etc., is insufficient to ensure that the potential 
integrals are satisfactory at their lower limits. The modified integrals given 
in (94), (95) are, however, satisfactory in this respect. The formal procedure 
employed in obtaining the above solution is justified if 


Ny, Uy, Ey, Eo, O(1), 01, %, 7}, Te o(u-*) as u—> 0, (106) 
although it is obvious that we may relax these conditions considerably. 
The solution simplifies considerably if we put w® = v° = 0, as was done 
by Hopkins (4). This is apparent in the first of the two special cases con- 
sidered below, viz the solution involving specified displacements along one 
boundary and specified stresses along the other. 
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14. The infinite strip with one boundary free from displacements 
and the other subjected to a Gaussian distribution of normal 
pressure 


In this case 


ae oe é 72/8? 
( . ~ 
= F = @, zy' = 0, yy = — =a (107) 
where dis a parameter. Hence 
Yy Xo €; Ey Ty To 0» 0, Oo; — eu? a (108) 
7 KO 
From (101) p :. g=r={= 0, 
b- 
ind, from (102), 
;, = @, guup, = x[2Ac+s(1—x)]eu"4, (109) 
V1 Q, JPLUy> «| 2As+-e(1+-«) |e Seu? /4 (110) 
Thus we have QD, = Q, = wa = w, = 0, 
2ui 2 - 
and (2(z) {2,(z) + 22,(z) | (Ae*"+ Be-“") du, (111) 
0 
where 
quua | (2A+-1)e4+ ne“ ]e-Ou'4, quuB [(2A l)e-A — xeAJe dPu*/4. 
also 
2u ; ; 
w(z) 15(z) + w(z) 2Q(z) 4 (Pe®"4+ Qe-@"+ R) du, (112) 
7 
where 
guu?P | (« 1)(2As+-ce+x«c)+(1—x)(2Ac+s8 Ks) |e SPui4 
guu*Y | («-+1)(2As+ce+x«c)+(«—1)(2Ac+<8 Ks) |e Put 
guu*R 2(1 L-#)*, 


15. The infinite strip with one boundary free from displacements 
and the other subjected to an isolated force. 

Neither the solution given by Hopkins nor that described above is rigorous 
when 6 = 0, i.e. when the specified stresses have a singularity corresponding 
to an isolated force. A satisfactory solution in this case may be obtained 
in the following manner. 

Potentials corresponding to an isolated force F = X+iY at the point 
icy on the boundary of elastic material are given by 


7 Q(z) 2 F log(z— ico), (113) 
mw(z) = [2F (z—icy)—2 Fic, |log(z—icy). 
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By the method given for the solution of problems of a half-plane (7) we find 
that the following potentials annul the displacements given by (113) along 
the boundary y = 0: 


mkQ(z) = 2F log(z-4 icy) +4ic, F/(z+iey), 


mko(z) = —mKzQ(z)+[2(F—1«2F)(z+icy)+4ic, F |log(z+-iey)+ (114) 


The latter potentials give the following stresses along y = cy: 
mKb(x) = mKay} “[Xe(1 Kk) + 2V e9(K?—1)]-4 
+- 4c? Q?(Xa— 2c, Y)+ 64c4 Q3(2c, Y— Xz), 
ments) =< wath £ [Yale 1) -+2X¢q(3-+x2)]— 


4c? Q?(Ya+-10c, X)+ 64c8 Q3(2c, X+-Ya), 
where Q = (x?+ 4c?)-1, 


Thus the sum of the above two pairs of potentials gives zero displacements 
along the line y = 0 but introduces the additional stresses given in (115) 
along y¥ = Cy. We therefore introduce further potentials which give zero 
displacements along y = 0, have no singularities inside or on the boundary 
of R, and give 

yy = p(x), zy = x(z). (116) 


Hence we require the potentials given in (94), (95) with f,, ete., given by 


(102) and x(x), x(#) as defined in (115). In the present case 


Xy Ne €} €y v, 


2xe%o, = Y(«?+1+44A+4)2), 


2Keé 2g, X (x2?—1+- 42), 
(117) 
2xerAz, X («2+ 1—4A+-4)?), 


2Ke Ar, Y(1—x?+4A?). 


The modified potentials given in (94), (95) are required, as (104) is not 
satisfied; since the transforms of the boundary stress functions satisfy 
(89) no further modifications are required. From (101) 


pp —K0}, eq =I. pr KO9, pt K 


4 


(118) 
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and, from (100), 


gPi © Y| («?+ 1+-4A+-4A?)(2Ac-+ (x?-+-4A2— 1)(2As—c—xc)], 


8—ks)-+ 


2 


ee: “ Y| («2+ 1+4A-+ 4A2)(2As+-e+x«c)+ («2+ 4A2—1)(2Ac+-«s—s)], 


gP2" X| («2+ 1—4A+ 4A?)(2As—c—xce)+(—x?+ 1— 4A?)(xs—s—2Ac)], 


dy,u > > > 9 9 

A8® X| («2+ 1—4A-+ 4A2)(s—«s—2Ac)+ (—x?+ 1—4A?)(e+«c-+ 2As) |. 
(119) 

The complete solution to the problem is given by the sum of the potentials 

given in (113), (114), (94), (95). 

terms than do those in the potentials obtained by treating the problem as 


Although the integrands contain more 


the limiting case of a Gaussian distribution, the above solution has the 
advantage of mathematical rigour, and, further, the integrals involved 
converge more rapidly. For example, the integrand in (111) is Of¢ A) at 
u = 0, uniformly in R, compared with O(e-*) in the present case. 

We now solve problems involving isolated forces in the interiors of infinite 
strips. Complex potentials are employed to determine the stresses and 
displacements given by an isolated force in the interior of an infinite plate, 
and the lines y = 0, ¢, are made either stress-free or free from displacements 
by the introduction of further potentials which have no singularities either 
inside or on the boundary of R and which annul the stresses or displacements 


given by the isolated force potentials along the lines y = 0, Cp. 


PROBLEMS OF TyPE IV 


16. Formulae connected with an isolated force in the interior of an 
infinite plate 


The potentials corresponding to an isolated force F X+iY at the 


point z ia in the interior of an infinite plate are 
Q(z vF log(z—b— ia), 
(120) 
w/(2) ve F(z—b—ia)log(z—b—ia)+-vF(b—ia)log(z—b—ia), 
where va(«k+1) 2. 
From (3), (4) we find that the stresses and displacements are given by 
$777) «—1)[X(x—b)—Y(y—a)]  4(y—a)*[| X(a—b)+Y(y—a)] 
(a—b)?+-(y—a)? (x b)?-+-(y a)?|? , 
10H (l—x)¥(a—b)—(x+3)X(y—a) 4(y—a)*[ Y¥(x—b)—X(y—a)] 
(x—b)?+-(y—a)? | (~7—b)?+ (y—a)? ; 


(121) 
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and 
Cpu —2X(y — 2} = b)(y—a) L (1—«)X —«X log| (w—b)?+ (y—a)®}, 
v (x—b)?+(y—a)? 


8uv =. 2Y (y—a)? +- 2X (a—b)(y—a) , . - . P 

. ——— a — (1+ — ol (~—b)2+-(y—a)? 

v (a—b)?+ (y—a)? ( x)Y¥ —K¥ log|(x—b)?+(y—a) |. 
(122) 


17. Image potentials suitable for removing the non-evanescent 

displacements given by the isolated force potentials 

The whole of the preceding theory is dependent on the validity of the 
Fourier integral expressions for the boundary stress and displacement 
functions. This weakness, which is common to all previous methods of 
treatment of the infinite strip, may be removed by the addition of elementary 
potentials which we shall refer to as ‘image potentials’. In the case of an 
isolated force in the interior of an infinite strip with stress-free boundaries 
the stresses which are to be annulled have representative functions which 
are expressible as Fourier integrals and the general method developed above 
needs no modification, but, if either boundary is to be free from displace- 
ments, we can see from (122) that the displacements to be annulled contain 
terms which are infinite at infinity. However, the author has shown 
previously (7) that the potentials 


Q)(z) vF log(z—b-+ia), 
= (123) 
w(z) 2Q(z)+( F —xF)(z—b- ia)log(z—b-+ ia) 


remove the finite and infinite terms in the displacements along y = 0 given 
by the isolated force. We see that these potentials also remove the non- 
evanescent displacements along y = cy. From (3), (4) the stresses and 
displacements given by (123) are 


4yy — (1—«)X(a—b)+-(14+«)¥(y+a)—2YVy | 
po (a—b)?+-(y+a)? 


. 4yf ¥(y+a)2+-X (a—b)(y+a)] 
[(w—b)?+-(y+a)?}? 
(124) 
404 = (Kk—1)¥(~—b)+(1+n)X(y+a)+2Xy , 
(x b)? {. (y r a)? 


_4yl—X(y-+-a)?+¥(e—b)(y+a)] 
[wb HytaPzP 
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and 
? 2y|_ X (ye "(a—b = » 
_ 21 X(y = ( = } —(1—x)X+-«X log] (a—b)?+(y+a)?], 
y (w—b)?+-(y+a)? 
Suv 24 y. 4 MP I "(4 > , ¢ 
id y[X (2 5 “ite 2 @)| L(1+«)¥+«} log| (a—b)? + (y+a)?]. 
l (x b)2 (y+a)° 


18. Isolated longitudinal force in the interior of an infinite strip 
with the straight boundaries stress-free 
We solve the problem of an isolated force F X+0i at the point 
z= b-+ia in the interior of the strip 0 < y <c, with the boundaries 
y = 0, Cy stress-free. The solution to a similar problem was first given by 
Howland (3). From (121) we see that the potentials additional to those 


given in (120) must satisfy the boundary conditions 


19° = =(1—x)(a—b) 4a2(x—b) 
vX (x—b)?+a? © [(x—b)?+a?}?’ 

4.ry/° (x+3)a 4a3 

vX — (x—b)?-+a? * [(a—b)? +a? (126) 
1yy} (1—x)(a—b) 4(c,—a)2(x—b) 

vX ~ (w—b)®+(co—a)® * [(@—b)?+ (cy—a)? 

1zy! (x+3)(¢y—a) 4(c,—a)? 


vX (~—b)*+-(cg—a)* [(x b)*? + (C9 a)*|? 


Thus we require the solution given in section 2 of this paper. With the 
notation given in (5), (7), (9), (10) we have 


fa, aX ¢ («— 1—2au)sin bu, 
4a, av .X e~*"(K— 1— 2au)cos bu, 
te, mvX e—""(« + 1— 2au)cos bu, 
fe, = mvXe-“"(« + 1— 2au)sin bu, id 
: ; (127) 
fo, = mv Xe ”!«—1—2u(cy—a) |sin bu, 
to, = mvXe—"o—| «— 1— 2u(cg—a) |cos bu, 
AT, mv X ¢ «+ 1—2u(ey—a) |cos bu, 
47, mv.X e~“Co—®)| « +- 1 — 2u(cg—a) |sin bu. 
Now 
X15 Xo, Ey, Ee O(ue-@"), 01; To; Ti, Te O(ue-“o-) as u > oo, (128) 
and, since 0 < a < ¢9, (29) is satisfied and all integrals are uniformly con- 


vergent with respect to z in the region R. Using the notation of (54)—(57) 
we can see from (25) that in this case 
A, = 3X(c,—2a)/c3, N, = X(2cy—3a)/c?, 


(129) 
B, = B, = M,= M, = 0. 
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The stresses at x = +0 are respectively 
Ty = ry = 0, TX = + X| 3Y(Cy— 2a) —C4(2¢)— 3a) | Che (130) 


each of which is statically equivalent to the force wrench 
f r /¢ 4 . « 
{—X/2, 0, Xa/2}. (131) 


Fig. 2 shows the external forces on the strip in the case 3a/2 < c 














0 < 2a. 
y 
4 eecec ccc ecccccccccccccscccces |cceee se edeabkitins 
1X < b 1 
ai < < a ™ 24 
Fic. 2 
The force resultants at x +00 may be modified by the superposition of 


simple tensions parallel to the x-axis. 


19. Isolated transverse force in the interior of an infinite strip 
with the straight boundaries stress-free 
If an isolated force F = 0+-iY is applied at the point z = b+-ia in the 
interior of an infinite strip 0 < y < cy, then we can see from (121) that the 
additional potentials required to make the straight boundaries stress-free 
must satisfy the boundary conditions 
477° (l—x)a 4a 
vY (~x—b)?+4 a [ (a b)?4 a*}?’ 
477/° («—1)(x—6) 4a?(x—b) 


vY (a—b)?+a2 | | (a—b)?- a?’ (132) 
477} (x —1)(cg—a) 4(c,—a)* 

vY (w—b)?+(cy—a)? © [(x—b)?+ (cg—a)?]?’ 

4ry/} (x—1)(x—b) 4(x—b)(c,—a)? 


vY (w—b)?+-(c,—a)? | [ (~—b)?+ (eg—a)? |?” 
We again use the solution given in section 2. With the notation of (5), (7), 
(9), (10) we have 

4a, —mvYe—*"(«+1-+ 2au)cos bu, 


4a, = mvYe-*"(«+1+ 2au)sin bu, 


4e, mvYe-*"(«— 1+ 2au)sin bu, 

4e, mvYe~*"(«—1-+ 2au)cos bu, _ 
4o, = mvYe-"o-®| «+ 1+ 2u(cy—a) |eos bu, ve 
4o, mY e~“eo—a) «-+-1-+4-2u(ey—a) |sin bu, 

47, -avY e~"o~| « — 1 +-2u(eg—a) sin bu, 


mY e~ “oY ¢— 1 +-2u(ey—a) |eos bu. 
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GENERALIZED 
The orders of magnitude of these functions at w = 0 are the same as in the 
previous case and are, therefore, satisfactory. With the notation of (25), 
(54)-(57) we have in the present case 


A, 6bY /c3 B, = 6Y/c, M, = —3Y/ce, 


; diag (134) 
N, = 36Y /c, B, = M, = 0. 
The stresses at x +-oo are respectively 
iy = 0, ry 1 3Y y(y—c,)/c3, Px lim [3Y(cg—2y)( |x| 6)/c3], 
. (135) 
which are statically equivalent to the wrenches 
lim {0, Y /2, Y(xa—b)/2! (136) 


> 


which are shown in Fig. 3. 








y 
“ ‘ a 
Y(l+t $Y (l—b) 

b 
a | 
b) LY 
Y Y 
Fic. 3 


The above solution reduces to that of Howland (3) when the constant 6 is 


equated to zero. 


20. Isolated force in the interior of an infinite strip with the 
straight boundaries free from displacements 
In this case we begin with the isolated force potentials given in (120) and 
the image potentials given in (123). From (122), (125) we see that, in order 
to free the straight boundaries from displacements, we require potentials 
which satisfy the following boundary conditions: 


= b) asin aX (x 5) a*y (137) 


fu? vy 


Auy®/y 
(a—b)*-+-a? oe (~a—b)' 
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We require the solution given in section 12 and, with the notation of 
(5), (7), (78), (84), we have 


4ue*"x, = var(X cos bu—Y sin bu), 


(140) 


4ue*"e, = van(—X sinbu—Y cos bu), 


(X 

4ue "xy = var(—X sin bu—Y cos bu), 
( 
_- 


4uee, van(—X cos bu-+-Y sin bw), 





duc, (cg — *|sinh an—aet X cos bu+-[2e,) sinh aw—ae*"|¥ sin bu, 
Vit 

4uc, eA ’ jac 2( cp. “)sinh an X sin bu+-[2ce)sinhau—ae*"|Y cos bu, 
oar u 

ate [ 2c) sinh au—ae*"|X sin bu jac 2(cg ae ;)sinh au | Y cos bu, 
u 

ae [ 2c) sinh au—ae®"|X cos bu 2(co r “)sinh au act Y sin bu. 
ViT ; 


(141) 

From (140), (141) we can see that (89) is satisfied at uw = 0, but that (90) is 
not. However, 

Xy—O, = aY (e*”—1)cos bu+ O(u) O(u), 

€9—T, = aX(e*"— 1)cos bu+ O(u) O(u) asu>o. 
Thus (91) is satisfied and the modifications developed above are not 
necessary in this case. Also, since (128) is satisfied, all integrals are 
uniformly convergent in the region R. Satisfactory additional potentials 
are therefore given by (79)-(82). 


21. Isolated force in the interior of an infinite strip with one 
boundary free from displacements and the other stress-free 
We again begin with the isolated force potentials and the image potentials. 

In order to free the boundary y = 0 from displacements and the boundary 

y = Cy from stresses we require potentials which satisfy the conditions given 

in (137) along y = 0 and, from (121), (124), give along y = c 


0 
Pon = «)[X (a—b) Y (ce )— a)| 4(C, -a)?[ X (~—b) + Y (cy a)) 
7 (~—b)?+ (cy—a)? [(~—b)?+ (Co— a)*|? 
(«—1)X(a—b)—(k+1)V(eg+a)+2c,Y  4¢9(cg+a)[ X(x~—b)+Y(c,+a)] 
T (w—b)?+-(e,+a)? : [ (~—b)?+ (cy +a)?}? 
on (x—1)Y(x —6) } (k+3)X (co -() 4 4(c, sist ok b) X (co: a)|_ 
(x—b)?+-(c¢,—a)? | (w—b)*?+-(cy+-a)? |? 
(«—1)Y¥(a—b)—(x+1)X(c,+a)— 2c, X __ 4e6(Co+ a)| Y(a—b)—X(c,+a )] 
: . (~x—b)?+-(c,+a)? ~ [(x—b)?-+(e,4 ae ; 


(142) 
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We require the solution given in section 13 and, with the notation of (5), (7), 
(10), (78), we see that a, ag, €, €, are given by (140), whilst o,, o5, 7,, 7, are 


given by 


2a, eh . — 
[awe*"+-(«—1—2A)sinh au]X sin bu 
ViT 
+-| —aue*™+-(1+«-+ 2A)sinhau]|Y cos bu, 
) \ 
200 € , , , - 
- [awe*” + («—1—2A)sinhau]X cos bu 
|awe*"—(1+«-+ 2A)sinhau]Y sin bu, 
27, ef . . 
i |awe*”+-(1+-«—2A)sinhau]|X cos bu 
ViT 
+-[aue*" + (1—x«—2A)sinhau]Y sin bu, 
a 
: | —aue*” —(1+-«—2A)sinh au |X sin bu 
ViT 


+-[aue*”—(«—1-+ 2A)sinhau]Y cos bu. 
(143) 
In this case (104) is not satisfied but, since (89) is satisfied, the modifications 
obtained above are necessary and sufficient to ensure that the integrals are 
satisfactory at wu = 0, whilst conditions at wu = oo are the same as in the 
previous case. The solution to the problem is given by the sum of the 
isolated force potentials, the image potentials, and the potentials given in 
(79), (80), (94), (95). 
The author wishes to take this opportunity of expressing his gratitude 
o A. C. Stevenson for the advice which he has freely given. 
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NEUTRAL HOLES IN PLANE SHEET—REINFORCED 
HOLES WHICH ARE ELASTICALLY EQUIVALENT TO 
THE UNCUT SHEET 
By E. H. MANSFIELD? (Royal Aircraft Establishment, Farnborough) 

[Received 24 January 1952] 


SUMMARY 

It is shown that, in a plane sheet under any particular loading system, acting in 
the plane of the sheet, certain reinforced holes may be made which do not alter th 
stress distribution in the main body of the sheet. These reinforced holes (hereafter 
called neutral holes) necessarily have exactly the same stiffness and at least thy 
same strength as the portion of the sheet that has been cut out. The weight of th: 
reinforcement is usually greater than the weight of the sheet that has been cut out, 
though there are cases where it is less. 

Mathematical formulae are developed to determine both the shape of a neutral 
hole and the variation along the hole boundary of the cross-sectional area of th 


reinforcement. 


List of symbols 


Ox, Oy Cartesian coordinate axes. 

t thickness of sheet. 

Oy, y;T, Stresses in the sheet. 

i) Airy stress function (see equation (1)). 

yb angle made by the tangent to boundary of hole with Ox. 
, load in reinforcing member bounding hole. 

A,, section area of reinforcing member bounding hole. 
A, section area of straight reinforcing member. 

E strain in reinforcing member bounding hole. 

E Young’s modulus. 

v Poisson’s ratio. 


1. General properties of a neutral hole 

CONSIDER first the sheet in the uncut state. The stresses are such that all 
elements of the sheet are in equilibrium which, for a typical element of the 
sheet, requires that 


CO, OT ry 0 
C2 cy 
CO OT x 

and v= 0, 
cy CX 


+ Acknowledgement is made to the Chief Scientist, Ministry of Supply, for permissio! 
to publish this paper. Crown copyright reserved. Reproduced by permission of thi 
Controller, H.M.S.O. 


[Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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Both these relations are automatically satisfied by introducing a stress 
function @ such that the stresses are to be derived from it by the equation 


oh 
Oy 9 
cy” 
9 
c 2h 
OC, 3 ( ] ) 
Car? 
a 
Try 





CXCY 
The complete state of stress in a sheet, therefore, can be described by the 


stress-function ¢ alone. This stress function is usually introduced as an 


iid to the determination of the stresses in a plate subjected to given 
boundary conditions. For example, for a plain sheet it can be shown that 
ssatisfies a particular equation (V4dé = 0) which, together with the boundary 
conditions, is sufficient to determine ¢, and hence the stresses. Here, 
however, it is assumed that the complete stress distribution is already 
known, but ¢ will still be used because it admits of great generality and 
because the properties of a neutral hole can be expressed simply in terms 
of d and its derivatives. The function ¢ itself assumes a new and special 


significance in determining the shape of a neutral hole. 
1.1. The shape of a neutral hole 
To fix ideas whilst considering the equilibrium of the reinforcing member 
ind the adjacent sheet Figs. 1 and 2 are given below. 
P sin w+ dP sinw) 
4 


A 
( —>P cosy +d (Pros y) 





K ~ Tx 
= A 
Pos y ~{ ™ry' 
G,; 
. -e 
P siny 
Fic. 1. Sheet wit! utral Fic. 2. Forees acting on the 


hole. element ABC. 


[t can be shown that, in a practical constfuction, the bending stiffness 
of the reinforcing member is negligible compared with the tensile stiffness 
so that the reinforcing member has the properties of a chain in that the 
line of action of the load P in the member is directed along the length of 
the member. For the element shown in Fig. 2 the conditions of equilibrium 
are therefore d(P sinus) t(o,,dx—r,,, dy) 


; (2) 
and di P cosa) t(r dx- a, dy) j 


‘ey 
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Equations (1) and (2) can be combined to give 


(1/t)d(P sings) — i dx + ht dy | 
ox? Cx 


xey 3 
oo ‘ 2 4 ) 
and (1/t)d(P cos) = — : ¢ dx —* Y ay 
Oxcy oy” 


These equations are in the form of total differentials and may therefore 
be integrated to give . / 

= - P sind Od | , 

‘ _ aoe 


t Cx (4 
P cos - ag 
t oy 


where a and 3b are arbitrary constants. 
The load P may be eliminated from equation (4) to give 





Od/ex+a 
tan — . 5 
ep/oy = 0) 
Substituting dy/dx for tanys, equation (5) may be integrated to give 
d+axr+byte = 0 (6) 


as the equation for determining the shape of a neutral hole. 
Terms of the type (aa+by+c) can be added to ¢ without altering the 
stresses and so there will be no loss of generality by writing equation (6) as 
d 0. (7) 
Hereafter equation (7) will be used instead of equation (6) as this appre- 
ciably simplifies the presentation of results. It is worth noting that equa- 
tion (7) is derived purely from considerations of staties and is independent 
of the elastic properties of the sheet. 


1.11. Hole bounded by arcs 

Since there are no restrictions on the constants a, b,c that may be included 
in ¢ it will be seen that there is a large variety of curves from which the 
hole shape may be chosen. Furthermore, the hole shape may be bounded 
by ares of curves, each of which is determined by a different set of values 
of a, b, c; in this case it will be necessary to apply balancing loads at the 
junction points of adjacent arcs to ensure equilibrium of the loads in the 
reinforcing members. It can be shown that such balancing loads can 
normally be produced by inserting a simple tension or compression member 
from one junction point to another. 

The calculation of such balancing loads is straightforward. Suppose 


are 1 is given by - 
' : y d a,x by 1e—Q 


and are 2 by d, = ¢' +a,x4+b,y+c, = 0; 
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then from equation (4) the vertical (or y) component of the load in ares 1 
and 2 is t(éd'/@a-+-a,) and t(é¢’/éa+-a,), respectively, and there are similar 
expressions for the a-components. It follows that the vertical and hori- 
zontal components of the balancing load are respectively 


t(a,—4@z) | 


. (8) 
and t(b,—b,) j 


1.2. Section area of the reinforcing member 
The section area of the reinforcing member can be determined from a 
knowledge of the load P and the strain ¢ in the reinforcing member by the 


relation A _ P/ Ee: 


m 


P is determined from equation (4) by eliminating %, giving 


rt . 


From elementary theory the strain « is related to the sheet stresses as 
follows: 
Ke cos*ys(o,,—vo,,)+sin*ys(o, —vo,)+ 2 sin ps cos #s(1+v)r,,. (10) 
Equations (9) and (10) can be combined with (1) and (7) to give 
(* fe 213 a 2 &(ed\? | 207d Od 0p 
}  Oy?\ey) — Oxy Ox Cy 
(e*h/ed\*  e*h/(ed\? 267d ed &p\]-} 
‘ 4 = Pl) (1d) 


oy Cy*\ ex exey Ox oy} 


9 
Cx OXL"\ CX 





oy 
\c se 

2. Examples 

2.1. General case of uniform stress distribution 

If the axes Ox, Oy are chosen parallel to the directions of the principal 
stresses f,, f, the stress function is given by 

d= 4f,y?+f.%7)+ar+by+e 

and so the neutral hole is in the form of a conic. Furthermore, if f, and f, 
are of the same sign the conic is an ellipse with lengths of axes in the ratio 
\ (fifo). The constants a, b, c merely determine the position and size of the 
hole. 


2.11. Equal stress in x- and y-directions 
Such a stress distribution occurs, for example, in a thin-walled sphere 
subjected to hydrostatic pressure. The most general form for ¢ is 
db 3 f(x?+y*)+axr+by+e, 
so that the shape of the neutral hole (46 = 0) is circular. The fact that the 
constants a, b, c are at present arbitrary means that the circle can be chosen 
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to have any radius and to be situated anywhere. If we take the centre of 
the circle at the origin and let it have a radius r we can substitute 
d = $f(a?+y?—r?) = 0 
in equations (11) and (9) to find 
rt 
Am: l "i (12) 

and P = fri 

These properties of a neutral circular hole could, of course, have been 
deduced simply from first principles, though it is interesting to note that 
the circular hole is the only possible hole; furthermore the circular hole, 
with this value for A,,,, will not be neutral for any other stress distribution. 

There is no need to limit the number of such holes to one, in fact there 
can be two or more overlapping, as mentioned in section 1.11, provided 
there is inserted a tension or compression member to ensure continuity of 
load P and compatibility of displacements. Two examples will demonstrate 
this. 


& om | 
t r trt ¢ 4 4 











\ Tot 
NG», ‘G . 2 aoe 
ins \ an —> An2= 
Amr fe~ Zh PN ae 
> A 
f ~ 
B 
— i se 
mw =) < - 





| | | | | 
+t tered 
Fic. 3. Two circular holes overlapping. 


In the first example, shown in Fig. 3, two unequal circular holes of radii 
r, and 7,, whose centres lie on the 2-axis, overlap, meeting at an angle 
(0,+6,). It follows from equation (8) that the horizontal components of P, 
and P, at the junction points A and B will cancel out, but the resultant 
vertical component (P, cos 6,+P, cos @,) will necessitate a tension member 
between A B with a section area determined by the relation 
4, — —' (r,c0s 0, + r,.cos 8 
A;=; (7, cos #, +- 1, cos 84). 
—y 
A different arrangement is shown in Fig. 4 (a); in this case the sign of the 
5 5 z 
balancing loads at A, B is reversed and so (because members with negative 
section areas are not practicable) the members AC and BD must be present. 
These will again have a section area of 
oS 


(7, cos 8, 


—yp 


1, COS 6.) 
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NEUTRAL 
and loads of magnitude f#(r, cos 6, + 7, cos 6,) must now be applied externally 


to them. 
Fig. 4(b) is a combination of the types represented in Figs. 3 and 4 (a). 
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Fic. 4. Holes bounded by circular ares. 


2.12. Stress in y-direction double that in x-direction 


Such a stress distribution occurs in a thin-walled cylinder subjected to 

















hydrostatic pressure. Suppose 4 t ¢ | 
C, 20, tp oe 
then ¢ is given by 
. - 
db 1 f(a?+ dy?)+ax+by+c, fr 
so that the shape for the neutral hole is an 
ellipse with major and minor axes in the ratio \ \ 
¥2:1. Ifa. b,c are chosen so that the centre of { { Fs { 


the ellipse is at the origin and the minor axis Fic. 5. Elliptical hole with 


. . a 
is 2r, then the equation for the neutral hole ance 3:1. 


be 9 9 2 2 
will be d © OE ales by” a 0. (13) 
Substituting equation (13) in (11) gives 

A V2(1-+-2?/r?)! 
m ( <f (14) 
rt ]—2v+ 322/r* 


2.13. Uniform stress in y-direction alone 
If o, = f, the most general form for ¢ is 
d = 3 fx?+axr+by+c, 
and so the neutral hole (or cut-out, since the hole is not closed) is any 
parabola of the type $x a2@—ry = 0, (15) 
where 7 is arbitrary and determines the size. 
This parabolic shape is well known in connexion with the design of 
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certain suspension bridges in which the weight per unit length of span, 
i.e. to,, is constant. 
The section area A,, is obtained by substituting equation (15) in 
equation (11); we find 
ani ute) A, _ (1+4y/r)! 


- =— a 6 
rt 2(4y r—v)' ae 





rT [r.7) 
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Fic. 6. Ideal parabolic cut-out. 


Pr ad 
\K L 
Y—A m given by 
equation (16) 
Origin of one 
of the parabolas_| 
\ | —A j= 2wt 
8 














1F 4 


Fic. 7. Cut-out bounded by parabolic arcs. 


This expression becomes negative over the range 0 < y < vr/4, from Q 
to @ (say) in Fig. 6, and so any practical design must utilize those parts 
of the parabola from @ to R. In fact, as A,, is very large in the immediate 
vicinity of Q it will also be inadvisable, from a weight point of view, to use 
those parts of the parabola near Q. The inset in Fig. 6 shows the position 
of @ to scale taking v = }. It will be seen that the ‘useless’ region is com- 
paratively small; it is due entirely to the Poisson contraction. 

The simplest symmetrical arrangement with positive A,, everywhere 
would be that shown in Fig. 7. 

At the junction point J the horizontal components of the P’s will cancel 
out, but to offset the resultant vertical component 2fwt, given by equation 
(8), it will be necessary to introduce the tension member as shown. If only 
vertical external forces can be applied at K, L the necessary horizontal 
components could be produced by inserting a compression member of 
section area rt/2v between K and L. 
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HOLES IN PLANE SHEET 


2.2. Equal bending stresses in two directions 
The most general form for ¢ is now 


dx x+y3+axr+by+c, 


so that the shape of a neutral hole is, in general, a cubic. 

Owing to the more complex form for ¢ the coefficients a, b, c no longer 
refer directly to the position and size of the hole. Here, however, the 
coefficients will be chosen specially so that the stress function may be 
factorized. The advantage of this is that it will be possible to find a ‘closed 
form’ for the shape of the neutral hole. We take 


dx 2+y—r(a+y) = (x+y)(a?—azy+y?—1’), (17) 


and choose the shape of the neutral hole from the second factor only. 
This factor represents an ellipse with major and minor axes in the ratio 
¥3:1 and inclined at 45° to the x-axis. 
The cross-section area of the reinforcement, obtained from equations 
(17) and (11), is given by 
t(5r?—3ay)! 


n q ¢ ’ (18) 
'  6{(4—v)r?—3ary} 


which remains practically constant at its mean value of 0-50rt. 


2.3. Stress distributions with circular neutral holes and con- 
stant A,, 

It was pointed out in section 2.11 that the circular hole with the value 
for A,, given by equation (12) would be neutral for one particular stress 
distribution and none other. But the circular shape itself is not confined 
to that one particular stress distribution. For example, confining attention 
to plain sheet, any function which satisfies 

V4d = 0, (19) 

gives rise to a possible stress distribution. If, therefore, we take 
d = (x?+y?—r?) F(x, y) (20) 
and choose F(x, y) so that equation (19) is satisfied, then a stress system 
will be formed for which the circle is a possible shape for a neutral hole. 
The section area of the reinforcing member will usually not be constant, 
but there is a set of possible functions of F(x, y) in which A,, is constant. 
This set is characterized by F(x,y) being homogeneous in x and y. It is 
more convenient to use polar coordinates in discussing this set. Thus, if ¢ 
is the radius of the circle and n is the order of homogeneity of F(x, y), 


? i i 20 | 2C S 9 9 P 
equation ( ) become d (r2—c?)r"f(8), (21) 
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and it is known that to satisfy equation (19) 
f(@) = asin né+b cos n8. 


Substituting in equation (11) it will be found that 


ct 


w (2n+1—p)’ 


A 


which is a constant; so that, except for the special case with n 0 con- 
sidered in section 2.11 the weight of reinforcement inserted is always 
appreciably less than the weight of material removed. The reason for this 
lies in the fact that the sheet which is to be removed is comparatively 
lightly stressed. 
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SUMMARY 
This paper is an extension of the method of mean strains introduced by Betti 
(1), many examples of which are given by Chree (2). The fundamental theorem 


of the method is proved more directly. It is used to determine the total end force 
on a cylindrical elastic body when the displacements on the plane ends are known 
and to solve the analogous problem in plane strain. 


1. Theorem 
In any elastic body of volume V bounded by a surface S, throughout which 


the components of strain are continuous, 





| td é‘ ] (n) ’ a iw (n) : 
= (U; n Uu; n;) dS 7 x; T’, dS _ pot Xp T;, dS — 
S 8 § 
lt+y ff ; ; re. £ R . 
a | px(f,—X;) dr+ oe px,(f,—X,) dr. (1) 
p V 


The quantities with suffixes are Cartesian tensors, the suffixes taking the 


values 1, 2, and 3. The notation is: 


u, displacement, 

nN, direction cosines of outward normal to element of surface dS, 
xX, coordinate, 

DP; stress tensor, 

é,, = strain tensor, 


T, = stress; force per unit area acting on a surface with normal n,, 


f acceleration, 

X, = body force per unit mass, 

é, Kronecker delta l ifs = 3, Oifi ~), 
v Poisson’s ratio, 

E Young’s modulus, 


p density. 


Proof. By Green’s lemma 
‘ . sy 
U;N; dS - u; dr. 
Cx, 


« . u 


S V 


{Quart. Journ. Mech. and Applied Math., Vol. VI, Pt. 3 (1953)] 
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Hence 


OX; OX; 


] ‘ " 1/eu; | eu; 
5 [ (ujme-turns) dS = | ( jy ) ar 


ry « 


= |e,,dr, by definition of strain, (2) 


“= | F a z Pind) dr, by the generalized Hooke’s law, 


at 


a 7 ad 
= E | Pij dr — Bot | Prk dr. (3) 
V 4 


(n) 
Again by Green’s lemma and using 7; = p;;7,;, 








F (n) r 
x; T, dS —= v; Pik Ni, dS — | z "% (2; Pix) dr 
ke iS iis 

a a, ‘ 

95, Dix t+ | Xj —— Pix At 

: K OX), 

Vv V 

: | pi; ar+ | px;(f,—X ,;) dr (4) 
Vv Vv 


since, by the equations of motion, 

C ’ . 

« tateX. = phe 

CXL). 
Substituting for | pi dz from (4) in (3) proves the theorem. Equating the 

;: 
right-hand sides of (1) and (2) gives the mean strain theorem. 
This theorem cannot be directly applied to general curvilinear orthogonal 

coordinates, since the components of the transformation matrix from the 


cartesian to the general tensor will in general depend upon the position of 


the element dS. However, if summation takes place over the indices i and 
j, the integrand of each integral is an invariant and the equation will be 
true for any system of curvilinear orthogonal coordinates.+ 

+ Chree proves the mean strain theorem from the reciprocal theorem 


| ej Pij AT = €4; Pi; at, 

j i 
where the accented and unaccented symbols refer to two different stress systems for a 
particular elastic body. This reciprocal theorem is true for any curvilinear orthogonal 
coordinate system since the integrands are invariants. To prove the mean strain theorem 
a particular pj; is chosen equal to 1 and the other pj; to 0. Such is a possible stress system 
for cartesian coordinates but it is not so in general. For example, in both polar and 
cylindrical coordinates the stress system defined by 70 1, other stress components zero, 
violates the equation of equilibrium in the 6-direction when neither body forces nor 
accelerations are present. 
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Equation (1) consists in fact of six equations, obtained by giving all 
possible values to i and j. In the examples to follow, body forces and 
accelerations are absent. Equation (1) then becomes 





, lt+v ff vf @ 
- | mm ds = | 7, ds—*. | x, T, dS, (5) 
3 " § “§ 
lf lt+y fp ® lt+v ff. @ 
V, . (Uy N3+Uz NM.) dS E | xT, dS ~ | x, T, dS, (6) 
"5 _ " § 


(3) | and four more equations obtained by permutation of the indices 1, 2, and 
3 in equations (5) and (6). 
The total force Q; over a part S of the boundary surface S is given by 
~ (n) 
Q; | 7, dS. (7) 
s 
2. Compression of a cylinder 
Consider a cylinder of any cross-section whose plane ends are normal to 
the generators. It is compressed between two smooth rigid surfaces, which 
are mirror images of one another in the central cross-section of the cylinder, 
(4) but otherwise are only limited in shape by having continuously turning 
tangent planes at all points and such that all normals to the surfaces 
make small angles with the generators. 
Let the central cross-section be taken as the plane x, = 0. If the height 


of the cylinder is 2a the end planes are 2, ta. When the ends are just 
the in contact with the rigid surfaces all over the ends, the displacement wu, at 
each point on the end is known from the geometry of the rigid surface.t 
On the end x, = a, n, 1; on the end 2, a,n, = —1; and on the 
na . ‘ 
t] curved surface of the cylinder, n, = 0. 
1e 1: ; : . . . . 
F Since the displacement wu, is odd in a, the value of the left-hand side of 
nO . 
(5) is known and is equal to 2 | u,dS. Ifno forces are applied to the sides 
und 
r a 
be Sod ; 7” 7 
of the cylinder, 7, = 0 on the sides. Since the surfaces are smooth and the 
(nm) (7) (7) 
ends are parallel to x, = 0, on the latter we have T, = T; = 0. T, is odd 
in x,. Consequently the right-hand side of (5) becomes 
i] 1 P (n (n) | 
2 aT, dS - aT, dS). 
. : oe E, 
nal 1=a 21=e 
= rt tatement is only approximately true, as the cylinder will expand laterally under 
em . ° - 
wit pressure. The value of the gral | w,dS will be underestimated by a fraction of the 
a 
nO order of Poisson’s ratio times the mean angle which the normal to the rigid surface over 
- its area of contact makes with the axis Ox,. The result obtained (equation (8)) is exactly 


hen the displacement wu, at the ends of the cylinder is explicitly given. 
















































382 D. R. BLAND 





(nm) 


Now the force F, on either end ( T, dS by (7). But, by (5), 


(1 ly ig (n) v re (n) | 
2 u, dS 4 ne | aT, dS —— | aT, dS 
E E 


| y 
. . “ 


ry=a ry a =a 
P <p E y 
or T, dS | u, dS. 
e a a 
rma ra 
= , 5 EK f¢ : 
lherefore the force F, on either end u, dS. (8) 
a 
re a 


[f the cylinder is compressed further, let the corresponding displace- 
ment be U. Then if v, is the displacement of any point on the surface, 


v; u,+U. By a similar analysis the force on the end is now 

Ee Re Er... 

Fk v; dS Uy dS + Uds 

a ‘i a J a A 
ry a ry a Ty a 
oe 
whence fk P+ — UA, (9) 
a 


where A is the area of the plane end. Equation (9) implies that once the 
surfaces are in complete contact with the cylinder, the additional force 
required for any further compression is equal to that force that would be 
required if the rigid surfaces were plane. 

Before the surfaces establish complete contact with the ends, equation (5) 
only gives an inequality. The integrand of { wu, dS is not known where 

rma 
the end is not in contact. It is less than the value necessary to bring it into 
contact with the plate. 

If known forces symmetrical about 2, = 0 are applied to the sides of the 
cylinder the above argument is unaltered except that these forces contribute 
to the second term on the right-hand side of (5). 

The final result is 

> ts) Ee i. tied 
T, dS u, dS (xv. T,+-4, T;) dS, (10) 
a a J 
m=a ra 
the last integral being taken over that part of the curved surface of the 
cylinder where x, > 0. 

As a particular example consider two rigid spheres each of radius R, 
compressing a right circular cylinder of radius R, and height 2a. The edges 
of the cylinder are just in contact with the spheres. 
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If r is the radius to any point on the surface of the cylinder, 


(r) R? rz 
u(r - a “ 
2R, 2R, 
Re , a 
so that u, aS | a a dr —he 


The compressive force on either sphere is, by (8), equal to 


_ hi R4 
7B RS (11) 
4ak, 


3. Applications to plain strain 
An equation (12) analogous to (1) exists in plane strain. Consider any 
cross-section of the body of area S whose normal lies in the direction of 
plane strain Ox,. Let the cross-section be bounded by the curve C. Then 
ms | 
u, dS 


by Green’s lemma in two dimensions; n; are the direction cosines of the 


i 


outward normal to C' and suffixes take the values 1 and 2 only. Hence, 


“l/eu; ou 
E all ‘) ds 
») . an 
: J 2\ex; Oxy 


e,,d8 








" ({l+yv v(l+y), Dai 
, 5,;D..} dS 
|g Po Pu Pury 


. . . . - | 
. (2; Pk ) dS | (5: Pik T x; 

‘ J Wy, J \ 

C ( S S 


c 





Pik dS 
OX, | 





a l »>f @ v(l+y). —-- : 
(uw. “u.nN-;) ds - : = 7 ds \ ~ Ms ry. T;. dS — 
2 : scala E _? D - sii as 
( c oy 
l+4 ; ,. wWi+yvr). f : ;. : 
u px ,(f,—X;) dS got px,(f;,—X,,.) dS. (12) 
S S 


limitation explained in the note on p. 381. 
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For i = 1, j = 1 and in the absence of body forces and accelerations, i.e, 
in equilibrium with no body forces, 


> 


] 
U,N, ds = - 


Y (n) | ae (n) (n) 
+ Vv Vv l+yv ; 
| tds" | (a T+ayT) de 
4 4 / 


C ¢ 
LF (n) (1+) - (n) 


=> J atids—"F | a Teds. (13) 
C C 

As an example of the use of (13), we consider the compression of a 
rectangle between two smooth rigid curves, which are symmetrical about 
the line x, = 0. The sides of the rectangle in contact with the curves are 
x, = +a. Assuming that contact just exists all along x, = La, u, is 
known on2, = +a from the shape of the rigid curves. Onz, = a,n, = +1; 
on x, = —a,n, = —1; and on the other two sides n, = 0. Now uw, is odd 


in x,; hence i r 
, | wr,ds=2 | uy, ds. 


Cc rm=a 
(n) (n) 
The compressive force F per unit width = T, ds; on C T,=0 
everywhere. Therefore, from (13), 


| —2(1—v?) 


2 | u, ds i 


r=a 


so that - (14) 


In particular for two right circular cylinders of equal radius R pressed 
symmetrically in plane strain into a strip of rectangular cross-section and 
breadth 2b, on xv, =a 


—U;(Xz_) 


2b 
3h 
2 Eb? 
and F = - ——— . (15) 
o( ] v?) Ra 
A further application of equation (13) is in the cold rolling of metals. 
It is used to calculate the compressive force acting on the strip in the 
regions where it is elastically compressed between the rollers (3). 
+ Subject to the limitation explained in the note on p. 381. 
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DISLOCATIONS AND PLASTIC FLOW IN CRYSTALS 


By A. H. coTTRELL. Illustrated. (International Series of Monographs on 
Physics.) 25s. net 


A dislocation is a certain kind of lattice defect which when present in crystals pro- 
foundly alters many of their properties, especially their plastic properties. Increas- 
ing confidence in the validity of the idea of dislocations and increasing respect 
for its power have led in recent years to a great proliferation of researches on 
dislocations. The fruits of this work are at present scattered widely throughout 
many scientific journals but when gathered together they display an impressive 
achievement. The purpose of this book is to bring together these various develop- 
ments in the theory of dislocations and crystal plasticity in a single co-ordinated 
account. 

The book is divided broadly into two main parts. The early chapters are con- 
cerned with elucidating the structures of dislocations and their geometrical and 
elastic properties, while in the later ones the results of this formal theory are 
applied to numerous problems of crystal physics, in particular those related to 
crystal growth, crystal boundaries, slip and twinning, mechanical hysteresis, de- 
formation of poly-crystals, precipitation hardening, yielding and strain ageing, 
work hardening, annealing, and creep. 
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